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ABSTRACT 


Tills  research  extends  an  observation  program  that  recorded  the  niglit 
airglow  from  a  Tucson  ground  station  using  an  imaging  spectrograph  known 
as  GLO.  GLO  was  designed  at  die  Lunar  and  Planetary  Laboratory  of  tlie 
University  of  Arizona  to  observe  auroral  and  airglow  emissions,  and  recorded 
midlatitude  airglow  data  near  equinox  during  Space  Shuttle  mission  STS-69 
in  September  1995. 

GLO  observations  from  the  shuttle  recorded  the  night  airglow  layer 
seen  edge  on  at  the  Eartli’s  limb.  These  observations  from  orbit  exliibit  a 
fundamentally  different  picture  of  the  niglit  airglow  compared  to  observations 
from  the  ground.  GLO  data  also  represent  tlie  first  simultaneous  optical 
measurements  of  airglow  emissions  over  tlie  spectral  range  from  1150  to 
9000  A,  showing  global  emission  variations  in  the  night  sky.  Intensity 
variations  are  not  correlated  among  emitting  species,  unplying  greater 
dynamism  and  more  complex  chemical  interactions  in  the  airglow  tlian 
previously  assiuned. 
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Altliough  other  researchers  have  described  observations  of  organized 
waves  and  tides  in  die  niglit  airglow,  tliese  prior  observations  are  sporadic  or 
averaged  over  long  time  periods.  The  niglit  sky  intensity  variations  recorded 
by  GLO  do  not  exhibit  any  obvious  relationship  to  atmospheric  tides.  They 
may  instead  be  the  result  of  a  chaotic  superposition  of  upward  and  downward 
vertical  motions.  This  upwelling  and  subsidence  may  cause  tlie  decoupling  of 
airglow  emissions  in  tlie  O2  Atmospheric  band,  the  OH  Meinel  band,  and  the 
atomic  oxygen  green  line  at  5577  A.  Emission  enliancements  with  maxima-to- 
minima  ratios  of  4  to  12,  depending  on  emitting  species,  have  been  observed 
in  die  GLO  data.  Emissions  in  die  O2  Atmospheric  band  system  and  die  01 
(5577  A)  green  line  show  a  greater  d3mamic  range  of  variation  than  die  OH 
Meinel  band  system.  The  chemistry  along  a  limited  line-of-sight  can  be 
explained  by  classical  airglow  chemistry,  but  only  over  a  limited  altitude 
range.  Dynamic  effects  in  die  80  to  100  km  region  are  sufficiently  chaotic  to 
present  mixed  results  when  inferring  chemical  processes  as  a  function  of 


altitude. 
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RESEARCH  OBJECTIVES 


The  Earth’s  mesosphere,  tliermosphere,  and  ionosphere  (MTI)  are 
coupled  by  a  complex  set  of  physical  and  chemical  processes.  Airglow 
emissions  are  a  key  signature  of  phenomena  in  the  MTI  region  (from  about  60 
km  up  to  several  hundred  km),  and  spectroscopic  measurement  of  these 
emissions  can  constrain  atmospheric  models.  Hydrodynamic  waves  have  been 
predicted  and  described  in  detail  for  the  MTI  region  and  have  been  detected 
by  ground-based  and  airborne  observations.  These  waves  include  thermal 
tidal  oscillations,  planetary  waves,  and  gravity  waves.  The  global  coverage 
provided  by  observations  from  orbit  may  help  validate  models  describing 
gravity  wave  saturation  and  dissipation  in  the  mesosphere.  Such  observations 
may  also  provide  evidence  of  dynamic  coupling  between  the  E  and  F  regions 
of  the  ionosphere,  predicted  by  modelers  but  not  yet  unequivocally  observed 
in  die  niglit  airglow.  From  GLO  observations,  useful  constraints  to  current 
chemistry  models  of  the  night  airglow  may  also  be  inferred. 
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Objectives  of  this  research  include  tire  identification  and  classification 
of  atmospheric  waves,  if  any,  in  the  niglit  airglow  and  the  evaluation  of 
airglow  chemistry  models  by  observing  variations  in  night  airglow  emissions. 
Wave  structure  will  be  inferred  directly  fi’om  temporal  and  spatial  intensity 
changes  in  the  night  airglow  layer  along  the  orbital  track  of  the  Shuttle. 
Chemistry  model  output  will  be  compared  to  spectral  signatures  recorded  by 
GLO  under  nearly  solar  minimum  conditions  for  tlie  O2  Atmospheric  (0,0) 
and  (0,1)  bands,  tlie  OH  Meinel  (6,2)  band,  and  die  01  (5577  A)  green  line, 
covering  the  altitude  range  fi-om  65  to  125  km  in  die  terrestrial  niglit  airglow. 

The  Jolmston-Broadfoot  night  airglow  model  (Jolmston  and  Broadfoot, 
1993)  will  be  used  to  predict  vertical  number  density  profiles  for  solar 
minimum  conditions.  Tliis  model  utilizes  vertical  profiles  of  atmospheric 
constituents  fi*om  the  1976  US  Standard  Atmosphere  and  seasonally  averaged 
upper  atmospheric  models  compiled  in  die  Air  Force  Geophysics 
Laboratory’s  Handbook  of  Geophysics  and  Space  Physics  (1985).  These 
profiles  are  consistent  with  those  widely  used  by  die  aeronomy  community 
(e.g,  Hedin’s  1987  MSIS  model).  Computer-assisted  spectral  analysis  using 
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GLO  observations  of  tlie  Earth's  liinb  will  provide  calibrated  data  for  emitting 
species  of  interest.  An  Abel  inversion  algorithm  for  cylindrically  symmetric 
geometry  (Hassani,  1991;  Dasch,  1992)  will  be  used  to  infer  volume  emission 
rate  profiles  fi’om  observed  integrated  line-of-sight  intensities  in  a  simple 
spherical  shell  model  of  die  atmosphere.  These  volume  emission  rate  profiles 
are  directly  proportional  to  vertical  density  profiles  for  die  emitting  species  in 
question.  Inferred  vertical  density  profiles  will  then  be  compared  to  die 
predictions  of  die  Jolinston-Broadfoot  model.  Tliis  procedure  will  help 
characterize  die  actual  spatial  and  temporal  behavior  of  atmospheric  vertical 
structure,  as  well  as  improve  die  utility  of  die  models. 

Today’s  niglit  airglow  models  have  been  developed  primarily  fi’om 
groimd-based  observations,  wliich  have  good  temporal  resolution  but  poor 
spatial  resolution  in  die  vertical.  Top  side  observations  fi'om  orbit  have  been 
sparse.  Comparison  of  model  output  with  profiles  inferred  from  GLO  data 
will  either  validate  die  model  physics  or  suggest  improvements  to  the 
dieoretical  descriptions.  Airglow  models  include  complex  interrelated 
subroutines  wliich  calculate,  for  example,  excitation  rates  for  atomic  and 
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molecular  species,  excited  state  production  and  transport,  difiusion  and 
temperature  effects,  number  densities  of  excited  species,  and  volume  emission 
rates.  The  Johnston-Broadfoot  model  was  developed  at  the  University  of 
Arizona  and  is  available  at  LPL  West.  It  is  representative  of  most  of  die 
models  in  use  today  (Garcia  and  Solomon,  1983;  Allen  et  al.,  1984;  Solomon 
et  al,  1988;  Rodrigo  et  al.  1986,  1991;  Hedin,  1987,  1991;  McDade,  1991; 
Mlynczak  et  al.,  1993).  Such  models  can  only  be  refined  and  vahdated  by 
considering  data  with  die  widest  possible  temporal,  spatial,  and  spectral 


coverage. 
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1.  INTRODUCTION 


Until  recently,  models  of  the  Earth’s  mesosphere,  tliermosphere,  and 
ionosphere  (MTI)  have  depended  on  observing  airglow  and  auroral 
emissions,  primarily  using  ground-based  remote  sensing,  short-duration 
rocket  fliglits,  and  a  limited  number  of  satellites.  Such  observations  are 
hmited  in  temporal,  spatial,  and  spectral  coverage.  To  address  tliis  difficulty, 
die  National  Science  Foundation  has  created  a  program  known  as  tlie 
Coupling,  Energetics,  and  Dynamics  of  Atmospheric  Regions  (CEDAR). 
CEDAR  is  a  key  part  of  initiatives  aimed  at  creating  a  global  observation  and 
modeling  capability  for  the  entire  depth  of  die  Earth’s  atmosphere.  However, 
until  permanent  orbital  platforms  exist  (such  as  die  International  Space 
Station)  dedicated  to  such  continuous  global  coverage,  data  collected  on 
Space  Shuttle  fliglits  can  help  to  fill  the  gaps  in  our  knowledge  and  improve 
our  models. 

Researchers  at  die  University  of  Arizona’s  Limar  and  Planetary 
Laboratory  West  (LPL  West)  have  designed  an  imaging  spectrograph  diat 
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usefiil  for  studying  die  MTI  region  from  orbit.  Tlie  Arizona  Airglow 
Instrument  (GLO)  is  a  modified  version  of  die  Arizona  Imager  Spectrograph 
(AIS)  previously  built  under  contract  to  the  U.S.  Air  Force  (Broadfoot  et  al, 
1992;  Viereck  et  al,  1990).  GLO  has  already  flown  on  several  Shuttle 
missions  and  additional  flights  are  planned  over  the  next  few  years.  Some  of 
die  future  missions  will  be  augmented  by  coordinated  ground-based 
observations  under  die  CEDAR  program.  GLO  data  can  be  accessed  by  die 
CEDAR  and  space  science  communities  via  die  Internet. 
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1.1  Historical  Background 

The  history  of  aurora  and  airgjow  physics  has  been  summarized  by 
Chapman  (1967).  Tire  aurora  has  been  known  since  antiquity,  but  a  scientific 
description  distinguisliing  aurora  fi-om  airglow  is  more  recent.  Angstrom  was 
the  first  to  study  the  aurora’s  spectrum;  in  1867,  he  found  a  prominent 
yellow-green  line  with  a  simple  spectroscope,  but  did  not  accurately 
determine  the  wavelength.  In  1901,  Newcomb  observed  tlie  faint  non-auroral 
light  of  die  night  sky  (tlie  night  airglow)  and  hypothesized  the  source  to  be 
multitudes  of  stars  too  dim  to  be  seen.  In  1913,  Vegard  identified  bands  of 
molecular  nitrogen  in  auroral  spectra  diat  were  also  found  to  be  characteristic 
of  airglow  spectra.  Tlie  faint  light  of  die  niglit  sky  identified  by  Newcomb 
was  studied  further  by  Slipher,  beginning  in  1915  and  for  many  years 
diereafter,  at  die  Lowell  Observatory  in  Flagstaff,  Arizona.  Babcock 
detennmed  die  5577  A  wavelengdi  of  Angstrom’s  prominent  auroral  yellow- 
green  line  in  1923.  Using  a  Fabry-Perot  interferometer,  Babcock  unknowingly 
observed  the  light  of  die  airglow  and  not  the  aurora,  as  he  had  intended!  In 
1924,  McLeiman  and  Shrum  proved  that  die  5577  A  line  was  caused  by  an 
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electronic  transition  between  tlie  metastable  *D2  and  ^So  levels  of  atomic 
oxygen.  Sliplier  identified  the  red  lines  of  atomic  oxygen  at  6300  A  and  6364 
A  in  airglow  spectra,  and  also  noted  that  tlie  5577  A  line  was  always  present 
in  such  spectra.  (The  6300  A  and  6364  A  lines  are  also  prominent  in  auroral 
spectra  and  correspond  to  tlie  ^D2  to  ^Pi  and  ^D2  to  ^Pi  forbidden 
transitions,  respectively).  Observations  of  tliese  tliree  lines  in  botli  the  aurora 
and  airglow  created  some  confusion  initially.  Researchers  eventually 
distinguished  airglow  from  auroral  emissions  using  the  fact  tliat  excitation 
potentials  exliibited  by  airglow  spectra  were  much  smaller  than  those  in 
auroral  spectra.  They  estabhshed  tliat  tlie  aurora  was  the  result  of  more 
energetic  processes  than  tlie  airglow,  but  was  otlierwise  very  similar. 

In  1931,  die  fourth  Lord  Rayleigh  described  tlie  airglow  as  the  “non¬ 
polar  aurora.”  Tliis  tenninology  was  gradually  superseded  by  “airglow,”  first 
introduced  by  Elvey  in  the  1960's  following  a  suggestion  by  Struve.  The 
absolute  intensity  or  surface  briglitness  of  tlie  airglow  5577  A  line,  expressed 
in  millions  of  photons  per  second  per  colmnn  cross-sectional  area  (in  cm  ) 
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along  the  line-of-siglit,  was  first  measured  by  Rayleigh  in  1930.  In  his  honor, 
the  rayleigli  (symbol  R)  has  been  adopted  as  tlie  unit  of  sky  briglitness. 

The  airglow  in  the  sunlit  sky  is  referred  to  as  the  dayglow.  In  1932, 
Dauvilher  proposed  that  dayglow  emissions  are  caused  by  a  continuous  flux 
of  higli-energy  solar  electrons  (10^®  eV)  interacting  with  the  Eartli’s 
atmosphere.  Chapman’s  earlier  tlieory  of  1931,  generally  accepted  today, 
states  that  airglow  photons  are  emitted  during  recombination  processes 
among  electrons,  ions,  and  dissociated  atoms  and  molecules  created  by  solar 
UV  flux  during  the  day.  Airglows  have  also  been  observed  in  tlie  atmospheres 
of  otlier  planets,  including  Mercury,  Venus,  Mars,  and  Jupiter,  and  perhaps 
recently  even  in  tlie  atnospheres  of  Jupiter’s  satellites. 

In  many  cases,  differentiating  between  die  aurora  and  die  dayglow  is 
not  a  simple  matter.  A  good  but  imperfect  working  definition  (O’Brien  et  al, 
1960)  is:  “If  excitation  results  from  direct  bombardment  by  energetic  charged 
particles,  it  is  aurora.  Otherwise,  it  is  airglow.”  (The  italics  are  mine).  Thus, 
the  aurora  is  primarily  caused  by  energetic  particles,  while  die  dayglow  is 
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mainly  caused  by  the  deposition  of  solar  energy.  In  the  dayglow, 
photoelectrons  can  compete  with  photoionization  and  UV  photons,  depending 
on  die  atnospheric  constituent  and  die  characteristic  energies  in  question.  On 
die  odier  hand,  the  night  airgjow  results  from  the  chemical  relaxation  of  die 
dayglow  atmosphere  when  the  solar  influx  is  removed.  Ultimately,  die 
production  of  atomic  oxygen  by  die  photodissociation  of  molecular  oxygen 
during  the  day  is  critical  to  creating  die  niglit  airglow. 
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1.2.  Da3^inie  Chemistry 


Hie  following  review  is  based  primarily  on  the  work  of  Tolimatsu 
(1990)  and  gives  tlie  background  for  the  discussion  of  niglit  airglow  chemistry 
in  the  next  section,  as  well  as  the  Jolinston-Broadfoot  model  outlined  in 
Appendix  F. 

The  photodissociation  of  molecular  oxygen  and  ozone  by  daytime  solar 
flux  produces  atomic  oxygen,  which  persists  after  sunset  and  dominates  the 
primary  reactions  in  night  airglow  chemistry.  Molecular  oxygen  is  dissociated 
during  the  day  by  tlie  absorption  of  solar  radiation  in  the  Herzberg  continumn, 
die  Schumann-Runge  bands,  and  die  Schumann-Runge  continuum.  As  a  result 
of  this  absorption,  die  following  photochemical  reactions  ensue,  creating 
various  states  of  atomic  oxygen  (Tolimatsu,  1990;  Chamberlain  and  Hunten, 


1987): 
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•  in  tlie  Herzberg  continuum: 

O2  (X^Sg-)  +  hv  (k  <  2424  A)  ^  O  (¥1)  +  O  (¥1); 

•  in  the  Schumann-Runge  bands: 

O2  (X’Xg-)  +  hv  (1750  A  <X<  2060  A)  ->  O2  V)  ^  O  (’Pi)  +  O  (’Pi) 
for  v'>  2  [V  are  the  vibrational  levels]; 

•  in  the  Schumann-Runge  continuum: 

O2  (X^Sg-)  +  hv  (>v  <  1750  A)  ^  O  (^Pi)  +  O  (^D2). 


Weak  absorption  occurs  in  tlie  Herzberg  bands  for  X  >2424  A, 
corresponding  to  the  forbidden  transition  O2  (X^Zg-)  O2  (A^Zu+)-  The 
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absorption  cross  section  of  O2  in  these  bands  is  typically  less  than  10"  cm  . 
Since  tliis  absorption  is  not  accompanied  by  dissociation,  it  is  neglected. 

Absorption  in  the  Herzberg  continuum  is  also  weak.  However,  with  O2 
abundance  in  our  atmosphere  being  so  liigli  and  diis  absorption  occurring 
where  tlie  solar  flux  is  large,  the  Herzberg  continuum  is  tlie  main  source  of 
atomic  oxygen  below  65  km  altitude,  where  it  is  important  in  the  formation  of 
tlie  ozone  layer  (Tolimatsu,  1990). 
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In  tlie  Schumann-Runge  bands  (1750  A<A,<2060  A),  dissociation  into 
atomic  oxygen  occurs  for  O2  vibrational  levels  v'  greater  than  or  equal  to  2. 
This  process  is  a  prime  source  of  atomic  oxygen  in  tlie  mesosphere  at 
altitudes  between  60  and  90  km. 

The  lower  thermosphere  absorbs  ahnost  all  incoming  solar  flux 
between  1350  and  1750  A  in  the  Schumann-Runge  continuum,  at  altitudes 
between  90  and  130  km.  This  absorption  is  the  major  source  of  atomic 
oxygen  in  the  ionosphere.  The  atomic  oxygen  so  produced  is  in  an  excited 
electronic  state  and  plays  a  very  important  part  in  the  chemistry  of  the  neutral 
atmosphere  at  these  altitudes.  Less  important  is  absorption  by  ozone  in  die 
Hartley  continuum  (A,  <  3100  A),  producing  excited  states  of  atomic  and 
molecular  oxygen  which  participate  in  the  chemistry  of  the  mesosphere. 
Additionally,  absorption  of  solar  Lyman-a  flux  at  1216  A  causes  about  ten 
percent  of  the  total  O2  dissociation  between  70  and  90  km  (Tohmatsu,  1990). 
The  net  effect  of  die  processes  described  here  is  to  estabhsh  an  “atmospheric 
cutoff”  near  3000  A.  All  solar  flux  shortward  of  diis  wavelength  is  ultimately 
absorbed  by  die  atmosphere. 
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1 .3  Niglit  Airglow  Chemistry 

Night  airglow  emissions  at  about  60  to  120  km  altitude  depend  on  the 
nighttime  vertical  distribution  of  chemically  active  atomic  and  molecular 
species  estabhshed  during  the  day.  Atomic  oxygen  dominates  the  chemical 
reactions  involving  NO,  NO2,  O2,  and  OH  that  create  the  features  of  the  night 
airglow  spectrum,  but  die  variabihty  of  atomic  oxygen  in  the  thermosphere 
has  yet  to  be  measured  systematically.  Ground  state  atomic  oxygen,  O  (^Pi), 
is  generally  believed  to  be  die  prime  agent  in  creating  die  niglit  airglow 
dirough  its  participation  in  tliree-body  association  reactions  forming  excited 
O2  states.  Atomic  oxygen  is  created  at  altitude  by  photodissociation  during 
the  day  and  diffuses  downward  into  die  thermosphere  and  mesosphere,  where 
it  takes  part  in  reactions  with  odier  species,  such  as  OH  and  O2.  Tlie  chemical 
reactions  coupled  in  the  niglit  airglow  are  schematically  depicted  in  Figure  1.1 
on  page  30  (Jolmston  and  Broadfoot,  1993). 
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Figure  1.1.  The  chemistry  of  the  night  airglow. 
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Tlie  two-part  Barth  process  creates  tlie  excited  states  of  O  and  O2  in 
the  night  airglow  (Chamberlain  and  Himten,  1987;  Rees,  1989). 

Step  One  of  the  Barth  process  is: 

O  (^Pi)  +  O  (^Pi)  +  M  ^  O2  (c  ^Zu-,  A'  ^Au,  A  ^Zu+,  b  ^Zg+,  a  +  M, 
wliich  is  followed  by  Step  Two: 

O2  (c  +  O  (’P,)  ^  O2  (X  ’Zg.)  +  O  ('So). 

(Some  losses  of  O  (’Pi)  to  the  Oj  and  O2  ground  state  occur,  see  Figure  l.I). 

Step  One  creates  the  O2  excited  states,  which  then  emit  as  follows: 

O2 (A ^Zu+)  ->  O2 (X  ^2g-)  +  hv  (Herzberg 1, 3000  A<X<  5000  A); 

O2 (c  ^Su-)  ->  O2 (X  ^Zg-)  +  hv  (Herzberg H,  3000  A<X<  5000  A); 

O2 (A'  ^Au)  ->  O2 (X  ^Zg-)  +  hv  (Herzberg HI,  3000  A<X<  5000  A); 

O2  (b  ^Zg+)  ->  O2  (X  ^Zg-)  +  hv  (Atmospheric,  X  >  7600  A); 

O2  (a  ‘Ag)  ->  O2  (X  ^Zg-)  +  hv  (IR  Atmospheric,  ^  >  1  pm); 

O2 (A'  ^Au)  O2 (a  *Ag)  +  hv  (Chamberlain,  3300  A<X< 4400  A). 
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Step  Two  is  quickly  followed  by  the  reactions 
0(»So)  O  (‘Di)  +  hv  (5577  A)  and  O(^So)  O  (^i)  +  hv  (2972  A). 

Tlie  metastable  transition  O  (*D2)  O  (^Pi)  +  hv  (6300  A)  is  possible, 
but  strongly  quenched  at  altitudes  below  200  km  due  to  its  long  radiative 
lifetime  of  110  seconds  (Rees,  1989).  The  6300  A  emission  is  insignificant  in 
the  airgjow  layer  at  100  km. 

Ozone  chemistry  dominates  die  night  sky’s  OH  Meinel  band  emissions 
(which  are  caused  by  vibrational  transitions  in  the  ground  electronic  state  of 
the  OH  radical): 

H  +  O3  ^  OH  (X  ^n3/2.i/2,  V<  9)  +  O2  +  3.3  eV. 

Tliis  vibrational  excitation  is  followed  by  the  reaction: 

OH  (X  ^n3/2,i/2,  V<  9)  OH  (X  ^U^|2M2)  +  hv  (Meinel,  5000  A  <  <  4  pm). 

No  OH  vibrational  band  emissions  fi:om  levels  above  ^=9  are  observed  in  the 
niglit  airglow  (Rees,  1989).  The  rotational  intensity  distribution  of  the  OH 
Meinel  system  can  be  used  as  a  measure  of  mesopause  temperature  and  its 
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variability,  but  the  vibrational  distribution  of  this  system  is  not  yet  fully 
quantified. 

The  Jolinston-Broadfoot  model  simulates  niglit  airglow  chemistry  from 
60  to  120  km  altitude,  dividing  the  atmosphere  into  30  finite  plane-parallel 
layers,  each  of  2  km  deptli,  and  using  a  time  step  of  300  seconds  (Johnston 
and  Broadfoot,  1993),  Each  layer  is  homogeneous  with  a  constant  particle 
number  density  for  tlie  major  constituents  O2  and  N2,  based  on  die 
assmnption  of  negligible  diurnal  variation  of  tiiese  constituents  at  the  altitudes 
in  question  (Rodrigo  et  al.,  1991).  Output  consists  of  vertical  number  density 
profiles,  integrated  vertical  column  densities,  and  altitudes  of  peak  emission 
for  the  lines  and  bands  of  O2,  OH,  and  the  various  species  of  O  (Johnston, 
1992). 


The  nmnber  densities  of  die  minor  constituents  are  calculated  using  die 


St = + P(Ni)-L(Ni) , 


equation  of  continuity; 
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where  tlie  number  density  of  the  minor  constituent  is  Nj,  the  chemical 
production  and  loss  terms  for  the  minor  constituent  are  P(Ni)  and  L(Ni), 
respectively,  t  is  time,  z  is  altitude,  and  the  tenn  didz  (Nj  Wj)  represents 
difiusion.  The  upward  diffusion  velocity  Wj  of  a  minor  constituent  is  given  by 


Wi 


^  Mjg  ^  1 


Ni  ^  z  k  T 


- )-k(^ 


A[Ni]) 


where  D  and  K  are  tlie  molecular  and  eddy  diffusion  coefficients, 
respectively,  T  is  temperature,  k  is  Boltzmann's  constant,  Mj  is  the  molecular 
or  atomic  mass  of  die  minor  constituent,  z  is  altitude,  g  is  die  acceleration  of 
gravity,  and  die  number  densities  of  the  major  and  minor  constituents  are  N 
and  Nj,  respectively  (Johnston,  1992). 


In  summary,  dissociation  provides  die  atomic  oxygen  needed  to  form 
excited  states  of  O2,  wliich  in  turn  help  create  the  excited  state  O  (^So).  hi  die 
niglit  airglow,  the  excited  states  of  O2  emit  in  die  Herzberg,  Atmospheric,  and 
Chamberlain  bands,  while  the  excited  state  of  atomic  oxygen,  O  (^So),  emits 
at  2972  A  and  5577  A.  It  is  noteworthy  that  all  these  emissions  of  the  niglit 
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sky  are  measured  simultaneously  by  GLO  (see  Figure  1 .2).  The  OH  radical 
reacts  with  atomic  oxygen  to  form  O2,  and  this  product  is  cycled  to  create 
ozone.  Since  ozone  is  formed  in  the  three-body  association  O  +  O2  +  M  -» 
O3  +  M,  OH  exists  in  a  complex  balance  widi  ozone  creation  and  destruction, 
mediated  by  the  presence  of  H,  HO2,  O,  and  O2.  Quenching  dominates  the 
rotational  and  vibrational  distributions  of  niglit  airglow  species  and  has  been 
qualitatively  modeled.  However,  a  quantitative  description  is  still  lacking. 
The  Johnston-Broadfoot  model  is  among  tlie  best  models  yet  available  to 
sunulate  tlie  complex  chemistry  of  the  night  airglow. 
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1 .4  Recent  Observations 

Weill  (1967)  summarized  tlie  early  liistory  of  ground-based  airglow 
observations  in  equatorial  regions,  wliich  had  already  established  the 
existence  of  an  equatorial  minimum  in  5577  A  emission  and  seasonal 
variations  in  sodimn  and  OH  emissions.  Many  other  researchers  in  tlie  last 
twenty  years  have  described  latitudinal  variations  in  the  5577  A  emission  of 
atomic  oxygen  inferred  from  botli  ground-based  and  satelhte  observations. 
Brenton  and  Silverman  (1970)  and  Takaliaslii  et  al.  (1989)  described  strong 
midnight  minima  in  5577  A  emission  measured  from  ground  stations  in 
equatorial  South  America.  In  the  1970s,  Donahue  et  al.  (1973)  and  Reed  and 
Chandra  (1975)  found  deep  equatorial  minima  in  5577  A  emission  from  their 
analyses  of  data  taken  by  the  OGO  series  of  satellites.  Cogger  and  Khaneja 
(1979)  also  reported  a  persistent  equatorial  minimum  in  5577  A  emission 
inferred  from  data  collected  by  the  ISIS-2  satellite.  Using  observations  of  tlie 
5577  A  emission  from  die  Atmospheric  Explorer  satellite,  Yee  and  Abreu 
(1987)  found  a  post-local  sunset  emission  maximum  and  a  post-local  midniglit 
minimum  for  latitudes  witiiin  20°  of  the  equator  in  the  spring,  and  tiiey 
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inferred  tliese  effects  were  connected  to  atmospheric  tides.  Modelers,  such  as 
Akmaev  and  Shved  (1980)  and  Forbes  et  al.  (1993),  have  seen  that  diurnal 
and  semidiurnal  atmospheric  tides  can  create  variations  by  factors  of  4  or 
more  in  5577  A  emission  and  induce  atomic  oxygen  density  reductions  on  the 
order  of  30%  at  airglow  altitudes. 

Observations  from  tlie  Upper  Atmosphere  Research  Satellite  (UARS) 
in  the  last  few  years  have  provided  interesting  clues  to  tlie  nature  of  the  niglit 
airglow.  Shepherd  et  al.  (1995a;  1995b)  and  Ward  et  al.  (1995)  have  used  tlie 
Wind  Imaging  Interferometer  (WINDII)  on  UARS  to  observe  nightglow 
emissions  in  5577  A  and  the  O2  Atmospheric  band.  Tliey  infer  vertical 
motions  in  die  MTI  region  and  spatial  variability  in  tlie  two  emissions.  These 
UARS  researchers  believe  tlieir  data  indicate  tliat  the  two  emissions  are  well 
correlated  in  space  and  time  (Ward  et  al.,  1995).  Tliey  also  intend  to  report 
findings  for  the  OH  Meinel  band  in  the  near  fiiture. 

The  persistent  emissions  of  the  night  sky  that  dominate  tliis  study  are 
the  O2  Atmospheric  bands,  die  OH  Meinel  bands,  and  die  01  (5577  A)  line  of 
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atomic  oxygen.  Figure  1.2  shows  tlie  averaged  night  airglow  spectrum 
obtained  by  the  GLO  instrument  dming  Space  Shuttle  mission  STS-53  in 
December  1992.  Note  tliat  tlie  complete  spectral  range  shown  was  recorded 
simultaneously  by  tlie  GLO  spectrographs. 

The  GLO  experiment  described  in  tlie  following  chapters  observed 
emission  intensity  variations  and  signatures  of  vertical  motions  in  the  niglit 
airglow.  Tliese  signatures  were  recorded  simultaneously  in  O2  (0,0),  O2  (0,1), 
O2  Herzberg,  O2  Chamberlain,  OI  5577  A,  and  OH  Meinel  emissions 
between  28.4°  N  and  28.4°  S  (geocentric  latitude)  under  equinox  conditions. 
In  contrast  with  the  conclusions  of  UARS  researchers,  the  GLO  data  exlubit  a 
distinct  lack  of  clear  correlation  among  niglit  sky  emissions. 

Advances  in  detector  and  computer  teclmology  in  the  last  decade  have 
made  comprehensive  observations  of  the  night  sky  possible.  The  night 
airglow  spectrum  by  Broadfoot  and  Kendall  (1968)  required  over  180  hours 
of  exposure  using  a  ground-based  spectrograph.  Their  spectrum  has  remained 
die  definitive  standard  for  the  aeronomy  community,  but  the  current  version  of 


spectrum  as  recorded  by  the  Arizona  Imaging  Spectrograph  (GLO)  on  STS-53 
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tlie  GLO  instrument  can  produce  an  airglow  spectrum  comparable  to 
Broadfoot  and  Kendall’s  product  in  about  10  to  15  minutes!  We  are  only 
beginning  to  systematically  study  the  niglit  sky  from  the  “top  side”  with  such 
a  capability.  No  previous  study  has  been  able  to  accumulate  continuous  and 
simultaneous  observations  of  die  various  airglow  emissions  on  a  global  scale. 

Tlie  data  used  in  tliis  research  were  obtained  during  Space  Shuttle 
mission  STS-69,  wliich  was  launched  on  7  September  1995  and  remained  on 
orbit  for  nearly  11  days.  The  orbital  geometry  of  STS-69  was  particularly 
iisefiil  for  low  latitude  niglit  airglow  studies  (Figure  1.3),  and  provided  a 
nearly  constant  observation  geometry  for  about  39  orbits  (orbits  66  tlirougli 
104)  during  solar  minimum  and  under  nearly  equinox  conditions.  Figure  1.4 
shows  representative  nighttime  orbit  tracks  (71  dirough  80)  to  indicate  STS- 
69  observational  coverage.  Note  tliat  each  orbit  track  is  22.5®  further  west 
(for  a  separation  of  about  2500  km  at  100  km  altitude)  tlian  die  track  diat 
preceded  it,  and  proceeds  from  the  southwest  toward  the  northeast  during  a 
niglit  pass. 
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Figure  1.3.  Orbital  geometry  for  Space  Shuttle  flight  STS-69. 
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2.  THE  GLO  EXPERIMENT  ON  STS-69 

Tlie  GLO  instrument  was  designed  to  observe  weak  emissions, 
including  aurora  and  airglow,  from  the  Space  Shuttle  and  satellites  (Broadfoot 
and  Sandel,  1992;  Broadfoot,  et  al.,  1992).  On  die  fourth  day  of  Mission 
Elapsed  Time  (MET)  for  STS-69,  the  Space  Shuttle  began  an  experiment 
wliich  required  a  constant  attitude  for  the  next  40  orbits.  The  LPL  West 
science  team  took  this  opportunity  to  continuously  observe  the  Earth’s  limb. 
The  GLO  experiment  proceeded  for  39  orbits  (orbits  66  tiirough  104),  and  die 
niglit  airglow  data  recorded  for  that  period  are  examined  in  this  research. 

Tlie  pointing  of  the  GLO  instrument  and  the  stability  of  the  Space 
Shuttle  platfonn  are  key  to  quantifying  diis  data  set.  GLO  can  be  pointed  with 
an  accuracy  of  ±0.005®.  GLO  tracking  errors  and  shuttle  axes  drifts 
contribute  to  altitude  errors  along  the  line-of-sight.  Knowledge  of  the  GLO 
observational  geometry,  the  shuttle’s  attitude  in  space,  and  die  pointing 
accuracy  of  the  GLO  instrument  is  necessary  before  the  data  recorded  on 
STS-69  can  be  analyzed  with  confidence. 
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2.1  Overview  of  tlie  GLO  Instrument 


Figure  2.1.  STS-69  payload  bridge  installed  in  the  shuttle  bay. 

GLO  spectrographs  and  imagers  are  housed  in  a  sensor  head  that  is 
mounted  on  a  motor-driven  scan  platfonn,  pennitting  motion  of  sensor 
boresiglit  in  elevation  (EL)  and  azimuth  (AZ)  relative  to  the  mounting  surface 
of  the  spacecraft.  The  two  stepping  motors  provide  rotation  in  AZ  and  EL  by 
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turning  worm  gears.  The  rate  of  sensor  head  movement  in  azimuth  is  3°  per 
second  with  0.0082®  resolution,  and  in  elevation,  3®  per  second  with  0.014® 
resolution.  The  total  range  of  sensor  head  motion  is  -60®  to  +135®  in  elevation 
and  -175®  to  +175®  in  azimuth  from  a  predefined  zero  position. 

Because  of  the  fit  of  tlie  worm  gear,  a  “backlash,”  or  position  error, 
occurs  whenever  the  stepping  motors  reverse  direction.  Backlash  control 
compensates  for  tliis  by  positioning  die  scan  platform  by  unidirectional 
rotation  in  the  positive  sense.  For  each  commanded  rotation  in  the  negative 
sense,  tlie  platform  is  overdriven  by  an  additional  0.5®  rotation,  followed  by  a 
positive  rotation  of  0.5°  before  stopping. 

The  GLO  sensor  head  consists  of  nine  spectrograph  channels,  as  well 
as  tln’ee  imagers  with  coaligned  optical  axes.  The  GLO  instrument  has  been 
previously  described  in  detail  by  Broadfoot  et  al.  (1992).  GLO  imager 
configuration  for  STS-69  was  changed  somewhat  from  tliat  description,  but 
die  spectrographs  remain  die  same  (see  Appendix  A).  GLO  records  airglow 
emissions  using  intensified  charge  coupled  devices  (ICCDs)  as  detectors  in 
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tlie  spectrographs.  The  apphcation  of  ICCDs  to  instrumentation  for  observing 
tlie  airglow  has  been  outlined  by  Broadfoot  and  Sandel  (1992).  The 
spectrographs  are  paired  so  tliat  two  gratings  are  illuminated  by  a  single 
entrance  slit.  Tlie  spectra  from  adjoining  spectral  regions  are  dispersed  side- 
by-side  onto  a  single  ICCD.  Tlie  nine  gratings  have  spherically  concave 
surfaces  that  have  been  holograpliically  ruled  to  correct  aberrations  and 
provide  a  slit  image  on  tlie  flat  ICCD  surface. 

The  ICCD  spectrographs  measure  spatial,  temporal,  and  spectral 
emission  signatures  at  wavelengtlis  from  1150  to  9000  A.  Spectrograph 
resolution  in  the  short  wavelengtli  region  is  about  4  A,  wliile  it  is  about  10  A 
at  tlie  long  wavelengths.  GLO  foreoptics  give  all  spectrographs  a  field-of- 
view  (FOV)  of  8.5°  by  0.15°,  and  permit  the  measurement  of  brightness 
variations  in  all  airglow  emissions  simultaneously.  Simultaneous  observation 
is  important  for  understanding  tlie  relationsliips  among  the  emissions. 

Spatial  infonnation  is  recorded  along  tlie  length  of  tlie  spectrograph 
slit.  The  ICCD  digital  “readout”  is  adjusted  electronically  to  give  24 
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contiguous  spatial  elements  or  “superpixels”  of  8  pixels  each  along  the  slit, 
per  spectral  unage.  An  IR  imager  assisted  the  observational  program  by 
tracking  the  airglo^v  layer  as  it  was  viewed  edge-on  at  the  Earth’s  limb. 


For  the  purpose  of  identifying  tlie  nine  spectral  images,  tlie 
spectrographs  have  been  labeled  as  A  through  E,  from  the  longest  wavelength 
to  the  shortest.  The  spectral  images  are  then  numbered  one  dirougli  nine  as 
outlined  in  Table  2.1. 


SPECTROGRAPH 

SPECTRAL  IMAGE 

SPECTRAL  RANGE 

A 

1 

9000-11000  A 

B 

2 

7500-9000  A 

3 

6000-7500  A 

C 

4 

5250-6000  A 

5 

4500-5250  A 

D 

6 

3750-4500  A 

7 

3000-3750  A 

E 

8 

2100-3000  A 

9 

1100-2100  A 

Table  2.1.  GLO  spectrograph  and  spectral  image  assignments. 
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For  spectrograph  A,  a  single  grating  covered  tlie  spectral  range  for 
channel  1.  This  spectrograph  was  not  sensitive  enougli  to  produce  useful  data. 
The  emissions  of  interest  here  were  recorded  by  spectrographs  B,  C,  and  E. 
The  ICCD  in  spectrograph  D  failed  shortly  after  launch.  Intensified  imager 
data  was  not  examined  in  this  work.  The  unintensified  IR  imager  was  used  to 
assist  in  spectrograph  pointing.  In  our  case,  it  was  also  used  for  tracking  tlie 
airglow  limb  and  for  GLO  boresigjit  calibration. 

Broadfoot  and  Sandel  (1992)  have  described  the  advantages  of  ICCDs 
over  bare  CCDs,  scanning  monocliromators,  and  Ebert-Fastie  spectrometers. 
ICCDs  maximize  tlvougliput  over  a  broad  spectral  range  and  exliibit  much 
lower  detector  noise  compared  to  historical  spectrometers.  The  small  size  of 
ICCDs  permits  much  smaller  instrumentation,  even  wliile  allowing  a  broader 
spectral  range  to  be  simultaneously  recorded.  The  reduction  in  instrument  size 
is  typically  an  order  of  magnitude,  and  tliis  reduces  the  sensitivity  to  interior 
scattered  light  by  tlie  same  factor.  GLO  spectrographs  contain  no  moving 
parts,  while  spectrometers  contain  rotating  gratings  and  stepping  motors. 
Witliin  a  GLO  spectrograph,  die  geometry  of  the  optics  and  ICCD  detector  is 
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fixed,  giving  a  well  defined  wavelength  scale,  whereas  in  a  spectrometer, 
accuracy  of  wavelengtli  determination  is  dependent  on  tlie  precision  of  tlie 
internal  stepping  motor  that  drives  the  rotating  diffraction  grating.  The 
reduction  in  die  size  of  the  foreoptics  and  detector  element  means  the  ICCD 
spectrograph  in  GLO  can  have  a  much  smaller  mass  and  volume  than  a 
comparably  perfonning  spectrometer,  a  significant  advantage  in  a  space  flight 


instrument. 
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2.2  Observational  Geometry 

The  niglit  airglow  is  known  to  originate  in  relatively  tliin  emitting 
layers  witli  thicknesses  on  tlie  order  of  10  km.  Therefore,  it  is  important  to 
determine  where  the  GLO  boresiglit  is  pointing  when  a  spectral  image  is 
recorded.  Figure  2.2  shows  tlie  airglow  viewing  geometry  for  a  typical  niglit 
pass.  Observations  were  taken  looking  into  the  shuttle  wake,  opposite  the 
direction  of  motion,  with  the  slit  tilted  about  30®  to  the  limb.  The  24  spatial 
elements  along  the  slit  are  depicted. 


Figure  2.2.  GLO  airglow  viewing  geometry  showing  spectrograph  slit  and  limb  tracking 
imager  Field  Of  View  (FOV). 
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At  a  slant  range  of  about  2000  km  from  the  airglow  layer,  the  field  of 
view  (FOV)  of  a  spectrograph  is  about  296  by  5  km.  The  slit  image  has  a 
lengtli  of  192  pixels  on  tlie  ICCD  and  is  divided  into  24  subgroups  (spatial 
elements  or  “superpixels”)  measuring  8  by  1  pixels  (Figure  2.2),  If  the  slit  is 
tilted  at  an  angle  to  tlie  limb,  vertical  resolution  when  sampling  die 
atmosphere  is  improved.  At  an  angle  of  10°  to  the  horizontal  and  a  slant  range 
of  2000  km,  a  51  km  depth  of  tlie  atmosphere  is  sampled.  With  a  slit  tilt  of 
30°  and  2000  km  slant  range,  GLO  surveys  a  148  km  depth,  wliile  at  90°,  a 
maximum  altitude  rauge  of  about  296  km  can  be  imaged  on  the  slit. 

Resolution  in  tlie  vertical  depends  on  superpixel  separation  along  the 
vertical  dimension  of  the  image.  A  horizontal  slit  will  create  no  vertical 
separation  for  the  spatial  elements;  this  gives  optimal  signal  strength,  but  only 
for  a  shallow  layer  about  5  km  tliick.  Althougli  vertical  resolution  cannot  be 
smaller  tlian  tlie  slit  widtli  (5  km  at  2000  km  slant  range),  signal  strength  is 
increased  by  minimizing  tlie  airglow  deptli  over  wliich  the  slit  is  imaged.  A 
vertically  oriented  slit  will  provide  maximum  vertical  spatial  resolution,  witli 
each  spatial  element  representing  an  altitude  increment  of  about  12  km,  but 
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since  this  covers  a  296  km  depth,  signal  strengtli  at  any  given  altitude  is 
minimized.  For  die  STS-69  altitude  of  400  km,  GLO  sht  tilt  of  30°,  and  slant 
range  to  the  tangent  point  of  2000  km,  the  vertical  separation  between  the 
centers  of  spatial  elements  is  about  6.8  km. 

These  factors  show  the  importance  of  limb  tracking  of  the  airglow. 
Shuttle  motions  cause  the  spectrograph  slit  to  sweep  over  the  limb 
periodically  by  ±1®  (±35  km  at  2000  km  slant  range).  The  GLO  limb  tracker 
follows  the  brightest  region  of  a  strip  through  the  imager  FOV  by  moving  the 
scan  platform  in  AZ  and  EL,  dius  compensating  for  tliese  motions.  Limb 
tracking  represents  an  important  enliancement  in  GLO’s  capabilities. 

Figure  2.3  shows  the  observation  slant  path  across  die  Earth’s  limb  to 
scale.  Angular  dimensions  given  are  for  die  maximum  vertical  extent  of  die 
spectrograph  and  imager  FOVs  widi  the  slit  tilted  at  30°.  The  dots  along  die 
airglow  layer  indicate  imager  exposures  at  8  second  intervals.  Spectrograph 
exposures  occurred  at  20  second  intervals,  so  there  were  two  spectrograph 
exposures  for  every  five  dots  indicated  in  the  figure. 
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second  intervals.  Compare  with  Figure  2.4. 
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Figure  2.4  below  shows  a  typical  view  of  tlie  niglit  airglow  layer  using 
GLO’s  IR  limb  tracking  imager. 


Spectrograph  Slit 


Tangent  Height  (km) 
160 


Slant  Range  (km) 


Figure  2.4.  The  night  airglow  in  GLO’s  limb  tracking  imager.  Tangent  height,  slant  range, 
GLO  AZ  and  EL,  and  vertical  extent  of  the  FOVs  of  the  imager  and  spectrograph  slit  from 
STS-69  are  indicated.  A  slit  tilt  of  30°  is  depicted.  Compare  with  Figure  2.3. 


Tlie  advantage  of  controlling  the  slit  orientation  is  well  illustrated  by 
Figure  2.4.  Note  that  half  of  tlie  spectrograph  slit  is  above  the  airglow  layer 
and  is  therefore  not  collecting  useful  data.  The  optimmn  orientation  would 
have  had  tlie  topmost  end  of  tlie  slit  just  above  tlie  airglow  layer.  However, 
this  first  use  of  limb  tracking  by  GLO  required  that  the  limb  be  near  the  center 
of  the  IR  imager  FOV.  This  will  be  modified  for  future  flights. 
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Tlie  geometry  used  to  collect  tliis  data  set  has  been  depicted  in  Figures 
2.2  tlu'ough  2.4.  It  is  also  important  to  understand  the  GLO  coordinate  system 
as  depicted  in  Figure  2.5  on  the  next  page.  The  GLO  azimuth  axis  is  along  the 
“+Z”  shuttle  axis,  where  “-Z”  points  toward  die  velocity  vector  (see  Figure 
2.6  in  the  next  section).  Thus  die  azimuth  axis  points  in  the  shutde  wake.  The 
Eardi  limb  seen  from  our  400  km  altitude  was  depressed  about  18°  from  dus 
axis.  Figure  2.5  shows  that  the  slit,  wliich  is  tangent  to  circles  of  constant 
elevation,  can  be  oriented  at  any  angle  to  die  limb  between  0°  and  90°  by  die 
proper  selection  of  GLO  AZ  and  EL. 
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GLO  ASEL  COORDINATE  SYSTEM 


Figure  2.5.  “Fish-eye”  view  of  the  GLO  scan  platform  coordinate  system  showing  the 
Earth  and  stars.  The  Earth  limb  is  depressed  about  18°  from  GLO  zenith.  Note  the  slit  tilt 
of  30°  at  the  limb,  the  orientation  of  the  Shuttle  axes,  the  star  in  the  IR  Imager  FOV,  and 
the  outline  of  the  shuttle  bay  as  seen  by  GLO. 
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2.3  Space  Shuttle  Attitude  and  GLO  Pointing  Accuracy 

The  shuttle  attitude,  and  hence  the  GLO  observation  line-of-siglit 
(LOS),  are  referenced  to  a  gyroscopic  platform  known  as  an  Inertial 
Measurement  Unit  (IMU).  Shuttle  attitude  is  maintained  by  the  IMU,  which 
provides  inertial  attitude  data  used  to  control  die  shuttle  flight  surfaces,  the 
tiirust  vector  for  die  main  engine,  and  the  reaction  control  system  verniers. 
The  IMU  permits  very  accurately  shutde  attitude  determination.  Aldiougli  die 
shuttle  pitch,  roll,  and  yaw  axes  drift  between  attitude  limits  called  die 
“deadband,”  shutde  instantaneous  attitude  is  well  known  (see  Figure  2.6  on 
the  next  page  for  definitions  of  shuttle  pitch,  roll,  and  yaw).  For  STS-69,  die 
deadband  was  ±1.0°  in  shuttle  rod  and  pitch,  and  ±3.0°  in  yaw. 

The  shutde  also  experiences  slight  attitude  variations  because  of  slow 
gyrocompass  drift  in  die  IMU.  Tlie  IMU  uses  accelerometers,  gyros,  and 
electromechanical  devices  to  maintain  attitude.  All  are  prefliglit  calibrated  to 
fine  tolerances,  but  the  gyros  require  periodic  on-orbit  alignment  for  drift. 
Typically,  die  IMU  is  recalibrated  daily  against  die  star  field,  but  this  is  not 
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Yaw  Pitch 


Figure  2.6.  Shuttle  pitch,  roll,  and  yaw  axes. 

always  possible.  On  the  STS-69  mission,  an  uncorrected  IMU  gyrocompass 
drift  of  0.14°  occurred  in  a  47  hour  period.  This  drift  represented  an  error  of 
about  5  km  in  some  direction  with  respect  to  the  Earth’s  limb  (aldiough  not  all 
of  die  drift  would  have  been  in  the  pitch  direction  to  constitute  an  altitude 
error).  The  IMU  was  reset  about  midway  through  our  39  orbit  run.  NASA 
mission  control  provided  the  times  of  realignments  during  the  mission,  so  the 
mavimiim  error  in  pointing  due  to  diis  drift  is  known  as  a  function  of  time. 
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Four  mounting  gimbals  isolate  the  IMU  from  shuttle  rotations  in  pitch, 
roll,  and  yaw,  in  order  to  maintain  inertial  attitude.  One  IMU  is  sufficient  for 
maintaining  attitude,  but  three  are  installed  for  redundancy.  They  are 
intentionally  not  coaligned  and  not  aligned  with  die  shuttle  pitch,  roll,  and 
yaw  axes.  This  skewed  arrangement  provides  redundancy  for  shuttle  attitudes 
in  wliich  a  single  IMU’s  gimbals  would  not  be  effective  and  pennits  cross¬ 
checking  of  IMU  failures.  However,  because  of  tlie  skewed  IMU  axes 
arrangement  and  handovers  among  IMUs  in  real  time,  the  direction  of 
gyrocompass  drift  with  respect  to  tlie  pitch,  roll,  and  yaw  axes  is  not  a 
parameter  readily  available  from  NASA.  Tlie  worst  case  scenario  would  be 
for  all  IMU  drift  to  occur  in  tlie  pitch  axis.  Appendix  C  discusses  the  effects 
of  shuttle  pitch  variations  on  GLO  pointing  and  limb  tracking. 

The  full  excursion  of  axis  drift  within  tlie  deadband  is  known  as  tlie 
shuttle  “limit  cycle.”  Figures  2.7  and  2.8  show  tlie  STS-69  limit  cycle  over 
two  full  orbits,  indicating  tlie  drift  in  roll,  pitch,  and  yaw  over  two  90  minute 
periods.  NASA  provides  post-fliglit  attitude  data  with  one  second  time 
resolution  for  precise  determination  of  shuttle  orientation.  NASA  data  give 
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shuttle  pitch,  yaw,  and  roll  in  degrees  to  at  least  ±0.0001°.  Previous  flight 
experiments  (Neupert  et  al.,  1992)  have  shown  that  the  overall  accuracy  of 
attitude  determination  along  any  axis  is  at  least  ±0.001°.  Attitude  variations 
on  one  second  time  scales  are  smootli  to  witliin  seconds  of  arc. 

During  the  collection  of  our  39  orbit  data  set,  the  shuttle  velocity  vector 
was  in  die  direction  of  tlie  -Z  axis  (“belly  to  ram”)  and  GLO  viewed  out  of 
tlie  bay  into  tlie  shuttle  wake.  The  fliglit  dynamics  are  shown  by  the  limit 
cycle  motion  indicated  in  the  previous  two  plots.  Recall  the  limit  cycle  is  the 
excursion  of  attitude  drift  with  respect  to  the  IMU.  Perturbing  forces  in  tlie 
pitch,  roll  and  yaw  planes  are  coimterbalanced  by  vernier  thruster  firings 
whenever  die  pitch  and  roll  axes  vary  ±1°  from  nominal  and  whenever  die 
yaw  axis  varies  ±3°  from  nominal.  Dynamic  imbalances  in  pitch  and  roll  are 
greater  dian  the  imbalance  in  yaw.  Cross-coupling  of  pitch  and  roll 
corrections  causes  die  yaw  motion  to  be  irregular. 
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With  the  GLO  scan  platfonn  motionless,  changes  in  LOS  pointing  can 
be  retrieved  from  knowledge  of  shuttle  attitude  to  an  accuracy  of  ±0.005°,  or 
about  ±175  meters  at  tlie  limb.  However,  removing  uncertainty  in  the 
alignment  of  tlie  shuttle  axes  and  the  GLO  scan  platform  AZTEL  axes  requires 
a  calibration  of  tlie  two  systems  with  respect  to  the  star  field.  The  offset 
between  actual  GLO  boresight  pointing  and  the  commanded  AZ/EL  position 
is  obtained  by  tlie  metliod  described  in  tlie  next  section. 
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2.4  Pointing  Calibration  Against  the  Star  Field 

Many  stars  pass  through  tlie  IR  imager’s  FOV  on  each  orbit.  For 
pointing  calibration,  an  IR  imager  frame  containing  a  known  star  is  chosen. 
Each  image  frame  recorded  by  the  IR  imager  contains  header  information 
giving  the  commanded  boresight  AZ  and  EL  at  die  time  of  observation.  Dr. 
Bill  Sandel  of  LPL  West  prepared  an  algorithm  to  compute  star  position  in 
“true”  GLO  AZ/EL  coordinates,  given  a  known  star’s  astronomical  riglit 
ascension  (RA)  and  declination  (Dec)  coordinates.  This  is  possible  because 
die  shuttle  coordinate  system  and  die  RA  and  Dec  system  are  related  by  a 
known  rotation. 

If  a  star  is  located  exactly  at  die  boresight  position  in  the  imager’s  pixel 
array.  Dr.  Sandel’s  program  will  immediately  give  an  offset  for  the  boresight 
position  with  reference  to  die  Shuttle  coordinate  system,  since  the  star’s 
position  is  well  known.  Generally,  a  star’s  position  will  not  be  so  fortuitous, 
and  angular  separation  across  the  FOV  must  be  calculated  in  both  coordinate 
systems.  If  die  FOV  of  die  unager  is  near  die  GLO  azimudi  axis  (or  RA  axis) 
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of  tlie  coordinate  system,  the  strong  curvature  of  lines  of  constant  GLO 
elevation  (or  Dec)  prevents  the  simple  detennination  of  angular  separation  by 
overlay  of  the  grid  system  on  the  imager’s  rectangular  pixel  coordinates.  Dr. 
Emmet  Anderson  of  LPL  West  provided  a  FORTRAN  program  to  account  for 
this  effect.  Using  the  algorithms  of  Dr.  Sandel  and  Dr.  Anderson,  the  “true” 
AZ/EL  coordinates  of  the  boresiglit  pixel  iu  the  IR  imager’s  CCD  array  can 
be  determined,  since  we  know  the  AZ/EL  pixel  scale  (57  pixels  per  degree  in 
tliis  case)  and  we  can  measure  the  angular  separation  between  the  star  and  tlie 
boresiglit.  Tlie  difference  between  die  “true”  AZ/EL  boresiglit  position  and 
tlie  commanded  AZ/EL  boresiglit  position  gives  the  pointing  error  between 
the  GLO  AZ/EL  coordinate  system  and  the  shuttle  coordinate  system.  The 
true  boresight  position  at  one  specific  time  was  determined,  based  on  an  IR 
image  with  the  stars  ^  Sgr,  o  Sgr,  tt  Sgr,  and  one  unnamed  star  in  Sagittarius 
(see  Figure  2.9  on  tlie  next  page;  note  that  NASA  uses  celestial  coordinates 
for  Julian  epoch  1950). 

However,  comparison  of  this  true  value  witli  the  commanded  boresight 
position  was  not  possible  because  of  uncertainty  in  the  precise  timing  of  the 
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Figure  2.9.  IR  image  used  in  boresight  calibration. 


exposure.  GLO’s  internal  clock  labels  spectrograph  and  imager  exposures, 
but  it  lags  behind  the  much  more  accurate  orbiter  clock.  Orbiter  mission  time 
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is  used  to  determine  the  position  of  tlie  Shuttle  and  tlius  observation  geometry 
and  the  tangent  point  height  along  the  LOS  direction.  Only  the  GLO  time 
label  is  knovm  for  a  given  imager  exposure,  not  the  orbiter  mission  time.  A 
five  second  offset  between  GLO  time  and  orbiter  time,  with  GLO  lagging 
behind,  was  typical  over  a  24  hour  period.  Since  die  tangent  point  heiglit 
changes  by  about  2.4  km  each  second,  a  discrepancy  of  only  5  seconds 
between  GLO  time  and  orbiter  time  would  amount  to  an  altitude 
determination  error  of  about  12  km.  The  GLO  clock  is  ideally  reset  to  mission 
time  every  24  hours,  but  this  was  not  always  done  during  mission  STS-69, 
and  errors  greater  than  5  seconds  were  possible. 

Thus,  the  most  serious  obstacle  to  determination  of  tlie  absolute 
pointing  accuracy  of  the  GLO  boresight  in  this  data  set  is  the  GLO  clock. 
With  some  effort,  tlie  drift  of  the  clock  can  be  estimated  and  propagated  from 
tlie  last  clock  reset.  In  the  future,  it  is  recommended  that  tlie  GLO  clock  be 
reset  to  orbiter  time  on  a  regular  basis  so  tliat  die  drift  never  exceeds  one 
second.  In  this  work,  airglow  altitudes  have  been  corrected  with  respect  to  the 
sodium  layer  at  90  km,  assumed  to  be  a  stable  reference  altitude. 
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3.  DATA  ANALYSIS  AND  OBSERVATIONS 

Since  GLO  observes  a  wide  spectral  range  simultaneously  from  a 
column  of  gas  along  die  line-of-sight,  a  large  amount  of  information  is 
recorded  in  each  exposure.  Juxtaposing  a  series  of  spectrograph  exposures  as 
the  instrument  FOV  sweeps  along  the  orbit  track  creates  a  three-dimensional 
data  array  of  column  emission  intensities  in  rayleighs  as  a  function  of 
wavelengtii,  altitude,  and  elapsed  time  (or  distance).  Such  an  array  is  called  a 
“monochromatic  image.”  These  arrays  can  be  studied  in  many  ways.  Tlie 
approach  taken  here  is  to  extract  monocliromatic  images  of  night  sky 
emissions  in  die  O2  Atmospheric  bands,  the  OH  Meinel  bands,  and  OI  (5577 
A),  showing  the  behavior  of  these  emissions  along  the  orbital  track. 

GLO  spectra  provide  column  emission  intensities  in  rayleighs  for 
various  emissions  in  die  airglow  layer  above  die  Eardi’s  limb.  An  Abel 
integral  inversion  algoritlun  (Appendix  D)  is  then  applied  to  these  column 
intensities  in  order  to  obtain  volmne  emission  rates  in  photons  per  cubic 
centimeter  per  second  and  to  construct  vertical  emission  profiles.  Such 
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profiles  show  tlie  behavior  of  the  different  airglow  emissions  as  a  fimction  of 
altitude  and  time.  From  such  profiles,  we  can  infer  how  tlie  airglow  layer 
changes  in  intensity,  heiglit,  and  tliickness,  and  by  knoAving  the  chemistry 
involved  in  the  different  emission  processes,  test  proposed  mechanisms  to 
explain  tlie  observed  variations. 

In  tlie  discussion  tliat  follows,  data  reported  or  plotted  in  rayleighs  has 
not  been  subjected  to  tlie  Abel  inversion  process,  wliile  quantities  reported  or 
plotted  in  photons  per  cubic  centimeter  per  second  represent  inverted  data. 
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3.1  Spectral  Imaging  and  Monocliromatic  Images 
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Figure  3.L  STS-69  GLO  spectrum  showing  O2  Atmospheric  and  OH  Meinel  bands. 


The  holograpliically  ruled  concave  grating  in  each  spectrograph  forms 
an  image  of  tlie  entrance  slit  dispersed  in  die  image  plane  of  the  detector.  Tlie 
ruling  corrects  aberrations  and  flattens  tlie  field  on  the  detector,  tlius 
preserving  spatial  information  along  tlie  slit.  For  the  data  recorded  here,  die 
slit  accepts  a  one-dimensional  portion  of  die  image  of  the  airglow  layer 
produced  by  die  foreoptics  and  fi-om  diis,  a  spectrum  is  fonned  showing 
spectral  intensity  versus  wavelength  and  altitude  (Figure  3.1). 
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Figure  3.2  below  shows  a  spectral  image  of  O2  Atmospheric  and  OH 
Meinel  band  emissions  and  a  plot  of  the  intensity  versus  wavelength.  The  two 
ciuves  in  Figure  3.2  are  the  same.  One  is  given  at  fiill  scale,  and  the  other  is 
reduced  by  a  factor  of  5  in  order  to  show  the  relative  intensity  of  tlie  O2  (0,0) 
band  compared  to  other  emissions. 


7^100  7600  7800  8000  8200  SW  8600  8800  9000 

ANGSTROMS 

Figure  3.2.  Plot  of  intensity  versus  wavelength  for  a  typical  spectral  image  from  GLO 
spectrograph  B  (see  Table  2.1). 
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Monocliromatic  altitude  profiles  of  selected  emissions  can  be 
constructed  fi-om  the  spectral  images.  Intensity  for  the  specific  spectral 
feature  of  interest  is  determined  from  the  calibrated  data.  A  spectral  sample 
and  background  sample  are  illustrated  in  Figure  3.3.  The  defined  spectral 
region  is  for  the  O2  (0,0)  band  emission.  The  intensity  is  summed  over  the 
region.  Then  the  background  intensity  region  is  defined,  scaled,  and 
subtracted  from  the  O2  (0,0)  intensity  region,  resulting  in  the  total  band 
intensity  in  rayleiglis  for  the  O2  (0,0)  band.  This  procedure  is  equivalent  to 
applying  a  monocliromatic  bandpass  filter,  except  in  tliis  technique,  the 
equivalent  filter  has  an  adjustable  bandpass  with  very  sharp  edges,  allows  all 
emissions  falling  within  tlie  bandpass  from  the  column  of  gas  in  view  to  be 
recorded  simultaneously,  and  permits  unwanted  emissions  and  backgroimd 
intensity  to  be  removed  easily  by  defining  the  bandpass  as  one  wishes. 
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Figure  3.3.  Digital  separation  of  signal  and  background  acts  as  a  bandpass  filter. 


Vertical  intensity  profiles  of  a  selected  emission  can  be  extracted  from 
each  spectral  image.  The  sampling  procedure  above  is  performed  on  each  of 
the  24  spectra  in  a  spectral  image,  resulting  in  24  monocliromatic  samples  at 
24  contiguous  altitudes.  These  vertical  intensity  profiles  of  the  original 
spectra  can  tlien  be  rearranged  in  a  side  by  side  fonnat  to  create  anotlier 
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array,  showing  intensity  of  tlie  selected  monochromatic  feature  versus  altitude 
and  horizontal  displacement  (or  equivalently,  elapsed  time)  along  the  orbital 
track.  This  monochromatic  image  is  effectively  “three-dimensional”  in 
altiUide,  intensity,  and  elapsed  time. 

Examples  of  monocliromatic  images  created  from  STS-69  spectra  are 
depicted  in  Figure  3.4,  and  in  Figures  3.5  tlirougli  3.7  in  negative  format.  The 
horizontal  dimension  of  these  monochromatic  images  is  time  along  the  orbit 
track,  fri  the  monochromatic  images  depicted  here,  each  vertical  intensity 
profile  represents  an  integration  of  20  seconds  on  the  detector.  The  relative 
briglitness  of  each  sample  is  represented  in  gray  scale,  forming  a 
monochromatic  image  showing  die  spatial  and  temporal  variations  in 
emissions  as  the  shuttle  moved  along  its  orbit  at  approximately  7.8  km/s. 
Since  the  gray  scale  of  the  images  is  not  quantitative,  the  emission  intensity 
along  the  orbit  track  at  constant  tangent  heiglits  is  plotted  in  a  graph  below 
each  image  in  Figures  3.5  through  3.7.  Note  tliat  representative  altitudes  for 
selected  rows  in  tlie  array  are  indicated  alongside  the  gray  scale  images  at  die 
top  of  Figures  3.5  dirough  3.7. 
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Figures  3.5  tlirough  3.7  show  the  line-of-siglit  column  brightnesses  of 
emissions  at  selected  constant  altitudes.  The  time  axis  of  each  briglitness  plot 
begins  at  the  moment  the  orbiter  crossed  the  terminator  into  shadow  along  the 
orbit  track,  relating  all  plots  in  local  solar  time.  Tlie  O2  (0,0)  and  01  (5577  A) 
plots  are  similar,  but  the  brightness  variations  of  the  two  emissions  do  not 
coincide.  Tlie  OH  Meinel  band  emission  is  not  correlated  with  eitlier  O2  (0,0) 
or  01  (5577  A).  Briglimess  changes  vary  witli  altitude  in  all  tiiiree  emissions. 

Signal-to-noise  ratios  for  O2  (0,0),  OH  Meinel,  01  (5577  A),  and  O2 
(0,1)  emissions  in  our  data  are  on  the  order  of  18,  10,  8,  and  6,  yielding 
briglitness  determination  errors  of  5%,  10%,  13%,  and  17%,  respectively. 
Siunming  N  altitude  profiles  into  one  improves  tlie  signal-to-noise  ratio  by  a 
factor  of  a  smmnation  of  9  slices  gives  a  factor  of  3  increase,  a 
summation  of  25  gives  a  factor  of  5,  and  so  on  (Broadfoot  and  Sandel,  1992). 
As  long  as  variations  among  the  slices  being  summed  are  minimized, 
information  loss  is  also  minimized.  The  “smearing”  tliat  does  occur  is 
equivalent  to  that  caused  by  increased  exposure  or  integration  time,  and 
represents  a  tradeoff  between  signal-to-noise  ratio  and  temporal  resolution. 
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Figure  3.6.  Monochromatic  image  and  observed  column  brightness  profile  for  the  01 
(5577  A)  emission. 


Tangent  Height 


Brightness 


79 


Minutes  after  tenninator  crossing  on  MET  day  4,  hour  13:41:05 
STS-69  Orbit  72 


Figure  3.7.  Monochromatic  image  and  observed  column  brightness  profile  for  the  OH 
Meinel  band  emission. 
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3.2  Signatures  of  the  O2  (0,0)  Band  in  the  GLO  Data  Set 

In  tliis  section,  we  will  discuss  O2  (0,0)  intensity  variations  seen  in  tlie 
night  airglow  for  29  orbits  (71  tlirough  99)  during  STS-69.  This  subset  of  die 
39-orbit  GLO  data  set  was  chosen  for  convenience  and  consistency,  since 
observation  geometry  was  less  stable  at  die  beginning  and  end  of  the  39-orbit 
data  set. 

The  atmosphere  is  optically  thick  to  emission  of  O2  (0,0),  so  die 
column  intensity  recorded  by  GLO  originates  between  die  tangent  point  and 
die  shuttle,  with  all  emission  from  the  far  side  of  the  line-of-si^t  being 
attenuated.  Tlius  GLO  observes  an  “O2  (0,0)  emission  surface”  restricted  by 
the  optical  depth  of  O2  (0,0).  Because  of  this  effect  and  because  it  is  the 
briglitest  airglow  emission,  the  O2  (0,0)  band  is  the  best  emission  to  use  in 
examining  wave  signatures.  On  the  odier  hand,  the  optical  depth  must  be 
considered  before  inverting  O2  (0,0)  intensities  to  obtain  volume  emission 
rates,  wliich  are  necessary  to  interpret  airglow  chemistry.  The  interpretation 
of  emissions  in  the  O2  Atmospheric  (0,0)  band  requires  radiative  transfer 
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theory  for  a  complete  description.  Solving  the  O2  (0,0)  radiative  transfer 
problem  will  not  be  attempted  here. 

The  variations  in  O2  (0,0)  intensity  are  of  interest  with  respect  to  their 
relationsliip  to  the  various  types  of  waves  tliat  may  exist  in  the  MTI  region. 
The  plot  of  O2  (0,0)  colimm  intensity  variations  along  the  track  of  orbit  72  in 
Figure  3.5  shows  typical  intensity  enliancements.  The  half-widtli  of  tlie 
intensity  peaks  is  about  six  minutes  of  orbiter  flight  time.  At  an  orbital 
velocity  of  7.8  km/sec,  this  half-width  represents  about  2800  km,  wliich  can 
be  taken  as  the  typical  extent  of  an  enhanced  airglow  region  along  the  orbit 
track.  With  each  ni^t  pass,  the  orbital  track  of  STS-69  shifted  22.5°  to  tlie 
west,  for  a  displacement  of  about  2500  km  (Figure  1.4).  If  tlie  intensity 
enliancements  are  part  of  a  long  hved  phenomenon,  tliey  may  be  detected  on 
consecutive  orbits.  With  a  characteristic  extent  between  2500  and  3000  km, 
tliey  would  just  span  the  displacement  between  night  passes.  An  isolated 
enliancement  miglit  be  observed  on  two  consecutive  orbits,  and  an  extensive 
wave  front  might  be  detected  over  several  orbits. 


82 


STS-69  ORBITAL  TRACK  COVERAGE  AS  A 
FUNCTION  OF  LATITUDE  AND  LONGITUDE 
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Figure  3.8.  Orbit  track  coverage  in  latitude  and  longitude  during  the  STS-69  mission. 
Orbit  71  is  the  first  orbit  at  the  right,  while  orbit  95  is  the  last  on  the  left. 


STS-69  OBSERVATION  POINTS  WITH  RESPECT  TO 
LATITUDE  AND  LOCAL  TIME 
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Figure  3.9.  GLO  observations  during  the  STS-69  mission  fell  along  an  orbital  track  fixed 
in  local  time. 
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Figure  3.8  shows  the  geograpliic  coverage  of  orbits  during  which  GLO 
data  was  recorded  and  illustrates  die  latitude  and  longitude  of  die 
observational  tracks.  Figure  3.9  shows  the  local  time  trace  of  GLO 
observations  during  night  passes  through  die  Eardi’s  shadow.  Observations 
are  seen  to  fall  on  a  single  orbital  track  fixed  widi  respect  to  local  solar  time 
(LSI)  in  Figure  3.9.  Note  that  die  orbital  track  is  also  fixed  in  latitudinal 
coverage  and  that  die  equatorial  crossing  occurs  at  0.5  LST  on  every  orbit. 

In  order  to  examine  die  data  set  for  evidence  of  organized  intensity 
enliancements  diat  may  be  related  to  tides,  the  intensity  along  die  orbital 
tracks  has  been  plotted  in  Figures  3.10  and  3.11  as  a  fimction  of  local  solar 
time.  These  two  Figures  show  O2  (0,0)  briglitness  in  rayleighs  for  29  orbits. 
Fifteen  orbits  are  plotted  in  each  figure,  with  orbits  71  dirough  85  in  Figure 
3.10  and  orbits  84  dirough  99  in  Figure  3.11.  The  orbits  overlap  between  the 
Figures.  Intensity  values  are  averaged  over  an  altitude  range  of  20  km, 
centered  at  90  km.  The  horizontal  direction  represents  about  80  observations 
per  orbit.  The  density  of  observations  in  this  plot  is  dius  much  higher 
horizontally  dian  vertically. 
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O2  (0,0)  BRIGHTNESS  IN  KILORAYLEIGHS 
AVERAGED  AT  LAYER  PEAK 


South  LATITUDE  North 


Figure  3.10.  Brightness  of  the  O2  (0,0)  emission  for  orbits  71  through  85  as  a  function  of 
orbit  and  local  solar  time.  There  is  no  evidence  of  tidally-induced  standing  waves. 
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South  LATITUDE  North 


Figure  3.11.  Brightness  of  the  O2  (0,0)  emission  for  orbits  84  through  99  as  a  function  of 
orbit  and  local  solar  time.  There  is  no  evidence  of  tidally-induced  standing  waves. 


SHY  (rayleighs) 


87 

AVERAGED  02(0,0)  INTENSITY  AS  A  FUNCTION 

OF  LOCAL  SOLAR  TIME  AND  LATITUDE 
FOR  29  ORBITS 
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Figure  3.12.  Intensity  of  the  O2  (0,0)  emission  averaged  over  orbits  71  through  99  as  a 
fiinction  of  local  solar  time.  Note  the  pronounced  intensity  minimum  near  1.0  LST,  a  few 
degrees  north  of  the  equator. 


88 


n 


If  tides  existed  in  Figures  3.10  and  3.11,  tliey  would  be  detected  as 
vertical  bands  of  constant  isopleths  fixed  in  LST.  However,  moving  vertically 
along  tliese  plots,  briglitness  values  change  significantly  firom  orbit  to  orbit, 
exhibiting  no  obvious  correlation  to  local  time.  The  implication  is  that  tidal 
variations  along  our  observation  track,  if  they  exist,  are  masked  by  other 
effects.  Intensity  variations  caused  by  short  duration  vertical  motions  may 
dominate  tlie  more  secular  changes  caused  by  tides. 

A  more  sensitive  test  was  performed  by  summing  and  averaging  all  of 
the  data  in  local  time  to  give  the  single  curve  in  Figure  3.12.  Tliis  is  also  a 
plot  of  tlie  data  as  a  function  of  latitude.  Figure  3.12  shows  tliat  on  average, 
tlie  O2  (0,0)  intensity  changes  by  a  factor  of  1.5  over  a  night  pass.  If  tidal 
effects  are  present,  they  are  masked  by  latitudinal  variations.  The  minimmn 
occurs  near  the  equator,  about  one  hour  after  local  midnight.  There  is  a 
suggestion  of  asymmetry  in  emission  from  south  to  north  (dusk  to  dawn),  with 
the  nortliemmost  (early  morning)  intensity  being  greater.  The  equatorial 
minimum  is  consistent  witli  other  observations,  but  tliere  is  still  no  signature 


of  a  tidal  effect. 
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A  measure  of  tlie  total  energy  emitted  in  the  O2  (0,0)  band  during  a 
night  pass  can  be  determined  by  averaging  the  intensity  in  each  orbit.  This  is 
shown  in  Figure  3.13.  There  is  a  general  upward  trend,  wMi  average  energy 
increasing  by  about  a  factor  of  2  in  tlie  2-day  period.  More  remarkable  is  the 
variability  in  energy  over  short  periods  of  time;  Figure  3.13  shows  a  change  in 
average  intensity  from  30  kR  to  70  kR  in  only  6  orbits  (9  hours  elapsed  solar 
time).  With  this  liigli  level  of  variability,  it  is  not  surprising  that  local 
briglitness  enhancements  change  significantly  from  orbit  to  orbit. 

Since  the  geographic  position  of  each  20  second  spectrograph  exposure 
is  known,  emission  intensity  can  be  examined  as  a  fimction  of  latitude  and 
longitude.  Tlie  observation  pattern  is  shown  in  Figure  3.8;  about  80  exposures 
were  recorded  along  each  orbit  track  in  this  Figure.  Figure  3.14  shows  the 
geograplucal  distribution  of  emission  activity,  but  there  is  no  obvious 
signature  of  organized  wave  activity.  Tlie  average  intensity  along  the 
meridians  is  shown  in  Figure  3.15.  The  intensity  peaks  in  Figure  3.15  appear 
to  have  a  characteristic  width  of  about  25®,  similar  to  the  angular  distance 
between  orbital  tracks  (22.5®)  shown  in  Figure  3.8.  Since  emission  activity  at 
the  equator  is  low,  this  suggested  emissions  might  exhibit  a  zonal  variation 
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moving  from  soutli  to  nordi.  The  intensity  data  were  divided  into  tliree 
latitudinal  zones:  28°  S  to  8°  S  (southern  bin),  8°  S  to  8°  N  (equatorial  bin), 
and  8°  N  to  28°  N  (northern  bin).  Then  the  meridional  averages  were  again 
plotted  in  Figures  3.16a,  3.16b,  and  3.16c  to  examine  the  emission  activity  in 
tliese  zones.  No  periodicity  is  apparent  in  these  Figures.  The  remarkable 
variations  in  local  activity  are  still  evident.  Notice  tliat  in  Figure  3.16a,  tlie 
average  intensity  changes  from  a  minimum  of  18  kR  to  a  maximum  of  80  kR 
in  90°  of  longitude.  Similar  excursions  are  apparent  in  all  three  latitude  bins. 
It  may  be  significant  that  a  deep  minimum  occurs  at  130°  longitude  in  all 
tliree  bms,  but  otlierwise,  tliere  seems  to  be  no  organized  pattern  of  activity. 
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Figure  3.16c.  The  northern  bin:  O2  (0,0)  intensity  averaged  in  longitude  over 
range  of  8°  N  to  28°  N. 
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3.3  The  O2  (0,0)  and  O2  (0,1)  Atmospheric  Bands  Observed  from  Space 

Up  to  this  point,  we  have  utilized  tlie  O2  (0,0)  band  exclusively  in  tliis 
discussion  of  GLO  observations  of  the  night  airglow.  On  the  other  hand,  the 
O2  (0,1)  band  has  characterized  the  O2  Atmospheric  band  system  for  ground- 
based  aeronomers  because  the  O2  (0,0)  band  is  self-absorbed  and  not 
observable  from  the  surface.  However,  since  the  O2  (0,0)  band  is  20.38  tunes 
brighter  than  the  O2  (0,1)  band,  it  is  die  most  intense  emission  seen  from 
orbit.  Although  the  intensity  of  tlie  O2  (0,0)  emission  makes  it  ideal  for  tlie 
work  described  in  the  previous  section,  the  optical  depth  problem  prevents 
tills  band  from  being  useful  for  die  study  of  airglow  chemistry.  Instead,  die  O2 
(0,1)  band  will  be  used  for  diis  purpose. 

On  die  following  pages.  Figures  3.17  and  3.18  show  plots  of  briglitness 
in  himdreds  of  rayleiglis  for  O2  (0,0)  and  (0,1)  during  orbit  72,  in  order  to 
indicate  the  effects  of  optical  depth  on  O2  (0,0).  The  briglitness  of  die  O2 
(0,0)  band  is  scaled  down  by  a  factor  of  20.38,  die  ratio  of  die  transition 
probabilities  of  O2  (0,0)  and  O2  (0,1)  from  die  same  upper  state,  in  order  to 
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compare  tlie  two  emissions.  Discrepancies  between  the  O2  (0,1)  and  scaled 
O2  (0,0)  briglitness  plots  indicate  the  difference  in  the  emission  source 
functions. 

Tlie  O2  (0,1)  emission  isopleths  in  Figure  3.18  is  shown  to  be  sliifled 
upward  compared  to  tlie  O2  (0,0)  isopletlis  in  Figure  3.17,  with  a  maximum  at 
an  altitude  of  about  85  km,  5  km  higher  tlian  tlie  O2  (0,0)  maximum.  Tliis  is 
not  surprising,  since  the  O2  (0,0)  band  originates  from  a  surface  close  to  tlie 
shuttle  orbiter  rather  than  from  the  limb,  where  the  O2  (0,1)  band  originates. 
The  column  brightnesses  for  both  bands  are  very  similar,  but  differences  are 
apparent.  Measurements  of  the  O2  (0,1)  band  intensity  will  be  used  in  die 
discussion  of  atmospheric  chemistry. 

Figures  3.17  and  3.18  show  orbit  72  exhibited  multiple  O2  Atmospheric 
band  maxima  (at  20  LST,  24  LST,  and  2.5  LST),  widi  a  secondary  minimum 
at  21.5  LST  and  a  deeper  minimum  at  1  LST.  Tlie  minimmn  at  1  LST  reflects 
a  feature  common  to  most  night  passes  on  average,  as  depicted  in  Figure 
3.12,  while  the  minimum  at  21.5  LST  is  unique  to  diis  orbit. 
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Figure  3.18.  Brightness  in  hundreds  of  rayleighs  for  the  O2  (0,1)  band  on  orbit  72  as  a 
fijnction  of  altitude  and  local  time. 
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3.4  Spectral  Signatures  of  Night  Airglow  Chemistry 

Tlie  GLO  data  set  is  unrivaled  in  tliat  it  provides  tlie  first  opportumty 
for  an  extensive  examination  of  airglow  emission  characteristics  recorded 
fi-om  a  space  platform.  The  unique  feature  of  GLO  observations  is  the 
capability  to  observe  all  the  related  emissions  fi-om  the  airglow  layer 
simultaneously  fiom  a  single  column  of  gas,  recorded  at  contiguous  altitudes 
across  tlie  Earth’s  limb. 

Figures  3.19a,  3.19b,  and  3.19c  on  the  following  pages  show  the 
averaged  intensity  per  orbit  of  the  O2  (0,0)  band,  5577  A,  and  the  OH  Meinel 
(6,2)  band,  respectively,  as  a  fimction  of  orbit.  Figure  1.19a  has  been  repeated 
fiom  Figure  3.13  for  convenience.  The  O2  (0,0)  band  and  5577  A  (Figures 
3.19a  and  b)  botli  show  a  trend  of  increasing  intensity  witli  successive  orbits 
during  die  period,  wliile  die  OH  Meinel  (6,2)  band  (Figure  3.19c)  does  not. 
Althougli  die  5577  A  and  O2  (0,0)  intensity  plots  are  similar,  5577  A  shows  a 
deep  minimum  at  orbit  91,  but  O2  (0,0)  does  not.  The  5577  A  emission 
increased  by  about  260%  between  orbits  71  and  94,  while  O2  (0,0)  increased 
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by  about  190%.  The  three  emissions  show  no  clear  correlation,  althougli  it 
may  be  significant  that  there  is  a  local  peak  in  all  emission  on  orbit  85. 


AVERAGE  INTENSITY  (rayleighs) 


AVERAGE  INTENSITY  OF  02(0,0) 
AS  A  FUNCTION  OF  ORBIT 
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Figure  3.19a.  Average  intensity  of  O2  (0,0)  for  each  orbit. 
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AVERAGE  INTENSITY  OF  01  (5577  A) 
AS  A  FUNCTION  OF  ORBIT 


ORBIT  NUMBER 


Figure  3.19b.  Average  intensity  of  01  (5577  A)  for  each  orbit. 


AVERAGE  INTENSITY  (rayleighs) 


105 


AVERAGE  INTENSITY  OF  OH  (6,2) 
AS  A  FUNCTION  OF  ORBIT 


ORBIT  NUMBER 


Figure  3.19c.  Average  intensity  of  OH  Meinel  (6,2)  for  each  orbit. 
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Column  intensities  recorded  by  GLO  can  be  inverted  to  determine 
volmne  emission  rate  profiles  for  the  various  species.  Volume  emission  rate  is 
directly  proportional  to  density  of  emitting  species,  and  tlius  is  a  indicator  of 
airglow  chemistry.  The  data  from  three  successive  orbits  (85  tlirougji  87) 
were  inverted  for  tliis  study.  The  inversion  required  application  of  the  Abel 
integral  (Hassani,  1991),  using  an  algoritlun  provided  by  Dasch  (Appendix 
D).  Orbits  85  tlirough  87  are  representative  of  tlie  data  set  and  were  chosen 
for  tlieir  clear  signatures  of  dynamic  behavior  in  the  niglit  airglow.  The 
profiles  for  emissions  of  O2  (0,1),  01  (5577  A),  and  OH  Meinel  (6,2)  during 
orbits  85  througli  87  are  shown  in  Figures  3.20  tlirougli  3.22,  each  Figure 
being  a  set  of  diree  (a,  b,  and  c,  respectively). 

The  shuttle  passed  over  the  South  Atlantic  Anomaly  between  20  and 
22  LST  during  orbit  86  (Figures  3.21a,  b,  and  c).  The  Anomaly  is  a  region  of 
weak  geomagnetic  field  above  die  coast  of  Brazil.  Energetic  charged  particle 
bombardment  of  GLO  spectrograph  B  (see  Table  2.1)  in  tliis  part  of  tlie  orbit 
created  sufficient  noise  to  render  O2  (0,1)  data  useless  until  just  after  22  LST. 
Data  for  OH  Meinel  (6,2),  also  recorded  by  tlie  same  spectrograph,  shows 
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evidence  of  noise  just  before  22  LST,  but  is  not  otherwise  affected.  Data  for 
5577  A  recorded  by  spectrograph  C  was  not  affected  during  this  period. 

Ground-based  observations  over  the  years  have  led  to  the  belief  that 
night  airglow  emissions  are  organized  and  coherent.  Such  orderly  behavior  is 
not  seen  from  orbit.  The  plots  for  orbits  85  through  87,  each  separated  by 
only  90  minutes  in  real  time,  are  certainly  evidence  of  greater  dynamism  in 
die  night  airglow  tlian  is  generally  recognized.  Figures  3.20  through  3.22 
show  tliat  tlie  voliune  emission  rates  for  O2  (0,1),  5577  A,  and  OH  (6,2)  are 
dynamic  and  poorly  organized.  They  exhibit  a  chaotic  pattern  of  emission 
variations  that  changes  significantly  on  time  scales  of  hours  or  less  in  real 
time.  Observations  from  orbit  provide  a  completely  different  view  of  the  night 
airglow  than  ground-based  observations. 

Table  3.1  below  summarizes  the  emission  variations  seen  in  Figures 
3.20  through  3.22.  The  Table  reports  maximum  and  minimum  volmne 
emission  rates  at  the  same  altitude  for  a  given  emission,  as  well  as  tlie  time  of 
occurrence  and  the  maximum-to-minimum  ratio. 
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02(0,1) 

01  (5577  A) 

OH  (6,2) 

Orbit  85 

Maximum 

1300 

22.2  LSI 

750 

2.6  LST 

5500 

23.7  LST 

Minimum 

<200 

0.7  LST 

100 

0.6  LST 

2000 

1.7  LST 

Max/Min  Ratio 

>6.5 

7.5 

2.75 

Orbit  86 

Maximum 

1000 

2.0  LST 

1150 

1.5  LST 

6500 

22.2  LST 

Minimum 

<200 

0.5  LST 

100 

0.3  LST 

2000 

1.2  LST 

Max/Min  Ratio 

>5.0 

11.5 

3.25 

Orbit  87 

Maximum 

900 

2.2  LST 

750 

2.2  LST 

6500 

23.6  LST 

Minimum 

<200 

0.5  LST 

100 

0.2  LST 

2000 

1.3  LST 

Max/Min  Ratio 

>4.5 

7.5 

3.25 

Table  3.1.  Observed  emission  enhancements  during  STS-69  orbits  85  through  87. 


Note  that  although  the  altitudes  for  the  O2  (0,1)  and  OH  Meinel  peak 
emission  only  varied  by  about  1  km  over  the  three  orbits,  tlie  peak  emission  in 
5577  A  rose  from  about  96  km  to  100  km.  Figure  3.23  and  3.24  give 
averaged  vertical  profiles  describing  Regions  I  and  11  as  indicated  on  orbit  86. 
Regions  I  and  II  were  chosen  to  emphasize  that  emissions  are  not  correlated. 
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01  (5577  A)  VOLUME  EMISSION  RATE 
(photons  cm'^  sec'*) 

ORBIT  85 


LOCAL  SOLAR  TIME  (hours) 


Figure  3.20b.  Volume  emission  rate  for  01  (5577  A)  on  orbit  85  as  a  function  of  altitude 
and  local  time. 
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Figure  3.20c.  Volume  emission  rate  for  the  OH  Meinel  (6,2)  band  on  orbit  85  as  a 
function  of  altitude  and  local  time. 
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Figure  3.21a.  Volume  emission  rate  for  the  O2  (0,1)  band  on  orbit  86  as  a  function  of 
altitude  and  local  time. 
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Figure  3.21c.  Volume  emission  rate  for  the  OH  Meinel  (6,2)  band  on  orbit  86  as  a 
function  of  altitude  and  local  time. 
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Figure  3.22a.  Volume  emission  rate  for  the  O2  (0,1)  band  on  orbit  87  as  a  function  of 
altitude  and  local  time. 
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Figure  3.22c.  Volume  emission  rate  for  the  OH  Meinel  (6,2)  band  on  orbit  87  as  a 
function  of  altitude  and  local  time. 
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Figures  3.23  and  3.24  illustrate  that  die  niglit  airglow  emissions  of  O2 
(0,1),  01  (5577  A),  and  OH  (6,2)  are  not  coupled.  Figure  3.23  shows  die 
vertical  profiles  of  volume  emission  rates  averaged  between  0.4  and  0.6  LSI 
on  orbit  86  (Region  I),  about  one  hour  after  the  emission  maximum  for  OH 
Meinel  (6,2),  but  when  O2  (0,1)  and  5577  A  emissions  were  near  dieir 
minimnms  Figure  3.24  shows  the  vertical  profiles  one  hour  later  on  orbit  86 
(Region  II),  averaged  between  1.4  and  1.6  LSI.  Here,  5577  A  emission  was 
at  mavirmim  O2  (0,1)  emission  was  nearly  at  maximum,  and  OH  Meinel  (6,2) 
emission  was  near  minimum.  Tliese  plots  are  not  smooth  because  of  discrete 
sampling  along  the  spectrograph  slit,  but  they  indicate  the  lack  of  temporal 
and  spatial  correlation  among  die  diree  emissions.  Over  an  hour,  5577  A  and 
O2  (0,1)  emissions  varied  by  about  a  factor  of  5,  while  OH  (6,2)  varied  by 
about  a  factor  of  2.  The  emission  maxima  of  01  (5577  A)  is  about  7  to  10  km 
higher  dian  that  of  O2  (0,1),  while  peak  emission  of  O2  (0,1)  is,  in  turn,  about 
7  km  liigjier  than  peak  emission  for  OH  (6,2). 


ALTITUDE  (km) 
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VOLUME  EMISSION  RATES 
AVERAGED  OVER  0.2  HOUR 
ORBIT  86 


Figure  3.23.  Volume  emission  rates  for  airglow  emissions  averaged  over  0.2  hour  between 
0.4  and  0.6  LST  on  orbit  86  (Region  I).  OH  Meinel  (6,2)  emission  was  moderate  while 
5577  A  and  O2  (0,1)  emissions  were  at  or  near  minimum  at  this  time. 
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VOLUME  EMISSION  RATE 
(photons  cm'^  sec'^) 


Figure  3.24.  Volume  emission  rates  for  airglow  emissions  averaged  over  0.2  hour  between 
1.4  and  1.6  LST  on  orbit  86  (Region  H).  OH  Meinel  (6,2)  emission  was  near  minimum 
while  5577  A  and  O2  (0,1)  emissions  were  at  or  near  maximum  at  this  time. 
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3.5  Comparisons  of  GLO  Data  with  the  Johnston-Broadfoot  Model 

The  Johnston-Broadfoot  model,  outlined  in  Appendix  F,  assumes 
statically  stable  conditions.  GLO  data  presented  here  shows  this  assumption 
is  far  from  reality.  However,  tlie  model  can  provide  a  baseline  against  which 
observations  can  be  compared. 

Sensitivity  studies  using  the  Johnston-Broadfoot  model  indicate  that 
emissions  at  5577  A  on  the  one  hand,  and  emissions  in  the  OH  Meinel  and  O2 
Atmospheric  bands  on  the  other  hand,  react  very  differently  to  changes  in 
ambient  number  density  [M]  (=  [N2]  +  [O2])  and  oxygen  density  [O]  in  the 
airglow  layer.  By  increasing  ambient  number  density  [M]  in  tlie  model  (as 
might  be  caused  by  upwelling  from  below),  emissions  at  5577  A  decrease  as 
a  result  of  enhanced  quenching  of  [O*],  wliile  emissions  in  the  OH  Meinel 
and  O2  Atmospheric  bands  increase  as  a  result  of  enlianced  production. 
Wlien  O  concentration  alone  is  increased  (as  might  be  caused  by  subsidence 
from  above),  emissions  at  5577  A  and  in  the  O2  Atmospheric  bands  increase, 
while  OH  Meinel  band  emissions  remain  ahnost  unchanged,  exhibiting  only  a 
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sliglit  increase.  Vertical  transport  is  assumed  to  be  important  in  the  creation  of 
tlie  localized  emission  activity  detected  in  tlie  GLO  data  set.  Upwelling  and 
subsidence  effects  are  examined  here  with  respect  to  the  predictions  of  tlie 
Johnston-Broadfoot  model. 

The  [O]  and  pvi]  inputs  to  the  Jolinston-Broadfoot  model  can  be  varied 
as  desired.  Simple  upward  sliifl  of  [M]  values  in  the  input  data  array  witli 
respect  to  altitude  simulates  an  upwelling.  Increasing  [O]  by  a  constant  factor 
simulates  the  accumulation  of  subsiding  O  from  above.  Varying  the  input 
profiles  of  eddy  viscosity  above  90  km  simulates  atmospheric  turbulence  and 
gravity  wave  “breaking”  (turbulent  dissipation),  which  is  presmned  to  be 
ubiquitous  at  altitudes  near  100  km  (the  so-called  “turbopause”). 

Table  3.2  summarizes  the  results  of  varying  key  parameters  in  tlie 
Jolmston-Broadfoot  model.  The  Table  indicates  that  only  upwelling  in  a 
relatively  shallow  layer  enliances  O2*  emissions,  wliile  upwelling  in  general 
increases  OH*  emissions  while  decreasing  O*  emissions.  Increased  [O] 
enliances  all  emissions.  Molecular  viscosity  changes,  even  by  a  factor  of  10 
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(whetlier  increases  or  decreases),  do  not  seem  to  sigmficantly  affect  any 
emissions.  Increasing  eddy  viscosity  enliances  both  O2*  mid  OH*  emissions 
while  leaving  O*  emissions  unchanged.  Decreasing  eddy  viscosity  by  a  factor 
of  10  produces  negligible  effects  for  all  emissions. 


Parameter 

Change  in 
Parameter 

Change  in 
[O2*] 

Change  in 
[OH*] 

Change  in 
[0*] 

[M] 

5  km  uplift 

X  1.35 

X  1.63 

xO.65 

[M] 

10  km  uplift 

Negligible 

x2.15 

xO.24 

[M] 

20  km  uplift 

x0.17 

Negligible 

xO.Ol 

[0] 

x2 

x3.17 

X  1.67 

x7.84 

[O] 

X  3 

X  5.98 

x2.11 

x24.7 

[0] 

x4 

x9.20 

x2.45 

X  54.8 

[M]  &  [0] 

[M]  x2& 
roix2 

x2.58 

X  1.86 

x2.80 

[M]  &  [0] 

[M]  x2& 

roix3 

x7.94 

x3.63 

X  18.3 

Eddy  Viscosity 

X  10 

X  1.17 

X  1.20 

Negligible 

Eddy  Viscosity 

x20 

X  1.26 

X  1.34 

Negligible 

X  0.1 

Negligible 

Negligible 

Molecular 

Viscosity 

X  10  or  X  0.1 

Negligible 

Negligible 

Negligible 

Molecular 

Viscosity 

5  km  uplift 

Negligible 

Negligible 

Negligible 

Table  3.2.  Results  of  varying  input  parameters  in  the  Johnston-Broadfoot  model. 


An  increase  in  [O]  by  a  factor  of  2  in  the  Johnston-Broadfoot  model 
changes  [O*]  and  [OH*]  by  factors  of  7.84  and  1.67,  respectively.  Tliis 
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closely  parallels  tlie  emission  changes  by  factors  of  7.5  and  2.75  actually 
observed  for  5577  A  and  OH  (6,2)  on  orbit  86,  as  shown  in  Table  3.1. 
However,  tlie  model  indicates  that  increasing  [O]  by  a  factor  of  2  is  not 
sufficient  to  raise  [O2*]  by  die  factor  of  6.5  seen  on  orbit  86,  suggesting  that 
upwelling  [M]  as  well  as  subsiding  [O]  caused  die  O2  (0,1)  emission  peak. 
The  model  predicts  diat  upwelling  decreases  [O*].  The  enhancements  on  orbit 
87  by  factors  of  5.0,  11.5,  and  3.25  for  O2  (0,1),  5577  A,  and  OH  (6,2) 
emissions  are  similar  to  model  output  factors  of  5.98,  24.7,  and  2.11  for  an 
[O]  increase  by  a  factor  of  3.  However,  a  blend  of  upwelling  and  subsidence 
is  again  suggested  because  of  die  discrepancy  between  predicted  and  actual 
5577  A  enhancement.  Model  output  for  an  [O]  increase  by  a  factor  of  2  again 
closely  matches  die  actual  observations  on  orbit  87,  and  more  closely 
matches  the  smaller  increase  in  O2  (0,1)  emission  on  this  orbit.  Since 
upwelling  from  below  and  subsidence  from  above  need  not  be  mutually 
exclusive  events,  we  can  conclude  that  the  two  effects  could  compete  to 
detennine  die  actual  emission  enliancements  observed  by  GLO.  Tliis  is 
consistent  with  die  scale  of  local  dynamics  inferred  during  die  search  for 
organization  on  a  more  global  scale  in  Section  3.2. 
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A  separate  test  for  temperature  sensitivity  in  the  Johnston-Broadfoot 
model  was  performed.  Using  an  extreme  temperature  increase  of  40  K  at 
peak  in  a  vertically  symmetric  layer,  changes  were  negligible  and  nearly  all 
species  emissions  were  sliglitly  decreased,  wMi  only  tlie  OH  Meinel  bands 
seemingly  unaffected  by  tlie  temperature  change.  Model  output  does  not 
appear  to  support  the  contention  that  temperature  variations  in  situ  can 
directly  force  die  brightness  changes  observed.  Model  output,  however, 
cannot  simulate  die  secondary  effect  of  temperature  variations:  namely,  die 
destruction  of  static  stability  and  the  creation  of  vertical  motions,  which 
undoubtedly  will  change  [M]  and  [O]. 

The  production  of  O2*  (excited  state  molecular  oxygen)  is  detennined 
by  die  diree  body  association  reaction  O  +0  +  M  — >  O2*  +  M  (Step  I  of  the 
Barth  process  described  in  Chapter  1).  We  see  diat  O2  emissions,  being 
proportional  to  die  number  density  [O2*],  should  vary  as  [M][0]^,  and  are 
strongly  dependent  on  die  number  density  of  ground  state  atomic  oxygen 
(Shepherd  et  al.,  1995).  Widi  O2  emissions  (and  thus  [O2*])  varying  by 
factors  of  diree  or  more,  we  see  that  [O]  must  vary  by  a  factor  of  about  two  to 
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explain  our  observations  during  STS-69,  assuming  PVI]  does  not  change. 
Similarly,  since  tlie  production  reaction  for  excited  state  atomic  oxygen  is 
Step  n  of  die  Bardi  process,  O2*  +  O  ->  O2  +0*,  we  see  tiiat  variations  in 
O*  (neglecting  quenching  effects)  must  follow  the  product  [02*][0],  or 
alternatively  [M][0]^.  This  product  gives  an  enhancement  factor  of  6,  if  [O2*] 
increases  by  a  factor  of  3  and  [O]  increases  by  2,  This  is  in  good  agreement 
witii  die  variations  of  O*  emissions  actually  observed  in  die  GLO  data  set. 

Quenching  caused  by  die  increases  in  ambient  density  [M]  and  atomic 
oxygen  [O]  must  be  taken  into  account.  02*  is  primarily  quenched  by 
collisions  with  N2,  which  has  a  concentration  rouglily  equal  to  0.78  [M], 
while  O*  is  quenched  both  by  O  and  O2.  The  quenching  rate  coefiticient  for 
O2*  and  N2  is  2  X  10’^^  cm^  s'^;  for  O*  and  O2,  3  x  10'^^  cm^  s'*;  and  for  O* 
and  O,  7.5  x  10'*^  cm^  s'*  (Vallance  Jones,  1974).  Altiiough  O  is  25  times 
more  efficient  tiian  O2  at  quenching  O*,  O2  is  so  much  more  numerous  than  O 
that  O*  losses  caused  by  0*-02  collisions  are  40  times  greater  tiian  die 
losses  caused  by  0*-0  collisions.  Upwelling  motions  will  increase  [O2]  and 
[N2],  enliancing  both  production  and  quencliing  of  O2*  and  OH*  by  [M],  but 
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increased  production  dominates  tlie  increased  quenching  in  this  case. 
However,  O*  quenching  by  [M]  will  dominate  O*  production  whenever  [M] 
is  increased,  as,  for  example,  in  upwelling  from  lower  altitudes.  Thus, 
chemistry  considerations  indicate  that  GLO  data  are  consistent  with  vertical 
transport  of  reactants  into  the  airglow  layer. 

A  measxire  of  die  severity  of  quencliing  is  the  altitude  at  which  an 
emission  is  reduced  to  half  its  unquenched  rate.  This  parameter  is  known  as 
the  quenching  height,  hq.  Above  hq,  quenching  quickly  becomes  neghgible, 
while  below  hq,  quencliing  is  extreme.  The  greater  hq  for  a  given  emission, 
the  greater  that  emission’s  sensitivity  to  quencliing. 

For  example,  the  atomic  oxygen  line  at  6300  A  has  a  quenching  height 
between  250  and  350  km,  and  it  is  not  observed  at  airglow  layer  altitudes.  For 
die  5577  A  emission,  hq  ranges  from  95  to  115  km.  For  die  O2  (0,0)  emission, 
hq  is  about  95  km  (Vallance  Jones,  1974).  The  tendency  here  is  for  O*  to  be 
more  sensitive  to  quenching  than  O2*.  Tlie  hq  for  the  OH  Meinel  band  is  not 
well  established.  Figures  3.20  dirougli  3.24  indicate  OH  emission  is  reduced 
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when  5577  A  emission  increases,  possibly  as  a  result  of  subsidence  of  O. 
Tin’s  would  imply  OH  quenching  by  O  widi  die  reaction  0  +  OH  ^  O2  +  H. 
Thus,  hq  for  OH  Meinel  emissions  could  be  quite  variable,  depending  on  the 
dynamics  of  die  downward  flow  of  O  into  die  MTI  region  at  night. 
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3.6  Comparisons  of  GLO  Data  witli  Previous  Observations 

GLO  data  from  STS-69  agrees  in  general  with  observations  reported  in 
both  liistorical  and  contemporary  research.  The  uniqueness  of  the  GLO  data 
set  provides  insight  into  the  nature  of  niglit  airglow  chemistry  and  dynamics. 
In  particular,  GLO  observed  evidence  for  tlie  dynamic  influence  of  vertical 
winds  on  niglit  airglow  emissions. 

Airglow  parameters  based  on  previous  research,  the  Jolmston- 
Broadfoot  model,  and  tlie  inversion  of  GLO  data  (Appendix  D)  are  shown  in 
Table  3.3.  Parameters  based  on  rocket  observations  are  marked  by  tlie  symbol 
“A  ”  orbital  observations  and  model  output  references  are  unmarked. 
The  volume  emission  rate  for  the  O2  (0,0)  band  marked  witli  the  symbol 
was  estimated  by  multiplying  the  rate  for  the  O2  (0,1)  band  by  20.  The  results 
of  Shepherd  et  al.  and  Ward  et  al.  referenced  in  Table  3.3  come  from  data 
recorded  by  die  WINDII  instrument  on  board  UARS.  Observations  from 
orbit,  such  as  those  by  GLO  and  WINDII,  can  capture  airglow  variations  on  a 
planetary  scale,  unlike  spatially  isolated  and  short  duration  rocket  data. 
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EMISSION 

FEATURE 

PEAK 

ALT 

km 

FWHM 

km 

SOURCE 

No  estimate 

89-95 

10 

Ward  et  al.,  1995 

2400-11000 

88-96 

5-12 

^Lopez-Gonzalez  et  al.,  1992a 

4800 

94 

13 

^Murtagh,  1995 

8000-26000* 

89-93 

11-13 

This  dissertation,  1997 

1600 

94 

10 

Johnston-Broadfoot  Model,  1992 

No  estimate 

88-96 

5-12 

'^Lopez-Gonzalez  et  al.,  1992a 

400-1300 

89-93 

11-13 

This  dissertation,  1997 

78 

94 

10 

Johnston-Broadfoot  Model,  1992 

5577  A 

100  -  400 

92-97 

5-7 

^Lopez-Gonzalez  et  al.,  1992a 

165 

97 

8 

^Murtagh,  1995 

140 

95 

8 

Shepherd  et  al.,  1995a 

370 

92 

12 

Ward  et  al.,  1995 

65-110 

95-98 

7-13 

^Kitaetal.,  1992 

250-1150 

96-100 

10-12 

This  dissertation,  1997 

78 

96 

7 

Johnston-Broadfoot  Model,  1992 

2500  -  6500 

85-87 

10-12 

This  dissertation,  1997 

■HlHi 

750 

89 

10 

Johnston-Broadfoot  Model,  1992 

Table  3.3.  Comparison  of  measured  and  modeled  airglow  parameters. 


UARS  researchers  (Burrage  et  al.,  1994)  have  indicated  tliat  they 
detected  tidal  signatures  in  their  data  averaged  over  periods  of  weeks  and 
months.  However,  they  reported  a  persistent  region  of  maxunmn  emission 
between  10°  N  and  10°  S  tliat  ran  from  sunset  to  23  LSI,  and  a  region  of 
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minimum  emission  beginning  just  after  midnight.  This  is  just  west  of  the 
observation  track  shovm  in  Figure  3.8.  Our  data  set  is  consistent  with  tliis 
observation  of  Burrage  et  al,  since  tlie  evening  tidal  enhancement  they 
describe  could  not  have  been  detected  by  GLO.  However,  Burrage  et  al.  also 
detected  a  short  term  enhancement  at  tlie  equator  that  persisted  for  a  few 
days.  GLO  did  not  observe  any  such  long-lived  region  of  organized  emission 
activity,  at  the  equator  or  elsewhere. 

Retrieved  altitudes  of  peak  emission  are  consistent  with  output  from 
the  Johnston-Broadfoot  model  and  previous  measurements.  The  Jolmston- 
Broadfoot  model  depends  on  an  Air  Force  liigh  altitude  winter  atmospheric 
model  that  may  not  be  well  suited  for  low  to  midlatitudes  where  the  STS-69 
data  set  was  recorded.  From  die  results  of  Table  3.3,  the  Johnston-Broadfoot 
model  appears  to  be  systematically  underestimating  reactant  concentrations. 
The  differences  between  model  predictions  and  GLO  observations  indicate 
tliat  modeled  reactant  concentrations  are  low  by  about  a  factor  of  5.  From 
GLO  data,  O2  Atmospheric  emissions  were  seen  to  peak  at  around  90  km, 
OH  Meinel  at  around  85  km,  and  01  (5577  A)  at  around  95  km,  altliough 
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these  altitudes  were  variable.  Fluctuations  in  tliese  altitudes  were  noted  in 
Figures  3.20  through  3.24  in  section  3.4. 

Swenson  et  al.  (1989)  observed  niglitly  altitude  shifts  of  about  8  km  in 
tlie  5577  A  and  O2  night  airglow  layers  using  the  Atmospheric  Emissions 
Photometric  Imager  (AEPI)  on  STS-9.  AEPI  observations  were  consistent 
with  earlier  findings  by  Wasser  and  Donaliue  (1979)  based  on  OGO  6  data. 
Similar  variations  have  been  confirmed  by  WINDII  and  now  by  GLO,  with 
variations  of  7  to  10  km  observed  on  orbits  85  through  87. 

Shepherd  et  al.  (1993)  have  observed  emission  in  5577  A  to  vary  by 
factors  of  4  or  more,  with  similar  patterns  exhibited  by  OH  emissions.  Ward 
et  al.  (1995)  observed  the  5577  A  emission  to  vary  by  a  factor  of  10  and  die 
O2  (0,0)  emission  by  a  factor  of  5.  These  findings  are  in  good  agreement  with 
GLO  data  that  indicate  enhancement  factors  between  3  and  12,  in  general,  for 
all  emissions.  They  are  also  in  accord  witli  the  Jolinston-Broadfoot  model. 
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Ward  et  al.  (1995)  have  inferred  from  tlieir  WINDII  observations  tliat 
emissions  in  01  (5577  A)  and  the  O2  Atmospheric  (0,0)  band  are  highly 
correlated.  GLO  observations  suggest  that  die  two  emissions,  when  averaged 
over  time  scales  of  days,  do  indeed  show  similar  trends  in  enhancement. 
However,  on  time  scales  of  hours  or  minutes,  diese  emissions  are  not 
coincident  in  space  and  time,  and  they  are  neither  simply  correlated  nor 
anticorrelated.  Chemistry  considerations  and  output  from  die  Johnston- 
Broadfoot  model  imply  that  the  dynamic  interaction  of  upwelling  from  below 
and  subsiding  atomic  oxygen  from  aloft  probably  drives  diese  briglitness 
changes  in  a  complex  fashion.  Upwellings  may  be  connected  to  upwardly 
propagating  gravity  waves.  The  associated  increase  in  eddy  viscosity  with 
such  propagating  gravity  waves  would  tend  to  increase  OH  and  O2  (0,0) 
emissions  wliile  leaving  those  of  5577  A  unchanged.  The  downward  flow  of 
atomic  oxygen  in  the  MTI  region  is  only  now  being  parameterized  in  general 
circulation  models  at  die  National  Center  for  Atmospheric  Research  (NCAR) 
and  is  not  yet  well  understood. 
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Roble  and  Shepherd  (1997)  describe  one  NCAR  general  circulation 
model  in  wliich  strong  upwardly  propagating  tides  disrupt  tlie  5577  A  airglow 
layer  near  95  km  and  reduce  its  emission  rate.  Tliis  same  upwelling  would 
increase  O2  (0,0)  briglitness  and  bring  its  peak  emission  rate  to  a  lower 
altitude.  Such  effects  are  consistent  with  GLO  observations.  Roble  and 
Shepherd  also  describe  a  weaker  tide  which  does  not  penetrate  high  enougli 
to  quench  tlie  5577  A  airglow  layer,  and  tliis  may  be  more  consistent  witli  the 
correlation  assumption  of  Ward  et  al.  (1995)  based  on  their  WINDII 
observations.  If  strong  tides  are  presumed  to  be  intennittent  and  of  brief 
duration,  averaging  over  long  time  scales  as  practiced  by  UARS  researchers 
may  reduce  tlieir  signatures  in  airglow  emission  data. 

Shepherd  et  al.  (1993)  reported  evidence  for  planetary  wave  influence 
on  5577  A  emissions  from  analysis  of  data  recorded  by  WINDII  on  UARS  in 
January  1992.  They  detected  possible  diurnal  and  semidiurnal  wave  patterns 
in  tliese  emissions  in  tlie  narrow  latitude  band  of  40°  to  42°  N.  Their 
conclusions  based  on  this  restricted  coverage  for  a  single  emitting  species  do 
not  agree  witli  tlie  phenomena  observed  by  GLO  on  a  global  scale  during  die 
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mission  of  STS-69.  Shepherd  et  al.  also  used  the  WINDII  experiment  to  infer 
die  effect  of  atmospheric  tides  on  5577  A  emission  in  the  equatorial  region. 
They  observed  tliat  5577  A  brightness  in  the  niglit  airglow  layer  near  100  km 
decreased  by  a  factor  of  2.7  in  a  6  hour  period  from  local  dusk  to  midniglit. 
GLO  typically  observed  variations  for  5577  A  much  greater  than  this,  and  on 
orbit  86  recorded  an  increase  by  a  factor  of  nearly  12  in  a  one-hour  period. 
Shepherd  et  al.  inferred  from  the  UARS  data  that  the  5577  A  decrease  by  a 
factor  of  2.7  could  caused  by  a  decrease  in  [M]  by  a  factor  of  1.7,  combined 
witli  a  decrease  in  [O]  by  a  factor  of  1.3.  They  reported  a  subsequent  gradual 
briglitening  at  5577  A  toward  morning  twilight,  ascribing  this  to  subsidence 
of  atomic  oxygen-rich  air  from  above  after  local  midnight,  as  predicted  by  the 
NCAR  Thenuosphere/Ionosphere  General  Circulation  Model  (TIGCM)  of 
Roble  (Roble  and  Shepherd,  1997).  The  briglitening  of  tlie  5577  A  emission 
in  response  to  the  subsidence  of  [O]  is  consistent  witli  die  GLO  observations 
summarized  in  Table  3.1,  as  discussed  previously. 
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4.  CONCLUSIONS 

GLO  recorded  midlatitude  night  airglow  data  near  equinox  during 
Space  Shuttle  mission  STS-69.  GLO  observations  from  tlie  shuttle  provided  a 
record  of  the  night  airglow  layer  seen  edge  on  at  the  Earth’s  limb.  These  top 
side  observations  show  a  fundamentally  different  picture  of  the  night  airglow 
compared  to  observations  from  the  ground.  They  complement  and  extend  the 
phenomenology  generated  by  otlier  global  observations  from  orbit.  GLO 
observations  also  represent  the  first  simultaneous  optical  measurements  of 
airglow  emissions  over  die  spectral  range  from  1150  to  9000  A,  showing 
global  emission  variations  in  the  night  sky.  When  primary  emissions  are  all 
measured  simultaneously,  intensity  variations  are  shown  to  be  uncorrelated 
among  emitting  species,  implying  greater  dynamism  and  more  complex 
chemical  interactions  in  the  airglow  tlian  previously  assumed. 

Altliough  otlier  researchers  describe  observations  of  organized  waves 
and  tides  in  the  night  airglow,  those  observations  are  averaged  over  long  time 
periods.  The  night  sky  intensity  variations  recorded  by  GLO  do  not  exhibit 
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any  obvious  relationship  to  atmospheric  tides.  They  may  instead  be 
associated  with  tlie  chaotic  superposition  of  upward  and  downward  vertical 
motions. 

Intensity  variations  in  the  niglit  airgjow  layer  are  hkely  to  be  driven  by 
upwelling  and  subsidence.  These  vertical  motions  may  cause  the  decoupling 
of  emissions  in  the  O2  Atmospheric  band,  die  OH  Meinel  band,  and  the 
atomic  oxygen  green  line  at  5577  A  in  die  night  airglow.  Emission 
enliancements  with  maxima-to-minima  ratios  of  4  to  12,  depending  on 
emitting  species,  have  been  observed.  Inferred  horizontal  extent  of  observed 
emission  enhancement  regions  is  on  the  order  of  3000  km.  Emissions  in  the 
O2  Atmospheric  band  system  and  die  01  (5577  A)  green  line  show  a  greater 
dynamic  range  of  variation  than  die  OH  Meinel  band  system. 

Briglitening  in  die  O2  and  OH  bands  and  dimming  of  the  01  (5577  A) 
emission  may  coincide  widi  die  lowering  of  the  altitude  of  peak  emission. 
Tliis  signature  is  consistent  widi  upwelling  from  below.  The  upwelling 
increases  die  concentration  of  reactants  for  die  production  of  O2*  and  OH*, 
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wliile  quenching  O*,  and  brings  peak  production  of  all  tliree  emitters  to  lower 
altitudes.  Brightening  of  5577  A  (and  occasionally  in  O2  and  OH  emissions) 
occurs  with  the  subsidence  of  air  from  above,  enriched  witli  O  (compared  to 
otlier  species).  Lifting  of  peak  emission  altitudes  may  be  associated  with  diis 
subsidence. 

The  chemistry  along  a  limited  line-of-siglit  can  be  explained  by 
classical  airglow  chemistry  as  modeled  by  Jolinston  and  Broadfoot,  but  only 
over  a  limited  altitude  range.  Dynamic  effects  in  the  80  to  100  km  MTI  region 
are  sufiBciently  chaotic  to  present  mixed  results  when  inferring  chemical 
processes  as  a  function  of  altitude. 

The  GLO  instrument  is  an  invaluable  resource  for  investigation  of  die 
mesosphere-thermosphere-ionosphere  (MTI)  region.  Tlie  ability  to  record 
atmospheric  phenomena  simultaneously  over  a  wide  spectral  range  on  a 
global  scale  is  critical  to  the  advancement  of  our  knowledge  of  the  night 
airglow.  GLO  can  therefore  make  a  significant  contribution  to  ongoing 
CEDAR  campaigns  designed  to  study  tlie  middle  atmosphere. 
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APPENDIX  A:  DESCRIPTION  OF  THE  GLO  INSTRUMENT 

The  Arizona  Airglow  Instrument,  or  GLO,  is  designed  to  study  auroral 
and  airglow  emissions.  From  the  Shuttle,  GLO  can  measure  spatial,  temporal, 
and  spectral  emission  signatures  at  wavelengths  from  1150  to  11000  A, 
typically  for  altitudes  ranging  from  about  80  to  350  km.  The  following 
description  refers  to  a  nominal  configuration;  tliis  may  be  (and  has  been) 
altered  depending  upon  the  needs  of  a  particular  mission.  GLO  consists  of 
nine  spectrographs  and  twelve  unagers  with  coaligned  optical  axes, 
incorporating  botli  unintensified  (“bare”)  charge-coupled  devices  (CCDs)  and 
intensified  charge-coupled  devices  (ICCDs).  As  a  result  of  optimal  design  to 
minimize  weight  and  volume,  the  spectrographs  share  four  image  intensifiers 
wliile  tlie  imagers  share  tliree  intensifiers  (Broadfoot  and  Sandel,  1992). 

The  CCDs  used  in  tlie  GLO  instrument  are  384  by  576  pkel  two- 
dhnensional  arrays.  Each  CCD  is  divided  into  two  equal  sections  measuring 
192  by  576  pixels.  Each  pixel  is  a  0.022  mm  square  and  the  overall  sensitive 
area  of  tlie  CCD  is  12.7  by  8.4  mm.  Spectral  dispersion  along  the  lengtli  of 
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tlie  CCD  uses  all  576  pixels,  but  the  heiglit  of  a  single  spectrum  only  requires 
8  pixels  of  the  CCD.  Tlius  two  sets  of  24  spectra  each  can  be  focused  side  by 
side  onto  a  single  detector,  retaining  spectral  resolution  wliile  using  only  192 
pixels  per  set  of  24  (8  pixels  per  spectrum  x  24  spectra  =192  pixels).  This 
reduces  the  size  of  the  spectrographs  and  uses  the  CCDs  more  efficiently.  Tlie 
detectors  operate  effectively  only  below  0®  C.  This  is  not  a  problem  for  tlie 
typical  mission  on  orbit,  but  a  cooling  system  can  be  applied  if  mission 
requirements  or  ground  laboratory  tests  require  it.  Active  cooling  can  be 
accomplished  by  tliermoelectric  or  circulating  refrigerant  techniques,  wliile 
passive  cooling  with  radiators  is  also  possible. 

Tlie  seven  intensifiers  used  in  GLO  each  consist  of  a  photocathode  to 
detect  incoming  photons  and  a  phosphor  screen  to  generate  tlie  image. 
Between  these  components,  a  series  of  microchannel  plate  electron 
multipliers  amplify  the  signal  to  produce  a  brigliter  image  on  die  phosphor 
screen.  The  voltage  across  these  microchannel  plates  determines  tlie  gam  of 
tlie  intensifier.  A  single  photon  entering  the  intensifier  and  hitting  die 
photocadiode  can  cause  10"^  photons  in  die  5000  to  7000  A  wavelength  region 
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to  be  emitted  from  the  phosphor  screen.  Exposure  is  controlled  by  setting  tlie 
voltage  between  tlie  photocathode  and  the  nearest  microchannel  plate.  Wlien 
tlie  potential  of  tlie  photocathode  is  +40V  with  respect  to  die  front 
microchannel  plate,  no  electrons  are  emitted  by  the  cadiode.  With  the  cathode 
at  -200V,  photoelectrons  are  accelerated  and  colhde  with  the  front 
microchannel  plate,  mitiating  die  electron  multiplication  cascade.  Each  photon 
interaction,  or  photoevent,  is  spatially  indexed  by  the  microchannel  plates  and 
die  fiberoptics  coupler  to  die  associated  CCD.  Since  the  spectral  response  of 
an  intensifier  depends  on  die  material  of  die  photocathode,  three  types  of 
photocatiiodes  are  used  to  fix  the  sensitivity  of  the  intensifiers  to  specific 
spectral  regions  of  interest  (Broadfoot  and  Sandel,  1992). 

Image  intensifiers  make  each  pixel  of  an  ICCD  fimction  as  a  photon 
counting  device  that  preserves  die  Poisson  (discrete)  statistics  of  arriving 
photons.  The  dynamic  range  of  an  ICCD  is  not  limited  by  the  number  of  free 
electrons  per  pixel  in  die  detector  (the  “full-well  capacity”),  but  rather  by  die 
dark  count  rate  and  available  digital  memory.  The  dark  count  rate  for  the 
ICCDs  used  in  GLO  is  quite  small,  at  about  10  ^  events  per  pixel  per  second. 
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At  m?^xim1lm  gain,  tlie  image  intensifier  generates  about  500  electrons 
per  pixel  per  pliotoevent.  A  typical  M-well  capacity  of  about  2.5  x  10^  free 
electrons  implies  a  dynamic  range  limit  of  about  500  photoevents  per  pixel. 
However,  ICCD  counts  can  be  read  out  frequently  and  quickly  accumulated 
in  memory.  Integrating  the  signal  in  this  way  involves  a  modest  statistical 
penalty  factor  of  2^^  in  going  from  pulse  counting  to  pulse  integration,  but  the 

7 

dynamic  range  can  be  dramatically  extended  to  as  high  as  10  photoevents  per 
pixel,  witli  a  detection  tlireshold  of  about  10'^  photoevents  per  pixel  per 

4 

second.  Compare  this  witli  a  dynamic  range  of  about  10  photoevents  per 
pixel  and  a  detection  threshold  of  about  50  photoevents  per  pixel  per  second 
for  a  bare  CCD!  Tire  photocatliode  of  an  ICCD  can  also  be  chosen  to 
eliminate  sensitivity  to  unwanted  wavelengtlis,  an  additional  advantage  over 
tlie  bare  CCD,  which  can  yield  spurious  signals  after  being  contaminated  by 
photons  outside  the  wavelengtli  region  of  interest  (Broadfoot  and  Sandel, 
1992). 


Complete  spectral  coverage  from  die  UV  to  tlie  near-IR  eliminates  the 
ambiguities  of  partial  spectral  coverage,  and  recording  the  entire  spectnim 
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simultaneously  allows  the  measurement  of  temporal  variations.  Resolution  at 
tlie  shortest  wavelengtlis  is  4  A,  while  resolution  is  10  A  at  tlie  longest 
wavelengths.  Eiglit  of  the  nine  spectrographs  use  ICCD’s,  with  the  nintli 
using  a  bare  CCD.  The  twelve  imagers,  fitted  with  interference  filters,  have 
varying  fields  of  view  and  spectral  sensitivities  in  order  to  measure  the  spatial 
extent  of  briglitness  variations  over  several  orders  of  magnitude. 

Tlie  brightness  of  the  target  image  on  the  detector  depends  only  on  the 
f-ratio  of  tlie  optics.  The  instrument’s  CCD  focal  plane  detectors  are  very 
small,  permitting  much  shorter  optical  patli  lengtlis  tlian  other  technologies 
while  maintaining  tlie  desired  f-ratios  and  sensitivities.  GLO  uses  concave 
holograpliic  gratings  that  correct  spherical  and  chromatic  aberrations,  as  well 
as  replace  die  collimating  lens,  plane  grating,  and  refocusing  lens  needed  in  a 
conventional  spectrograph.  These  factors  combine  to  reduce  the  size  and 
weiglit  of  tlie  instrument.  The  sensor  head  includes  all  optics,  the  CCDs  and 
ICCDs,  and  some  of  the  electronics  for  GLO;  it  measures  about  40  by  40  by 
30  cm  and  weiglis  less  tlian  35  kg. 
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The  nine  nearly  identical  spectrographs  can  each  observe  a  long  thin 
swath  of  tlie  target  region.  Eiglit  spectrographs  are  grouped  in  four  pairs  (pair 
#1  includes  two  UV,  pair  #2  includes  one  visible  and  one  near-UV,  pair  #3 
includes  two  visible,  and  pair  #4  includes  one  visible  and  one  near-IR 
spectrograph).  Each  pair  shares  a  sht,  foreoptics,  and  an  optical  bench 
mounting  surface.  Tlie  foreoptics  for  the  near-IR,  near-UV,  and  visible  liglit 
optical  benches  consist  of  a  cover  with  hgjit  baffles,  two  cylindrical  focusing 
lenses,  a  spectral  cutoff  filter  to  eluninate  higher  order  spectra,  and  two  path¬ 
folding  prisms.  Uie  UV  foreoptics  are  similar,  but  utilize  cylindrical  mirrors 
instead  of  lenses  to  maintain  focus.  Each  slit  has  dimensions  0.045  by  4.5 
mm.  This  arrangement  gives  side  by  side  image  sets  at  tlie  CCD,  with  each 
slit  image  set  measuring  2  pixels  wide  and  192  pixels  long.  (Recall  that  each 
line  of  a  spectrum  is  an  image  of  the  slit).  The  192  pixels  are  divided  into  24 
segments  of  8  pixels,  and  a  separate  spectrum  is  recorded  in  each  segment. 
The  ninth  spectrograph,  operating  in  the  IR,  has  its  own  sht,  foreoptics,  and 
optical  bench.  The  IR  spectrograph  foreoptics  include  a  negative  focal  length 
lens,  and  this  optical  bench  has  a  slit  with  dimensions  0.045  by  8  mm.  For  all 
nine  spectrographs,  tlie  distances  between  die  gratings  and  die  slits  are  die 
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same,  and  tlie  fields-of-view  (FOV)  determined  by  the  foreoptics  are  about 
8.5®  along  die  slit  and  0.15®  perpendicular  to  the  slit.  Each  grating  was 
designed  for  a  specific  spectrograph  to  give  dispersion  and  spectral  coverage 
to  match  the  size  of  the  detector  (Viereck  et  al.,  1990). 

On  the  horizon,  at  typical  ranges  of  about  1500  km  fi-om  a  Shuttle 
orbiting  at  about  340  km,  the  field-of-view  of  a  spectrograph  encompasses 
about  a  260  by  4  km  swath  of  die  atmosphere.  With  die  sht  tangent  to  die 
limb,  a  4  km  depdi  of  atmosphere  is  dius  observed  by  the  instrument.  If  the 
slit  is  tilted  at  an  angle  to  die  limb,  a  deeper  scan  of  the  atmosphere  is 
possible.  For  example,  at  an  angle  of  10®  to  die  horizontal,  a  45  km  depth  of 
die  atmosphere  is  sampled,  while  at  90°,  a  maximum  altitude  range  of  about 
260  km  can  be  imaged  on  the  slit.  Since  the  slit  image  is  divided  into  24 
segments,  spectra  from  24  contiguous  atmospheric  layers,  each  of  about  11 
km  depth,  may  be  obtained  simultaneously  in  the  latter  case. 

Each  imager  system  consists  of  lenses  to  focus  die  source  onto  an 
image  intensifier,  an  interference  filter  to  limit  spectral  response,  an  imager 
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intensifier  to  increase  sensitivity,  and  a  CCD  to  collect  die  images  and 
transform  diem  into  electronic  signals.  (Note,  however,  diat  image  intensifiers 
are  not  used  for  die  infrared  imagers).  Tlie  spectrographs  and  imagers  both 
use  flexible  fiberoptic  cable  bundles,  or  “image  conduits”,  to  carry  images 
from  die  intensifiers  to  die  CCDs.  Each  fiberoptic  cable  consists  of  many 
strands  of  0.004  mm  diameter,  giving  a  total  cross  section  of  2.8  by  3.2  mm, 
defining  die  size  of  an  image.  There  are  about  20  strands  per  pixel,  and  in  the 
imager  systems,  diey  form  cables  diat  are  60  cm  long.  At  die  face  of  a  CCD, 
six  fiberoptic  cables  form  a  single  bundle  widi  dimensions  8.4  by  6.4  mm, 
wliich  covers  half  of  die  CCD.  The  image  conduit  bundle  allows  six  images  to 
be  collected  by  one  CCD.  In  die  spectrographs,  however,  the  fiberoptic 
cables  are  formed  into  short  rigid  blocks,  connecting  each  intensifier  to  a 
single  CCD  and  transferring  the  phosphor  unage  directly  to  die  CCD. 

Tlie  twelve  imagers  include  four  narrow  field-of-view  ultraviolet 
(NUV)  imagers  (FOV  1.8°  azimuth  x  1.6°  elevation),  four  wide  field-of-view 
ultraviolet  (WUV)  imagers  (FOV  25°  az  x  21°  el),  two  medium  field-of-view 
visible  (VIS)  imagers  (FOV  6.0°  az  x  5.3°  el),  one  wide  field-of-view 
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infrared  (WIR)  imager  (FOV  21°  az  x  19“  el),  and  one  narrow  field-of-view 
infrared  (NIR)  imager  (FOV  2.3°  az  x  2.0°  el). 

Tlie  eiglit  UV  imagers  have  narrower  bandpass  filters  than  the  four 
non-UV  imagers  in  order  to  provide  spectral  isolation  and  two-dimensional 
spatial  information  for  selected  UV  emission  features.  The  NUV  imagers 
must  use  focusing  mirrors  witli  relatively  long  focal  lengths,  requiring  larger 
diameters  than  desired  to  maintain  an  adequate  f-ratio.  The  large  diameter 
conflicts  with  tlie  design  requirement  to  keep  the  four  optical  axes  close 
enougli  to  use  a  single  intensifier.  This  problem  is  overcome  by  cutting  a 
single  mirror  into  equal  quarters  wliich  are  then  sliglitly  separated  to  provide 
four  distinct  optical  axes.  A  commercially  available  Cassegrain  mirror  system 
with  a  UV  coating  was  used  for  this  purpose.  Primary  and  secondary  mirrors 
were  radially  cut  into  four  pieces  and  remounted  with  9  mm  gaps  between 
their  straiglit  edges.  The  resulting  optical  axes  are  thus  mutually  separated  by 
9  mm,  yielding  foiu*  images  on  a  single  intensifier,  as  desired.  The  original 
Cassegrain  system  with  an  f-ratio  of  1.1  and  a  single  large  aperture  tlius 
becomes  a  system  with  an  f-ratio  of  2.2  and  four  smaller  apertures.  The  four 
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WUV  imagers  use  hemispherical  lenses  and  flat  spacers  made  of  magnesium 
fluoride  (MgF2)  with  3.5  mm  diameter  apertures  and  f-ratios  of  nearly  2.  Tliey 
also  share  a  single  intensifier.  The  NUV  and  WUV  imagers  use  interference 
filters  with  spectral  bandpasses  centered  at  1600,  2000,  2350,  and  2600  A, 
each  witli  a  250  A  spectral  half-width  (Viereck  et  al.,  1990). 

Tlie  VIS,  WIR,  and  NIR  imagers  obtain  pointing  and  target  briglitness 
information.  The  two  VIS  imagers  share  an  image  intensifier.  With  tlieir 
approximate  5®  fields-of-view,  tliey  coUect  broad-spectrum  data  fi'om 
extended  sources  and  star  fields  for  post-fliglit  pointing  analysis.  Each  has  a 
9  mm  diameter  lens  with  a  35  mm  focal  length,  yielding  a  3.8  f-ratio.  They 
are  identical  except  for  their  interference  filters:  one  filter  has  a  bandpass 
centered  on  5500  A  wliile  the  other  is  centered  on  7000  A,  and  both  filters 
have  spectral  half  widths  of  2000  A.  Tlie  WIR  and  NIR  imagers  are  not 
intensified  and  are  therefore  much  less  sensitive  tlian  the  otlier  imagers;  they 
can  observe  extremely  briglit  sources  without  suffering  damage  to  their 
detectors.  Images  are  formed  directly  on  tlie  ends  of  the  fiberoptic  cables  for 
these  much  simpler  IR  imagers.  Tlie  WIR  imager,  witli  an  approximately  20® 
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field-of-view,  monitors  target  briglitness,  while  the  NIR  imager,  with  an 
approximately  2®  field-of-view,  provides  fine-pointing  information. 
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APPENDIX  B:  CALIBRATION  OF  THE  GLO  INSTRUMENT 

Calibration  of  GLO  occurs  before,  during,  and  after  the  mission.  Pre- 
and  post-fliglit  calibration  consists  of  measuring  the  coalignment,  fields-of- 
view,  and  sensitivity  as  a  function  of  wavelength  of  both  the  imagers  and  the 
spectrographs.  Wavelength  calibration  and  spectral  line  shape  (instrument 
function)  determinations  are  also  accomplished.  Faulty  (“hot”)  pixels  and 
average  backgroimd  signal  for  specific  CCDs  are  identified  and  their  effects 
subtracted  during  post-flight  data  cahbration.  Observing  stars  of  known 
briglitnesses  in  specific  spectral  regions  provides  in-flight  and  post-fliglit 
calibration  and  die  best  check  of  instrumental  pointing  precision  and  spectral 
response. 

Laboratory  intensity  calibrations  are  referenced  to  a  calibrated 
photodiode  traceable  to  tlie  National  Institute  of  Standards  and  Technology 
(NIST)  in  Boulder,  Colorado.  The  spectral  luminosity  Bx  of  a  liglit  source  is 
detennined  using  diis  photodiode,  a  series  of  bandpass  filters,  and  optical 
elements.  Convolution  of  tiie  photodiode  response  with  tlie  filter  transmission 
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function  gives  the  spectral  luminosity  at  tlie  center  wavelength  of  the  filter 
bandpass:  one  filter  yields  a  single  value  for  luminosity  at  a  single 
wavelength.  (The  spectral  luminosity  is  assumed  constant  over  die  bandpass 
range,  typically  less  tiian  a  few  hundred  angstroms).  By  choosing  many 
narrow  overlapping  bandpass  filters,  a  spectral  luminosity  curve  for  the  liglit 
source,  measured  in  rayleiglis  per  angstrom,  can  be  constructed  over  any 

wavelength  range  of  interest.  Recall  that  one  rayleigli  (IR)  equals  10^  photons 

-1  -2 
sec  cm  . 

The  response  of  an  optical  detector  is  given  by  die  photometric 
equation: 

P  (in  photoevents  sec  *)  =  Ad  Q  T  e, 

-1  -2  .  -1 

where  Bx  =  source  luminosity  in  photons  sec  cm  steradian  at  wavelength 

2 

A.,  Ad  =  detector  area  in  cm  ,  Q  =  solid  angle  of  the  source  subtended  at  the 
detector  in  steradians,  T  =  optical  transmission  efficiency  (a  dimensionless 
ratio  of  die  nmnber  of  photons  transmitted  tlirougli  die  detector  divided  by  the 
number  of  photons  incident  on  the  detector),  and  8  =  detector  quantum 
efficiency  in  photoevents  photon*.  A  photoevent  registered  by  die  detector 


152 


results  in  an  electronic  signal  wliicli  can  be  quantified  by  an  analog-to-digital 
converter  (ADC).  The  ADC  produces  one  digital  Data  Number  (DN)  for  a 
specific  input  signal  fi'om  the  detector.  The  response  P  can  be  converted  to 
digital  response  R  by  tlie  formula  R  (in  DN/sec)  =  P  (in  photoevents/sec)  x  S, 
where  the  conversion  factor  S  depends  on  the  type  of  detector.  For  an 
unintensified  CCD,  Seed  is  approximately  0.02,  where  approximately  fifty 
photoelectrons  yield  a  signal  of  one  DN,  wliile  for  an  intensified  CCD,  Siced  is 
approximately  2.0,  where  two  DN  are  typically  registered  at  maximum  gain 
for  each  photoevent. 

The  solid  angle  Q  of  an  imaging  system  can  be  written  in  terms  of  tlie 

f-ratio  of  die  optics,  where  f-ratio  =  focal  length  divided  by  aperture;  Q.  =  n/4 

-2 

X  (f-ratio)  .  Using  this  definition  of  Q,  die  fact  diat  detector  pixels  each 
measure  22  mm  x  22  mm,  and  that  die  normalized  spectral  luminosity  equals 
imager  bandpass  in  angstroms  times  1  R/A,  we  obtain,  for  a  source  of 
briglitness  N  (R/A),  a  normalized  digital  response  per  pixel  given  by  the 
fonnula; 

R  (DN/pixel)  =  ^normalized  (DN/sec)  x  N  X  exposure  (sec). 
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Calibration  for  wavelengths  longer  tlian  4000  A  is  reproducible  to 
within  ±20%,  implying  an  overall  precision  of  about  80%  for  the  GLO 
spectrographs.  This  level  of  precision  is  representative  of  modem  optical 
instmmentation  standards.  No  degradation  of  response  has  been  detected  as  a 
result  of  detector  exposure  to  the  space  environment. 
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APPENDIX  C:  LIMB  TRACKING  VALIDATION 

Tlie  GLO  limb  tracker  follows  the  brightest  region  of  a  strip  through 
tlie  imager  FOV  by  moving  the  scan  platform  in  AZ  and  EL,  tliereby 
removing  die  shuttle  limit  cycle  motion.  Post-flight,  the  airglow  layer  heiglit  is 
calculated  by  combining  the  AZ/EL  coordinates  (corrected  to  shuttle 
coordinates)  of  die  GLO  scan  platform  and  the  shuttle  attitude  (given  by  the 
pitch,  roll,  and  yaw  axes). 

Position  data  is  available  describing  bodi  the  shuttle  attitude  and  GLO 
scan  platfonn  tracking  motion.  A  data  analysis  software  package  known  as 
SYBIL,  developed  by  Kalynnda  Berens  at  LPL  West,  allows  die  automatic 
adjustment  and  interpolation  of  this  position  data  to  automatically  extract 
spectrograph  and  imager  photometric  data  at  constant  altitudes.  In  diis 
section,  diat  process  is  validated. 

GLO  limb  tracking  capabihty  was  used  througliout  diis  data  set.  GLO 
pointing  control  adjusted  die  elevation  and  azimuth  of  die  GLO  boresight  in 
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order  to  track  the  airglow  layer.  Shuttle  pitch  information  and  actual  GLO 
pointing  direction  were  combined,  tlien  die  result  used  to  calculate  the  LOS 
elevation  angle  (LOS  0,  measured  up  from  nadir),  LOS  0  leads  to  tlie  tangent 
point  altitude  and  range  by  simple  riglit  triangle  relations.  The  calculated 
tangent  point  altitude  is  tlius  a  result  of  known  shuttle  pitch  variations  and 
unknown  vertical  motion  of  the  airglow  layer. 

Our  problem  is  to  separate  the  two  effects,  removing  the  influence  of 
shuttle  pitch  changes  to  discover  the  true  change  in  altitude  of  the  airglow 
layer. 


If  tliere  is  no  timing  error  between  the  GLO  clock  and  the  orbiter 
clock,  tlien  the  tangent  point  height  calculation  is  assumed  to  be  accurate.  We 
can  calculate  the  GLO  AZ/EL  coordinates  for  a  constant  tangent  point 
altitude  by  using  die  shuttle  attitude  record  provided  by  NASA.  The 
difference  between  diese  AZ/EL  coordinates  and  actual  AZ/EL  coordinates 
recorded  while  tracking  would  give  the  variation  in  die  altitude  of  the  airglow 
layer.  Consequently,  GLO  data  could  be  adjusted  to  true  altitude. 
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In  Chapter  2,  the  systematic  error  due  to  the  discrepancy  between  tlie 
GLO  clock  and  the  orbiter  clock  was  described.  The  correction  to  the  tangent 
point  altitude  is  calculated  by  SIBYL.  Tliis  calculation  uses  plane  geometry 
and  does  not  take  into  account  the  effect  of  small  angles  near  the  pole  of  GLO 
scan  platform  azimuth  axis.  Unfortunately,  the  niglit  sky  observation  geometry 
for  STS-69  required  data  to  be  collected  near  the  GLO  azimuth  pole.  This 
problem  has  been  referred  back  to  the  LPL  West  GLO  science  team.  For  tlie 
purposes  of  tliis  work,  the  airglow  layer  altitude  has  been  corrected  to  tlie 
sodimn  layer  at  90  km,  assumed  to  be  a  stable  reference. 

Figures  2.7  and  2.8  showed  tliat  the  maximum  excursion  in  shuttle 
pitch  was  typically  1.3°  during  a  niglit  pass.  Tliis  caused  the  LOS  to  sweep 
tlirougli  a  45  km  range  of  altitude  at  the  Earth’s  limb  at  a  range  of  2000  km. 
Tliis  would  have  severely  compromised  observation  time  on  the  limb  if 
tracking  had  not  been  used.  Althougli  GLO  scan  platform  pointing  accuracy  is 
known,  tliis  pointing  must  be  calibrated  in  terms  of  the  shuttle  coordinate 
system  to  acliieve  optimum  tracking. 
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In  order  to  correct  die  airglow  observational  altitude  to  the  absolute 
tangent  point  altitude,  tlie  scan  platform  was  calibrated  with  respect  to  die 
shuttle  coordinate  system  as  described  previously.  The  recorded  scan 
platform  motion  in  AZ  and  EL  must  now  also  be  reconciled  against  the  limit 
cycle  motion  of  the  shutde. 

The  shuttle  pitch  motion  causes  die  primary  change  in  GLO  pointing. 
Since  the  slit  is  at  an  angle  to  die  limb,  both  AZ  and  EL  are  affected.  Tracking 
compensated  for  pitch  changes  by  rotating  AZ  and  EL  in  the  opposite  sense 
of  those  changes.  Shuttle  pitch  was  recorded  in  real  time  by  NASA  ground 
control  at  one  second  intervals.  Thus  the  shuttle  pitch  changes  AP  (Figure 
C.l)  are  known  widi  high  temporal  resolution.  The  imposed  “sawtoodi 
function”  ripples  in  Figure  C.l  are  unresolved  as  yet,  but  are  apparendy  the 
result  of  data  truncation  in  die  NASA  recording  process. 
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Attained  AZ  and  EL  positions  (Figure  C.2)  were  also  updated  at  two 
second  intervals  in  response  to  the  error  determined  by  the  tracker  telescope 
on  GLO.  Attained  and  commanded  AZ/EL  positions  dififered  by  ±0.005° 
because  of  the  limits  of  scan  platform  pointing  accuracy. 

Spherical  geometry  gives  tlie  angular  offset  caused  by  AZ/EL  tracking 
(accounting  for  pitch  and  airglow  altitude  changes)  with  the  formula 

AF  ’  =  cos'^  [sin(ELi)sin(EL2)+cos(ELi)cos(EL2)cos(AZi-AZ2)], 

where  ELi  and  EL2  are  consecutive  elevation  values,  AZi  and  AZ2  are  die 
corresponding  azimutli  values,  and  AF  ’  is  die  compensating  angular  variation 
caused  by  tracking  (Figure  C.3).  The  sum  of  actual  tracking  variations  AF’ 
reflected  in  Figure  C.3  and  pitch  variations  AF  shown  in  Figure  C.l  should 
yield  a  small  residual  representing  die  vertical  motion  of  die  airglow  layer. 
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STS-69  AZIMUTH  AND  ELEVATION  FOR  ORBIT  72 


Figure  C.2.  STS-69  azimuth  and  elevation  while  tracking  during  orbit  72. 
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Figure  C.4  shows  the  tangent  point  altitude  calculated  from  LOS  0  as  a 
frmction  of  time.  Tracking  begins  at  4.4  minutes  into  tlie  orbit  and  ends  at 
34.2  minutes.  Figure  C.4  should  be  a  reflection  of  Figure  C.3,  and  yet  the  two 
plots  are  very  different.  LOS  0  is  a  secondary  product  derived  by  the  LPL 
West  science  team  from  NASA  mission  data.  The  discrepancy  between 
Figures  C.3  and  C.4  is  caused  by  the  need  for  spherical  geometry  calculations 
near  the  pole  of  tlie  GLO  scan  platform,  as  described  previously. 
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TANGENT  POINT  ALTITUDE  FOR  ORBIT  72 


MINUTES  AFTER  MET  4/13:45 


Figure  C.4.  Tangent  point  altitude  as  a  function  of  time.  This  plot  reflects  a  superposition 
of  shuttle  attitude  drift  and  GLO  scan  platform  motion. 
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Note  that  Figure  C.4  exhibits  tlie  sawtootli  variation  described 
previously.  The  sawtooth  has  a  period  of  25  seconds  and  an  amphtude  of 
0.05°.  It  represents  an  undetectable  change  in  line-of-siglit  direction  and  is 
not  an  artifact  of  tracking.  After  subtraction  of  AP the  new  plot  of  airglow 
altitude  variation  would  still  possess  the  sawtooth.  The  sawtooth  is  easily 
removed  by  smootliing  witli  a  running  point  average.  Averaging  over  50 
seconds  of  data  produced  the  best  results  (Figures  C.5  and  C.6).  The 
smoothed  plots  on  tlie  next  two  pages  were  expected  to  represent  tlie  relative 
variation  in  altitude  of  tlie  airglow  layer  as  GLO  followed  limb  brightness. 
However,  tliey  depict  remnants  of  tlie  shuttle  limit  cycle  and  not  true  airglow 
height  variations  because  of  the  limitations  of  the  SIBYL  analysis  package. 
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ALTITUDE  INFERRED  FROM  THE  DIFFERENCE 
BETWEEN  CALCULATED  TANGENT  HEIGHT 
AND  PITCH  COMPENSATION  DURING  ORBIT  71 


Figure  C.5.  Tangent  point  altitude  minus  pitch  drift  effects  for  orbit  71. 


ALTITUDE  INFERRED  FROM  THE  DIFFERENCE 
BETWEEN  CALCULATED  TANGENT  HEIGHT 
AND  PITCH  COMPENSATION  DURING  ORBIT  72 


Figure  C.6.  Tangent  point  altitude  minus  pitch  drift  effects  for  orbit  72. 
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APPENDIX  D:  RETRIEVAL  OF  VOLUME  EMISSION  RATES  FROM 
LOS  BRIGHTNESS  PROFILES 

Dascli  (1992)  described  the  technique  used  here  to  infer  vertical 
profiles  of  emitting  species.  His  method  is  an  algorithm  based  on  the  known 
properties  of  the  Abel  integral  (Hassani,  1991).  The  algorithm  assumes 
cylindrical  symmetry  in  the  plane  of  observation  and  involves  an  integration 
over  observed  data  P(rj)  to  obtain  values  of  the  originatmg  field  distribution 
F(ri).  Here,  the  P(rj)  are  the  line-of-sight  surface  brightness  intensities, 
integrated  along  tlie  path  fi'om  tlie  tangent  point  to  die  GLO  instrument, 
arising  fi’om  the  actual  vertical  profile  of  volume  emission  rates  F(ri).  The 
teclinique  is  appropriate  here  because  lateral  spherical  curvature  is  negligible 
in  our  fields  of  view,  making  cylindrical  geometry  highly  accurate. 

Tlie  integral  for  a  given  F(rj)  can  be  expressed  analytically  as  the  sum 
F(r,)  =  -il[I*(l)-I,(0)]P(r,-l)-2I,P(rj)+[I,(l)+I,(0)]P(r,+l), 

i 

where  ly  (0)  and  ly  (1)  are  symmetry  weighting  factors  and  Ar  is  tlie  distance 
increment  between  points  where  F(ri)  is  evaluated  (Dasch,  1992). 
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The  symmetry  weigliting  factors  ly  (0)  and  lij  (1)  are  defined  as  follows: 
lij(0)  = 

0  forj  =  i  =  0  orj  <i; 

forj=i^0; 

2>r  [  2j  J  •’ 

1  ,..ff(2i  +  l)^-4i^r+2j+l 
l[(2j-l)^-4i^f^+2j-l 


and 

= 

0 

forj  <  I; 

^[(2j+l)^-4i2f2-2jLj(0) 

2k 

for  j  =  i; 

^{[(2j+l)2-4i2]''2-[(2j-l)2-4i2]'«}-2jlij(0) 

for  j  >  i. 

2ji 


Dasch  confirmed  the  cylindrical  symmetry  of  the  observation  geometry 
is  completely  captured  by  the  weighting  factors.  Retrieval  of  the  original 


for  j>i; 
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geometry  is  independent  of  the  spatial  increment  Ar  between  tlie  observations 
P(rj),  as  long  as  oversampling  (i.e.,  redundant  observation)  is  avoided. 

Tlie  geometric  weighting  factors  lij  (0)  and  lij  (1)  need  only  be 
computed  once  and  can  be  applied  to  any  set  of  observations  witli  cylindrical 
symmetry.  The  calculated  values  tlien  fill  a  square  transformation  matrix  ly, 
whose  dimension  is  appropriate  to  the  nmnber  of  observations  performed. 
The  observations  P(rj)  fonn  a  column  vector;  premultiplying  this  column 
matrix  by  the  square  transformation  matrix  ly  yields  the  originating  column 
vector  F(ri)  describing  tlie  vertical  profile  we  seek. 

The  square  matrix  of  these  values  acts  as  tlie  operator  for  the  Dasch 
algoritlun,  teclinically  known  as  a  “three-point  Abel  inversion.”  Roble  and 
Hays  (1972)  described  a  related  Abel  integral  technique  used  to  process 
occultation  data  for  recovering  vertical  density  profiles  of  atmospheric 
constituents.  Others  have  discussed  similar  techniques  since  at  least  the  1960s 
(e.g.,  Barr,  1962).  Dasch  compared  several  similar  techniques  and  concluded 
the  tliree-point  Abel  inversion  suffers  tlie  least  from  noise  and  is  the  most 
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effective  in  deconvolving  sharp  data  discontinuities.  It  is  also  easy  to 
calculate  using  spreadsheet  software,  such  as  Lotus  123  or  Quattro  Pro,  For 
tliese  reasons,  the  tliree-point  Abel  inversion  was  adopted  for  this  research. 

Two  tests  were  performed  of  tliis  method’s  effectiveness.  Manually,  the 
integrated  brightnesses  along  line-of-siglit  chords  in  concentric  spherical 
shells  were  computed  as  altitude  was  decreased  along  a  radial.  These  LOS 
values  arose  from  a  constant  volume  emission  rate  of  2000  photons  cm’^  sec'^ 
in  a  10  km  tliick  uniform  layer  centered  at  100  km  altitude.  Applymg  the 
Dasch  algoritlun  using  a  Ar  of  1  km  to  tliis  LOS  input  profile  provided  a  fairly 
accurate  retrieval.  The  volume  emission  rate  profile  was  a  “box  function,” 
sharp-edged,  discontinuous,  and  entirely  nonphysical  (Figure  D.l).  The  LOS 
brightness  profile  resulting  fi’om  this  volume  emission  rate  profile  should  have 
been  very  difficult  for  Dasch’s  algorithm  to  process,  yet  tlie  metliod  provided 
output  very  close  to  tlie  original  profile  (Figure  D.2).  The  output  returned  a 
skewed  box-type  profile,  while  overestimating  the  original  value  by  about 
20%  at  the  bottom  of  the  layer  and  10%  at  tlie  top  of  the  layer.  Based  on  tins 
result,  tlie  metliod  was  deemed  very  effective.  To  test  tlie 
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INTEGRATED  LOS  PROFILE  USED  AS  INPUT 
BY  THE  DASCH  INVERSION,  MANUALLY  COMPUTED 
FOR  THE  ORIGINAL  "BOX  FUNCTION”  PROFILE 


INTEGRATED  LINE  OF  SIGHT  AMPLITUDE 


Figure  D.l.  The  integrated  line-of-sight  profile  for  the  box  function  computed  manually 
and  used  as  input  to  test  the  Dasch  algorithm. 
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COMPARISON  OF  ORIGINAL  "BOX  FUNCTION"  PROFILE 
TO  PROFILE  RETURNED  BY  THE  DASCH  INVERSION 
OPERATING  ON  THE  INTEGRATED  LOS  PROFILE 


Figure  D.2.  Comparison  of  the  original  box  function  profile  with  the  profile  retrieved  by 
the  Dasch  algorithm  using  the  previous  plot  as  input. 
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algorithm  fiirtlier,  more  realistic  input  profiles  were  required.  Tliese  were 
computed  using  a  FORTRAN  program. 

Dasch’s  metliod  was  tested  on  models  of  spherical  shell  layers  witli 
gaussian  vertical  profiles.  A  FORTRAN  program  (Appendix  E)  calculated  tlie 
resulting  vertical  line-of-siglit  profiles  associated  witli  these  layers,  given  a  Ar 
of  1  km.  Layers  were  modeled  by  a  simple  gaussian  fimction/  (.x)  =  c  exp 
[(a/2)  x^],  witli  layer  “thickness”  defined  as  the  full  widtli  at  half  maximmn 
(FWHM)  and  the  step  size  for  x  (representing  tlie  increment  Ar)  taken  to  be 
one  unit.  Here,  FWHM=2a  and  values  beyond  ±3a  were  set  to  zero.  Radial 
altitude  of  the  peak  was  varied,  and  resulted  in  different  line-of-siglit  path 
lengths  tlirough  tlie  layer  by  changing  tlie  curvature  of  the  layer.  The  a  values 
were  varied  to  change  layer  thickness  and  tlie  gaussians  were  multiplied  by 
constants  c  to  give  different  peak  values  for  the  various  tests.  Knowing  tlie 
original  gaussian  profile  and  tlie  integrated  LOS  profile  derived  from  it 
allowed  a  comparison  with  output  from  the  Dasch  algoritlim. 
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INTEGRATED  LOS  PROFILE  USED  AS  INPUT 
BY  THE  DASCH  INVERSION,  COMPUTED  IN  FORTRAN 
FOR  THE  ORIGINAL  GAUSSIAN  PROFILE 


INTEGRATED  LINE  OF  SIGHT  AMPLITUDE 


Figure  D.3.  The  integrated  line-of-sight  profile  computed  for  the  gaussian  profile  by  a 
FORTRAN  program  and  used  as  input  to  test  the  Dasch  algorithm. 
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COMPARISON  OF  ORIGINAL  GAUSSIAN  PROFILE  TO 
PROFILE  RETURNED  BY  THE  DASCH  INVERSION 
OPERATING  ON  THE  INTEGRATED  LOS  PROFILE 


Figure  D.4.  Comparison  of  the  original  gaussian  profile  with  the  profile  retrieved  by  the 
Dasch  algorithm  using  the  previous  plot  as  input. 
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Figures  D.3  and  D.4  show  the  results  of  this  test  for  a  gaussian  layer 
witli  c  =  5  km  and  c  =  1,  centered  at  an  altitude  of  100  km,  and  its  resulting 
integrated  line-of-siglit  profile.  The  test  parameters  to  be  retrieved  were  the 
altitude  of  peak  amplitude,  tlie  peak  amplitude,  and  tlie  layer  thickness 
(FWHM).  Performance  of  the  Dasch  algoritlim  is  excellent  between  ±2.5c, 
wliile  outside  these  limits,  tlie  metliod  produces  negative  values  which  are 
ignored.  The  retrieval  of  the  layer  geometry  is  excellent;  the  altitude  of  the 
amplitude  peak,  layer  thickness,  and  shape  of  the  vertical  gaussian  profile  is 
accurately  preserved.  For  profile  values  at  the  peak,  error  is  10%.  At  the  top 
of  the  layer  (+2.5ct),  amphtude  error  is  less  than  1%,  while  at  the  bottom  of 
the  layer  (-2.5ct),  errors  are  on  tlie  order  of  5%,  although  this  depends  on 
layer  thickness.  Error  increases  with  decreasing  thickness,  as  miglit  be 
expected  due  to  resolution  limits  caused  by  a  discrete  step  size.  The  Dasch 
algoritlim  appears  to  produce  accurate  estimates  of  volume  emission  rate 
profiles,  as  long  as  Ar  is  much  smaller  than  tlie  FWHM. 

Performing  tests  on  gaussians  witli  the  ratio  i?  =  [Ar  ^  FWHM]  ranging 
from  0.10  to  0.90,  two  of  tlie  key  parameters  proved  very  sensitive  to  R. 
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Retrieval  of  tlie  altitude  of  peak  amplitude  is  consistently  accurate,  but  not 
tliat  of  tlie  peak  amplitude  or  layer  thickness.  The  algoritlim  is  best  for  small 
R  (liigli  resolution)  and  effectiveness  deteriorates  with  increasing  R. 

A  value  of  R  =  0.10  gives  retrieved  amphtude  peak  values  that  are 
elevated  by  10%,  while  producing  minimal  errors  in  peak  altitude  or 
thickness.  Witli  R  =  0.25,  layer  tliickness  and  peak  altitude  are  retrieved  witli 
minimal  error,  but  peak  amphtude  is  increased  by  a  factor  of  3.  For  R  =  0.50, 
retrieved  layer  thickness  is  increased  by  1.25  and  peak  amphtude  by  4.8, 
wliile  for  R  =  0.90,  thickness  is  increased  by  the  factor  1.5  and  peak 
amplitude  by  5.8.  For  values  of  R  greater  than  1.0,  diis  method  is  best  used  to 
detennine  peak  altitude  only.  For  GLO  data  analyzed  here,  Ar  varies  from 
about  6  to  8  km.  If  tlie  airglow  layer  tliickness  is  assumed  to  be  about  10  km, 
tlien  R  is  about  0.70,  and  the  Dasch  algoritlim  could  broaden  tlie  retrieved 
layer  tliicknesses  by  up  to  35%  while  increasing  tlie  peak  emission  by  about  a 
factor  of  5.  Only  peak  emission  values  presented  here  have  been  adjusted. 
Potential  layer  broadening  of  3  to  5  km  is  within  the  current  uncertainty  in 
accepted  values  of  airglow  tliickness  and  is  not  systematically  corrected  here. 
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APPENDIX  E:  THE  GAUSSIAN  LAYER  ALGORITHM 


♦*♦**  program  to  calculate  INTEGRATED 
***♦*  LAYER  EMISSIONS 

*****  USING  GAUSSIAN  VERTICAL  DISTRIBUTIONS 
*****  FOR  EMITTERS 


*****  Paul  Bellaire,  LPL  West,  University  of  Arizona,  Apr  95  ***** 


DOUBLE  PRECISION  A(1000,I00) 

DOUBLE  PRECISION  C(1000,100) 

DOUBLE  PRECISION  COR(1000,100) 

DOUBLE  PRECISION  TOP(IOO) 

DOUBLE  PRECISION  H(IOOO) 

DOUBLE  PRECISION  P(IOO) 

DOUBLE  PRECISION  V(IOO) 

DOUBLE  PRECISION  E(IOOO) 

DOUBLE  PRECISION  z,W,SIGMA,BOT,QMID,T,THK,PI 

INPUT  LAYER  MEAN  ALTITUDE  AND  THICKNESS 

write  (*,♦)  (TfflS  PROGRAM  ASSUMES  A  GAUSSIAN 
DISTRIBUTION  FOR') 

write  {**)  ('AN  EMITTING  SPECIES;  YOU  MUST  KNOW  THE 
FULL-WIDTH’) 

write  (*,*)  ('HALF  MAXIMUM  (FWHM)  OF  THE  LAYER 
DISTRIBUTION!') 
write  {**)  ('  ') 

write  (*,*)  ('Input  layer  mean  altitude  in  km  (to  nearest  0.1  km  AGL):') 
read  C",*)  QMID 

write  (*,*)  ('Input  layer  tliickness  (to  nearest  0.1  km;  max  is  100);') 
write  (*,*)  (NOTE;  tliickness  =  FWHM  of  die  emitter  distribution!') 
read(*,*)THK 


QMID=QMID+6378D0 
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CALCULATE  TOP 'T',  BOTTOM  'BOT',  AND  RMS  FACTOR 
***  'SIGMA',  NOTING  THAT  OUR  GAUSSIAN  IS  GIVEN  BY; 

♦♦♦  F(X)  =  EXP(-[SIGMA/2]X**2) 

**♦  SIGMA  COMES  FROM  THE  FORMULA: 

♦**  In  0.5  =  -[SIGMA/2]  [X**2],  WHERE  FWHM  =  2X  =  THK 

NOTE:  (3  X  FWHM)  =  'W  ACCOUNTS  FOR  >99.7%  OF  LAYER 

SIGMA=1 .3863D0*(2D0/THK)**2 
T=QMrD+(l  .5D0*THK) 

BOT=QMID-(l  .5D0*THK) 

W=3D0*THK 

PI=ACOS(-1DO) 

OPEN  (55,FILE='a.out') 

INITIALIZE  ARRAYS  AND  START  DO-LOOP  FOR 
99  1-KM  DECREMENTS 


A(1,1)=T 

TOP(l)=T 

j=l 

DO  100  jd=l,99 


ij=j 

j=j+l 

TOP(j)=TOPaj)-lDO 

A(lJ)=TOPa) 

♦***  DIRECT  PROGRAM  FLOW  FOR  L.O.S.  PATHS 
WITHIN  OR  BELOW  LAYER 


IF  (A(l  J).LT.BOT)  THEN 


50 


GOTO  210 
ELSE 
GOTO  50 
END  IF 

CALCULATE  GAUSSIAN-WEIGHTED  L.O.S.  PATHS 
WITHIN  LAYER 


n=INT((ij)*10) 

i=l 

DO  200  id=l,n-l 


ii=i 

i=i+l 

A(ij)=A(iiJ)+0.1D0 

C(iJ)=2DO*(DSQRT(A(iJ)**2-TOP(j)**2)) 


IF  (i.EQ.2)  THEN 
COR(ij)=C(iJ) 

ELSE 

COR(iJ)=C(iJ)-C(iiJ) 

END  IF 

EVALUATE  EXPONENTIAL  AT  THE  MIDPOINT 
OF  THE  ALTITUDE  INCREMENT 

H(i)=((((A(iJ)+A(iiJ))/2D0)-QMID)**2)*(SIGMA/2D0) 

E(i)=DSQRT(SIGMA/(2D0*PI))*DEXP(-H(i)) 

Pa)=2*COR(iJ)*E(i)+Pa) 

V(j)=DEXP(-(SIGMA/2DO)*(TOP(j)-QMID)**2DO) 

IF  (i.EQ.2)  THEN 
WRITE  (55, 12)  PG) 

FORMAT  (lx,fl6.8) 

END  IF 
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200 


210 


He*** 


END  DO 
GOTO  100 

CALCULATE  GAUSSIAN-WEIGHTED  L.O.S.  PATHS 
BELOW  LAYER 


m=INT(W*10) 

A(lJ)=BOT 

k=l 

z=0D0 


DO  300  kd=l,m-l 


kk=k 

k=k+l 

z=z+lD0 

A(kJ)=BOT+(z*0.1D0) 

C(kj)=2DO*(DSQRT(A(kJ)**2-TOPG)**2)) 

IF  (k.EQ.2)  THEN 

COR(kJ)=C(kJ)-DSQRT(BOT**2-TOP(j)**2) 

ELSE 

COR(kJ)=C(kj)-C(kkj) 

END  IF 

EVALUATE  EXPONENTIAL  AT  THE  MIDPOINT 
OF  THE  ALTITUDE  INCREMENT 

H(k)=((((A(kJ)+A(kkJ))/2D0)-QMID)**2)*(SIGMA/2D0) 

E(k)=DSQRT(SIGMA/(2D0*PI))*DEXP(-H(k)) 

P(j)=2*COR(kJ)*E(k)+P(j) 
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IF  (k.EQ.2)  THEN 
WRITE  (55, 11)  PG) 
FORMAT  (lx,fl6.8) 


END  IF 


300  END  DO 

100  END  DO 


*sjc:jet* 


WRITE  TO  OUTPUT  FILE 


QMID=QMID-6378D0 

OPEN  (5,FILE='pgauss.out') 

WRITE  (5,9)  QMID,  THK 

9  FORMAT  (lx, 'Layer  altitude  =’,f8.2,4x,'Layer  tliickness=’,f8.2) 


WRITE  (5,10)  P 
10  FORMAT  (lx,fl6.8) 


OPEN  (8,FILE- vgauss.out') 
WRITE  (8,15)  V 
15  FORMAT  (lx,fl6.8) 


STOP 

END 
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APPENDIX  F;  THE  JOHNSTON-BROADFOOT  MODEL 

Following  this  introduction,  the  reader  will  find  tlie  code,  written  in 
“C,”  of  tlie  Jolinston-Broadfoot  model  (Jolinston  and  Broadfoot,  1993).  Tliis 
program  models  the  night  atmosphere  fi:om  60  to  120  km  altitude,  dividing  it 
into  30  finite  plane-parallel  layers,  each  of  2  km  deptli,  and  uses  a  time  step 
of  300  seconds.  Each  layer  is  homogeneous  with  a  constant  particle  number 
density  for  the  major  constituents  O2  and  N2,  based  on  the  assumption  of 
neghgible  diurnal  variation  of  these  constituents  at  die  altitudes  in  question 
(Rodrigo  et  al,  1991).  Output  consists  of  vertical  number  density  profiles, 
integrated  vertical  column  densities,  and  altitudes  of  peak  emission  for  the 
lines  and  bands  of  O2,  OH,  and  the  various  species  of  O  (Johnston,  1992). 

The  number  densities  of  the  minor  constituents  are  calculated  using  the 


=  -|j(Niwi)+ P(Ni)-UNi) , 


equation  of  continuity; 
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where  the  number  density  of  the  minor  constituent  is  Nj,  the  chemical 
production  and  loss  terms  for  the  minor  constituent  are  P(Nj)  and  L(Nj), 
respectively,  t  is  time,  z  is  altitude,  and  tlie  term  didz  (Nj  w,)  represents 
diffusion.  The  upward  diffiision  velocity  wj  of  a  minor  constituent  is  given  by 


Wi  =  -D(^ 


<?Ni 


,  Mig  ,  1  ^T\  y/m  d  rNi-]\ 
Ni  kT  Tdz^  '^N 


where  D  and  K  are  tlie  molecular  and  eddy  diffusion  coefficients, 
respectively,  T  is  temperature,  k  is  Boltzmann's  constant,  Mj  is  the  molecular 
or  atomic  mass  of  the  minor  constituent,  z  is  altitude,  g  is  the  acceleration  of 
gravity,  and  the  number  densities  of  the  major  and  minor  constituents  are  N 
andNj,  respectively  (Johnston,  1992). 


The  HOx  cycle  reactions  occur  sufficiently  swiftly  that  HOx  species 
concentrations  are  assumed  to  reach  equilibrium  “instantaneously”  witliin 
each  time  step.  The  reactions  of  the  HOx  cycle  depend  on  the  number 
densities  of  O,  O2,  O3,  and  N2,  but  the  reactions  do  not  significantly  alter  tlie 
number  densities  of  tliese  non-HOx  reactants  in  one  time  step.  Tlie  reaction 


185 


rates  in  tlie  HOx  cycle  can  then  be  assumed  to  be  proportional  only  to  the 
amount  of  HOx  species  present.  Given  the  concentrations  of  HOx,  O,  O2,  O3, 
and  N2,  the  equilibrium  concentrations  of  H,  OH,  and  HO2  are  calculated. 
Below  an  altitude  of  86  km,  atomic  oxygen  in  tlie  0(^P)  state  is  lost  during 
tlie  production  of  ozone.  This  loss  is  approximated  by  a  negative  exponential 
function  describing  O  concentrations  during  each  time  step.  Eleven  reactions 
for  production  of  O2,  O3,  and  OH  and  twenty-one  reactions  involving 
quenching  or  radiative  losses  for  excited  states  of  O,  O2,  and  OH  are  included 
in  tlie  model  (Johnston  and  Broadfoot,  1993). 

The  standard  input  profiles  for  constituent  nmnber  densities  used  by  the 
model  are  based  on  the  1976  US  Standard  Atmosphere  and  the  Air  Force 
Geophysics  Laboratory’s  1985  atmospheric  model  (see  Appendix  G  for  tlie 
numerical  data  in  tabular  form).  Below  90  km,  this  input  profile  is 
representative  of  middle  latitudes,  wliile  above  90  km,  it  is  more  typical  of 
averaged  winter  conditions.  Thus,  we  expect  this  input  profile  to  be  more  of  a 
general  guide  tlian  a  precise  model  for  low  latitude  equinox  conditions. 
Reactant  concentrations  in  the  actual  atmosphere  will  be  different,  but  are 
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assumed  to  follow  the  altitude  dependencies  and  vertical  structure  indicated 
by  the  Johnston-Broadfoot  model  output. 

Figures  F.l  tlirough  F.3  provide  volume  emission  rate  and  briglitness 
profiles  for  the  airglow  emissions  of  interest  in  this  research,  namely,  the  O2 
Atmospheric  band,  the  oxygen  line  at  5577  A,  and  tlie  OH  Meinel  band. 
These  profiles  were  derived  fi'om  tlie  output  profiles  for  excited  state  emitter 
number  densities  produced  by  die  Jolinston-Broadfoot  model,  shown  in 
Figure  F.l.  Figure  F.2  indicates  the  related  volmne  emission  rate  profiles, 
wliile  Figure  F.3  shows  the  resulting  line-of-sight  (LOS)  integrated  brightness 
profiles  as  seen  fi-om  the  topside  (neglecting  tlie  optical  deptli  of  O2  (0,0) 
band).  Figure  F.3  was  developed  fi'om  ouQiut  of  the  Dasch  algoritlnn 
described  in  Appendix  D.  Known  transition  probabiUties  were  used  to 
compute  tlie  volmne  emission  rates  and  brightnesses  in  Rayleiglis  from  die 
output  number  densities.  The  probabilities  for  die  O2  (0,0)  and  (0,1)  bands 
(0.0793  sec**  and  0.00391  sec**,  respectively)  were  obtained  firom  Vallance 
Jones  (1974).  Turnbull  and  Lowe  (1989)  provided  the  diermally  averaged 
probability  of  3.48  sec**  for  die  OH  Meinel  (6,2)  band  at  200  K.  The 
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transition  probability  of  1.34  sec*’  used  for  the  01  (5577  A)  line  was  given  by 
Tohmatsu  (1990). 

Note  tliat  tlie  profiles  of  emitter  concentrations  and  volume  emission 
rates  can  be  closely  approximated  by  gaussians.  The  Johnston-Broadfoot 
model  predicts  full  widths  at  half  maximum  (FWHM)  for  O2  Atmospheric  and 
OH  Meinel  emissions  to  be  about  10  km,  wliile  for  tlie  5577  A  line,  it  predicts 
a  FWHM  of  about  7  km.  The  model  predicts  the  altitude  of  peak  emission  at 
approximately  94  km  for  tlie  O2  Atmospheric  band,  at  96  km  for  the  atomic 
oxygen  5577  A  line,  and  at  about  89  km  for  the  OH  Meinel  band  (see  Table 
3.3  for  a  comparison  of  the  model  with  actual  observations). 
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CONCENTRATION  PROFILES  PRODUCED  BY  THE 
JOHNSTON-BROADFOOT  MODEL 
FOR  EXCITED  STATE  EMITTERS 


Figure  F.l.  Profiles  of  02(b  %+\  0(‘So),  and  OH  (X  produced  by  the  Johnston- 

Broadfoot  model.  These  excited  states  emit  in  the  night  airglow. 


ALTITUDE  (km) 


VOLUME  EMISSION  RATE  PROFILES 
FOR  AIRGLOW  EMISSIONS 
BASED  ON  THE  JOHNSTON-BROADFOOT  MODEL 
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VOLUME  EMISSION  RATE  (photons  cm’^  sec  *) 


Figure  F.2.  Volume  emission  rate  profiles  for  O2  (0,0),  O2  (0,1),  01  (5577  A),  and  the  OH 
Meinel  (6,2)  band,  based  on  the  Johnston-Broadfoot  model  output  shown  in  Figure  F.l. 


ALTITUDE  (km) 
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INTEGRATED  LINE  OF  SIGHT  BRIGHTNESS  PROFILES 
FOR  AIRGLOW  EMISSIONS 
BASED  ON  THE  JOHNSTON-BROADFOOT  MODEL 


Figure  F.3.  Integrated  line-of-sight  (LOS)  brightness  profiles  for  the  emissions  of  O2  (0,0), 
O2  (0,1),  01  (5577  A),  and  the  OH  Meinel  (6,2)  band,  based  on  the  Johnston-Broadfoot 
model  output  shown  in  Figure  F.l. 
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/*  NIGHTGLOW.C 
/♦ 

/*  AN  ATMOSPHERIC  MODEL  FOR 

/*  CALCULATING  DENSITIES  OF  NIGHTGLOW  SPECIES 

/* 

/*  MODEL  BY  DR.  JEFF  JOHNSTON  (1992) 

/*  WITH  THE  GUIDANCE  OF  DR.  A.  LYLE  BROADFOOT 

/* 

/*  MODIFICATIONS  AND  UPDATES  TO  RATE  COEFFICIENTS 
/*  AND  INPUT/OUTPUT  FORMAT  BY  PAUL  BELLAIRE  (1995-96) 
I* 

/*  MODEL  ASSUMPTIONS: 

/* 

/*  1)  THE  DOWNWARD  FLUX  OF  O  IS  CONSTANT 

/*  THROUGHOUT  THE  NIGHT 

/♦ 

/*  2)  CONCENTRATION  OF  02  AND  N2  IS  CONSTANT  WITH 

/*  RESPECT  TO  TIME 

/* 

/*  3)  DIFFUSION  CAUSED  BY  TEMPERATURE  GRADIENTS  IS 

/*  NEGLIGIBLE 

I* 

/*  4)  QUENCHING  COEFFICIENTS  FOR  ALL  O*  STATES  ARE 

/*  THE  SAME 

/* 

/*  5)  THE  REVERSE  CHEMICAL  REACTIONS  ARE  NEGLIGIBLE 

/♦ 

/*  6)  THE  HOx  REACTIONS  REACH  LOCAL  EQUILIBRIUM 

/* 

/*  7)  0(3P)  DISAPPEARS  IMMEDIATELY  AFTER  SUNSET 

/*  BELOW  70  KM 

/* 

/* 

/* 

/*  PRODUCTION  REACTIONS:  REFERENCE: 

/* 
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/♦  1)  O  +  O  +  M  ->  02*  +  M 

/*  2)  O  +  02  +  M  ->  03  +  M 

/*  3)  O  +  03  ->  02  +  02 

/*  4)  H  +  03  ->  OH*  +  02 

/*  5)0  +  0H->02  +  H 

/*  6)  O  +  H02  ->  OH*  +  02 

/*  7)  OH  +  03  ->  02  +  H02 

/*  8)  H  +  02  +  M  ->  H02  +  M 

/*  9)  H02  +  03  ->  OH  +  02  +  02 

/*  10)OH  +  HO2->H2O  +  O2 
/*  ll)H  +  H02->H2  +  02 
/* 

/*  QUENCHING  REACTIONS: 
/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 


Rodrigo  et  al.  1985 
Rodrigo  et  al.  1991 
Rodrigo  et  al.  1985 
Rodrigo  et  al.  1991 
Rodrigo  et  al.  1991 
Rodrigo  et  al.  1991 
Chamberlain  &  Hunten  1987 
Rodrigo  et  al.  1991 
Rodrigo  et  al.  1991 
Rodrigo  et  al.  1991 
Rodrigo  et  al.  1991 

REFERENCE: 

Rodrigo  et  al.  1991 

Rodrigo  et  al.  1991 

Radzig  &  Smirnov  1985 

Slanger  1972 

A.V.  Jones  1974 

Krauss  &  Neumann  1975 

Radzig  &  Smirnov  1985 

Radzig  &  Smirnov  1985 

Bates  1988 

Bates  1988 

Bates  1988 

Bates  1989 

Bates  1988 

Bates  1988 

Bates  1989 

Bates  1989 

Bates  1988 

Bates  1988 

Bates  1988 

Bates  1989 

Bates  1988 

Bates  1988 


12)  0(1D)  +  02  ->  0(3P)  +  02 

13)  0(1D)  +  N2  ->  0(3P)  +  N2 

14)  0(lD)->0  +  hv 

15)  0(1  S)  +  02  ->  0(3P)  +  02 

16)  0(1S)  +  N2  ->  0(3P)  +  N2 

17)  0(lS)  +  0->0(3P)  +  0 

18)  0(1  S)  ->  0(1D)  +  hv 

19)  0(lS)->0(3P)  +  hv 

20)  02(A3)  +  02  ->  02(X3)  +  02 

21)  02(A3)  +  N2  ->  02(X3)  +  N2 

22)  02(A3)  +  O  ->  02(X3)  +  O 

23)  02(A3)  ->  02(X3)  +  hv 

24)  02(A'3)  +  02  ->  02(X3)  +  02 

25)  02(A’3)  +  N2  ->  02(X3)  +  N2 

26)  02(A'3)->02(X3)  +  hv 

27)  02(A'3)->02(al)  +  hv 

28)  02(cl)  +  02  ->  02(X3)  +  02 

29)  02(cl)  +  N2  ->  02(X3)  +  N2 

30)  02(cl)  +  0(3P)  ->  02(X3)  +  0(1S) 

31) 02(cl)->02(X3)  +  hv 

32)  02(bl)  +  N2  ->  02(X3)  +  N2 

33)  02(bl)  +  O  ->  02(X3)  +  O 


/*  34)  02(bl)  ->  02(X3)  +  hv  Bates  1988 

/*  35)  02(al)  +  02  ->  02  +  02  Bates  1988 

/*  36)  02(al)  +  O  ->  02  +  O  Bates  1988 

/*  37)  02(al)  ->  02  +  hv  Bates  1988 

#mclude  <stdio.h> 

#include  <matli.h> 

#mcliide  <string.h> 

#mclude  <stdlib.h> 

/♦  DEFINE  CONSTANTS  FOR  SCALE  HEIGHT  DETERMINATION 
#defme  AMUG_K  1 . 1 5e-5 

/*  (AMUG_K,  in  Kelvin/cm,  is  given  by  tlie  following  formula: 

/*  gravity  *  mass  per  particle  *  particles  per  mole  /  gas  constant) 

/*  DEFINE  ATOMS,  MOLECULES,  AND  INITIAL  CONCENTRATIONS 
(RODRIGO  et  al.) 

/*  LET  ALTITUDE  [i]  BE  (i+60) 
void  set_varO; 

double  T[61],  D[61],  K[61],  N2_conc[61],  02_X3_conc[61],  03_conc[61]; 
double  0_3P_conc[61] ,  H_conc[61],  H02_conc[61],  OH_conc[61]; 
double  0_lD_conc[61],  0_lS_conc[61] ,  02_al_conc[61], 
02_bl_conc[61]; 

double  02_cl_conc[61],  02_A3_conc[61],  02_A3_cli_conc[61], 
OH_oz_conc[61]; 

double  03_new[61],  02_X3_new[61],  0_3P_new[61],  H_new[61], 
H02_new[61]; 

double  0_lD_new[61],  0_lS_new[61],  02_al_new[61],  02_bl_new[61], 
02_cl_new[61]; 

double  02_A3_new[61],  02_A3_ch_new[61],  0H_new[61]; 

/*  INITIALIZE  COLUMN  DENSITIES 


194 


double  03_col  =  0.,O2_X3_col  =  0.,O_3P_col  =  0.,H_col  =  0.,HO2_col  = 

0.; 

double  0_lD_col  =  0.,O_lS_col  =  O.,02_al_col  =  O.,02_bl_col  = 
O.,02_cl_col  =  0.; 

double  02_A3_col  =  0.,O2_A3_ch_col  =  0.,OH_col  =  0.; 

char  garb[20]; 

double 

diffI7],alt,pres_time,delta_t,H_total,kay[4][4],kappa[3],0_const,max[7]; 
int  ij,k,iter,peak[7],lwii[3]; 

FILE  *outfp,*eiiiitj^,*timel|),*init^; 
mainQ 
{ 

set_var(max,7,0.); 

set_var(0_l  D_conc,6 1,0.);  set_var(0_l  S_conc,6 1,0.); 
set_var(02_al  _conc,6 1,0.); 

set_var(02_b l_conc,6 1 ,0.);  set_var(02_c l_conc,6 1 ,0.); 
set_var(02_A3_conc,6 1,0.); 

set_var(02_A3_ch_conc,6 1,0.);  set_var(OH_oz_conc,6 1,0.); 
set_var(H_new,6 1 ,0.);  set_var(OH_new,6 1 ,0.);  set_var(H02_new,6 1 ,0.); 
set_var(03_new,6 1 ,0 .);  set_var(02_X3_new,61 ,0.); 
set_var(0_3P_new,6 1,0.); 

set_var(0_l  D_new,6 1,0.);  set_var(0_l  S_new,6 1,0.); 
set_var(02_al  _new,6 1,0.); 

set_var(02_b  1  _new,6 1,0.);  set_var(02_c  1  _new,6 1,0.); 
set_var(02_A3_new,6 1,0.); 
set_var(02_A3_cli_new,6 1 ,0.); 

prmtf("Please  enter  the  filename  of  your  initial  conditions  profile  :\t"); 

scanf("%s",garb); 

printf("\n"); 

initfp  =  fopen(garb,"r"); 

fscanf(initfjp,"%s  %s  %s  %s  %s  %s  %s  %s  %s  %s  %s",garb,  garb,  garb,  garb, 
garb, 

garb,  garb,  garb,  garb,  garb,  garb); 
for(i=0;  i<=60;  i++){ 
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fscanf(imtl^,"%d  %lf  %lf  %lf  %lf  %lf  %lf  %lf  %lf  %lf  %lf 

&j,  &T[i],  &D[i],  &K[i],  &N2_conc[i],  &02_X3_conc[i], 
&0_3P_conc[i],  &03_conc[i],  &H_conc[i],  &OH_conc[i], 
&H02_conc[i]); 

} 

fclose(init^); 

initfp  =  fopen("mit_conc'’,"w"); 

fprintf(ini^;'%8s\t%8s\t%8s\t%8s\t%8s\t%8s\t%8s\n'', 

"Altitude”,"[0(3P)]  ",”[02]  ","[03]  ","[0H]  ", 

"[H]  ","  [H02]\n"); 
for(i=0;i<=60  ;i++)  { 

^rintf(imt1^,"%8d\t%8.2e\t%8.2e\t%8.2e\t%8.2e\t%8.2e\t%8.2e\n", 

i+60,0_3P_conc[i],02_X3_conc[i],03_conc[i], 

OH_conc[i],H_conc[i],H02_conc[i]); 

} 

fclose(imtft)); 

printf("Please  enter  the  maximum  number  of  time  steps  desired  (dt  =  5 
min);\t"); 

scanf("%d",&iter); 

/*  THE  TIME  INCREMENT  IS  5  MIN  =  300  SEC 

I*  MODEL  STARTS  AT  SUNSET  AND  ASSUMES  MIDNIGHT  6  HRS  = 
21600  SEC  LATER 

deltat  =  300.0; 
pres_time  =  0.0; 
printf("\n\n"); 

timefp  =  fopen("time_conc","w"); 

J5)rintf(timej^,"%8s\t%8s\t%8s\t%8s\t%8s\t%8s\t%8s\t%8s\t%8s\n", 

"Time 

(lirs)","[0(lS)]","[0(3P)]","[02(A3)]","[02(A'3)]","[02(cl)]","[02(bl)]","[ 

02(al)]","[0H]\n"); 


/*  START  TIME  LOOP 
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for(j=Oj<itery++){ 

pres_time  +=  deltat; 

if(pres_time  >=  21600.0){  printf("Midmght!\n"); } 

/*  START  ALTITUDE  LOOP 
for(i=0  ;i<=60  { 

/*  CALCULATE  MOLECULAR  AND  EDDY  DIFFUSION  EFFECTS 
/*  AND  APPLY  BOUNDARY  CONDITIONS  FOR  DIFFUSION 

if(i==0){ 

di£EIl]  =  0.; 
di£f[2]  =  0.; 

difiI3]  =  (((4.0*D[i+l]+K[i+l])/2.0e5)* 
(H_conc[i+l]-H_conc[i]+ 

(H_conc[i+l]/T[i+l])*(T[i+l]-T[i]))+H_conc[i+l]* 
(4.0*D[i+l]*(1.0*AMUG_K/T[i+l])+K[i+l]* 
(28.*AMUG_K/T[i+l]))-3.9e6)/2.0e5; 
diff[3]  =  0.; 
diffI4]  =  0.; 
diff[5]  =  0.; 

} 

if(i==60){ 

diff[l]  =  (1.4el2-(((D[i]+K[i])/2.0e5)* 
(0_3P_conc[i]-0_3P_conc[i- 1  ]+ 
(0_3P_conc[i]/T[i])*(T[i]-T[i-l]))+0_3P_conc[i]* 
(D[i]*(16.*AMUG_KA’[i])+K[i]* 
(28.*AMUG_K/T[i]))))/2.0e5; 
diffl2]  =  0.; 
diff[3]  =  0.; 
difiI4]  =  0.; 
dijff[5]  =  0.; 
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/*  0(3P)  STATE 

if(i>0  &&  i<60){ 

diffll]  =  (((D[i+l]+K[i+l])/2.0e5)*(O_3P_conc[i+l]- 
0_3P_conc[i]+(0_3P_conc[i+l]/T[i+l])* 
(T[i+l]-T[i]))+0_3P_conc[i+l]* 
(D[i+1]*(16*AMUG_K/T[i+1])+K[i+1]* 
(28.*AMUG_K/T[i+1]))- 
(((D[i]+K[i])/2.0e5)* 
(0_3P_conc[i]-0_3P_conc[i-l]+ 
(0_3P_conc[i]/T[i])»(T[i]-T[i-l]))+0_3P_conc[i]* 
(D[i]*(16.*AMUG_K/T[i])+K[i]* 
(28.*AMUG_K/T[i]))))/2.0e5; 


/*  03  MOLECULE 

diff[2]  =  (((D[i+l]+K[i+l])/2.0e5)* 
(03_conc[i+l]-03_conc[i]+ 

(03_conc[i+l  ]/T[i+l])* 

(T[i+l]-T[i]))+03_conc[i+l]* 

(D[i+1]*(48*AMUG_K/T[i+1])+K[i+1]* 

(28.*AMUG_K/T[i+1]))- 

(((D[i]+K[i])/2.0e5)* 

(03_conc[i]-03_conc[i-l  ]+ 

(03_conc[i]/T[i])*(T[i]-T[i-l]))+03_conc[i]* 

(D[i]*(48.*AMUG_K/T[i])+K[i]* 

(28.*AMUG_K/T[i]))))/2.0e5; 


/♦  HATOM 


diffI3]  =  (((4.0*D[iH-l]+K[i+l])/2.0e5)* 
(H_conc[i+l]-H_conc[i]+ 

(H_conc[i+l]/T[i+l])*(T[i+l]-T[i]))+H_conc[i+l]* 

(4.0*D[i+l]*(1.0*AMUG_K/T[i+l])+K[i+l]* 

(28.*AMUG_K/T[i+1]))- 

(((4.0*D[i]+K[i])/2.0e5)* 
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(H_conc[i]  -H_conc  [i- 1  ]+ 
(H_conc[i]/T[i])^(T[i]-T[i-l]))+H_conc[i]* 
(4  0*D[i]*(l .0*AMUG_K/T[i])+K[i]* 
(28.*AMUG_K/T[i]))))/2.0e5; 


/*  H02  MOLECULE 

diff[4]  =  (((D[i+l]+K[i+l])/2.0e5)* 

(H02_conc[i+ 1  ]-H02_conc[i]+ 

(H02_conc[i+ 1  ]/T[i+ 1  ])* 

(T[i+l]-T[i]))+H02_conc[i+l]* 

(D[i+1]*(33.*AMUG_K/T[i+1])+K[i+1]* 

(28.*AMUG_K/T[i+1]))- 

(((D[i]+K[i])/2.0e5)* 

(H02_conc[i]-H02_conc[i- 1  ]+ 
(H02_conc[i]/T[i])*(T[i]-T[i-l]))+H02_conc[i]* 
(D[i]*(33.*AMUG_K/T[i])+K[i]* 
(28.*AMUG_K/T[i]))))/2.0e5; 


/*  OH  MOLECULE 

diff[5]  =  (((D[i+l]+K[i+l])/2.0e5)* 

(OH_conc  [i+ 1  ]-OH_conc[i]+ 
(OH_conc[i+l]/T[i+l])*(T[i+l]-T[i]))+ 

OH_conc[i+l]*(D[i+l]*(17.*AMUG_K/T[i+l])+K[i+l]* 

(28.*AMUG_K/T[i+l])H((D[i]+K[i])/2.0e5)* 
(OH_conc[i]-OH_conc  [i- 1  ]+ 
(OH_conc[i]/T[i])*(T[i]-T[i-l]))+OH_conc[i]* 
(D[i]*(17.*AMUG_K/T[i])+K[i]* 
(28.^AMUG_K/T[i]))))/2.0e5; 


/*  CALCULATE  CHEMICAL  SOURCES/SINKS 
/*  ASSUME  HOx  CYCLE  EQUlLffiRATES  IMMEDIATELY 


I*  AND  CALCULATE  CONCENTRATIONS 

/*  (TfflS  IS  THE  EFFECT  OF  REACTIONS  4  TO  1 1) 

/*  SOURCE:  "COMPREHENSIVE  CHEMICAL  KINETICS" 

(BAMFORD),  P.32 

Htotal  =  H_conc[i]+OH_conc[i]+H02_conc[i]+ 
((diff[3]+diff[4]+diff[5])-(2.0*OH_conc[i]* 
H02_conc[i]*1.7e-ll*exp(416.0/T[i])+2.0*H_conc[i] 
H02_conc[i]  *6 .66e- 1 2))*delta_t; 

kay[l][2]  =  03_conc[i]*1.4e-10*exp(-470.0/T[i]); 

kay[l][3]  =  02_X3_conc[i]*(02_X3_conc[i]+N2_conc[i])* 
5.5e-32*pow(T[i]/300.0,-l  .6); 

kay[2][l]  =  O_3P_conc[i]*2.2e-ll*exp(117.0/T[i]); 

kay[2][3]  =  O3_conc[i]*1.6e-12*exp(-940.0/T[i]); 

kay[3][l]  =  0.; 

kay[3][2]  =  0_3P_conc[i]*3.0e-ll*exp(200.0/T[i])+ 
O3_conc[i]*1.4e-14*exp(-580.0/T[i]); 

kappa[l]  =  (kay[l][2]*kay[3][2]H-kay[l][2]*kay[3][l]+ 
kay[3][2]*kay[l][3])/(kay[2][l]*kay[3][2]+ 
kay[2][l]*kay[3][l]+kay[2][3]*kay[3][l]); 

kappa[2]  =  (kay[l][2]*kay[2][3]+kay[2][l]*kay[l][3]+ 
kay[2][3]*kay[l][3])/(kay[l][2]*kay[3][2]+ 
kay[l][2]*kay[3][l]+kay[3][2]*kay[l][3]); 

H_new[i]  =  H_total/(l+kappa[l]+kappa[l]*kappa[2]); 

OH_new[i]  =  H_new[i]*kappa[l]; 


/*  IN  THE  CALCULATIONS  THAT  FOLLOW, 

/*  CERTAIN  RATE  COEFFICIENTS  HAVE  BEEN  MODIFIED;  SEE 
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/♦  JOHNSTON  &  BROADFOOT  (J&B),  1993 

/*  AND  J&B  CORRECTION,  1994 

/*  CASCADE  COEFFICIENT  FROM  McDADE  &  LLEWELLYN,  1987 

OH_oz_coiic[i]  =  H_new[i]*kay[l][2]*0.32/(276.0+ 
(02_X3_conc[i]+N2_conc[i])*  1 .5e-l  1  *0.934+ 

1 .5e-10*O_3P_conc[i]); 

H02_new[i]  =  H_new[i]*kappa[l]*kappa[2]; 

/*  DETERMINE  OZONE  CONCENTRATION 

/*  SOURCES:  REACTION  2  LOSSES:  REACTIONS  3,4,7, 9 


if(i<26){ 

03_new[i]  =  delta_t*(diff[2]-((0_3P_conc[i]* 

8.0e-12*exp(-2060.0/T[i])+H_conc[i]*1.4e-10* 

exp(-470.0/T[i])+OH_conc[i]*  1 .6e-12* 

exp(-940.0/T[i])+H02_conc[i]*  1 .4e-14* 

exp(-580.0/T[i]))*03_conc[i]))+ 

O_3P_conc[i]*(l-exp(-6.0e-34* 

pow((T[i]/300.),-2.3)*02_X3_conc[i]* 

(N2_conc[i]+02_X3_conc[i])*delta_t)); 

}else{ 

03jiew[i]  =  delta_t*((diff[2]+0_3P_conc[i]* 

02_X3_conc[i]*(N2_conc[i]+  02_X3_conc[i])* 

6.0e-34*pow((T[i]/300.0),-2.3)H(0_3P_conc[i]* 

8.0e-12*exp(-2060.0/T[i])+H_conc[i]*1.4e-10* 

exp(-470.0/T[i])+OH_conc[i]*1.6e-12* 

exp(-940.0/T[i])+H02_conc[i]*1.4e-14* 

exp(-580.0/T[i]))*03_conc[i])); 

} 

/*  DETERMINE  0(3P)  CONCENTRATION 
/*  SOURCES:  NONE  CHEMICAL 
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/♦  LOSSES:  REACTIONS  1,2,3,5,6 

if(i<26){ 

0_const  = 

02_X3_conc[i]*(02_X3_conc[i]+N2_conc[i])* 

6 .0e-34*pow((T[i]/300 . 0),-2 .3)+ 
OH_conc[i]’''2.2e-l  l*exp(l  17.0/T[i])+ 
H02_conc[i]*3  .Oe-1 1  *exp(200.0/T[i]); 

0_3P_new[i]  =  delta_t*((diff[l])- 

(2.0*O_3P_conc[i]*O_3P_conc[i]*(N2_conc[i]+ 
O2_X3_conc[i])*9.4e-34*exp(484.0/T[i])+ 
0_3P_conc[i]  *O3_conc[i]*8.0e- 12* 
exp(-2060.0/T[i])))+ 

0_3P_conc[i]  *(exp(-0_const*delta_t)- 1 .); 

}else{ 

0_3P_new[i]  =  delta_t*((diffll])- 

(2.0*O_3P_conc[i]*O_3P_conc[i]*(N2_conc[i]+ 
O2_X3_conc[i])*9.4e-34*exp(484.0/T[i])+ 
0_3P_conc[i]  *02_X3_conc[i]  *(N2_conc[i]+ 
02_X3_conc[i])*6.0e-34*pow((T[i]/300.0),-2.3)+ 
0_3P_conc[i]  *O3_conc[i]*8.0e- 12* 
exp(-2060.0/T[i])+0_3P_conc[i]*OH_conc[i]* 
2.2e-l  l*exp(l  17.0/T[i])+O_3P_conc[i]* 
H02_conc[i]*3.0e-ll*exp(200.0/T[i]))); 


/*  DETERMINE  0(1S)  CONCENTRATION  (STEADY  STATE) 
/*  SOURCE:  REACTION  30  LOSSES:  REACTIONS  15-19 

0_lS_conc[i]  =  (02_cl_conc[i]*0_3P_conc[i]*3.7e-12)/ 
(02_X3_conc[i]  *4.0e- 1 2  *exp(-87 0 .0/T[i])+ 
N2_conc[i]  *5  .Oe-1 7+0_3Pconc[i]*2  .Oe-1 4+ 
1.28+0.078); 


/*  DETERMINE  0(1D)  CONCENTRATION  (STEADY  STATE) 
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/*  SOURCE:  REACTION  18  LOSSES:  REACTIONS  12-14 

0_lD_conc[i]  =  (0_lS_conc[i]*1.28)/(02_X3_conc[i]* 

3.2e-l  1  *exp(67.0/T[i])+N2_conc[i]*l  .8e-l  1  * 
exp(107.0/T[i])+9.1e-3); 

/♦  THE  FACTORS  0.07, 0.03,  0.04, 0.18,  AND  0.06  THAT  FOLLOW  ARE 
/*  THE  APPROPRIATE  RATIOS  FOR  THE  VARIOUS  STATES  OF  02 

/*  DETERMINE  02(al)  CONCENTRATION  (STEADY  STATE) 

I*  SOURCES:  REACTIONS  1,27,32,33  LOSSES:  REACTIONS  35-37 


02_al_conc[i]  = 

(0.07*0_3P_conc[i]*0_3P_conc[i]*(N2_conc[i]+ 

O2_X3_conc[i])*9.4e-34*exp(484.0/T[i])+ 
02_A3_ch_conc  [i]  *  5 .2e- 1  +02_b  l_conc[i]  * 
N2_conc[i]*2.2e-15+O2_bl_conc[i]*O_3P_conc[i]*8.0e- 
14)/ 

(02_X3_conc[i]*1.5e-18+0_3P_conc[i]*1.6e-16+1.9e- 


4); 


/*  DETERMINE  02(bl)  CONCENTRATION  (STEADY  STATE) 

/*  SOURCES:  REACTION  1  LOSSES:  lUEACTIONS  32-34 


02_bl_conc[i]  = 

(0.03*0_3P_conc[i]*0_3P_conc[i]*(N2_conc[i]+ 

O2_X3_conc[i])*9.4e-34*exp(484.0/T[i]))/ 

((N2_conc[i]*2.8e-15+O_3P_conc[i]*1.0e-13)+8.7e-2); 


/*  DETERMINE  02(cl)  CONCENTRATION  (STEADY  STATE) 

/*  SOURCES:  REACTION  1  LOSSES:  REACTIONS  28-31 
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02_cl_conc[i]  = 

(0 .04  *0_3P_conc  [i]  *0_3P_conc[i]  *(N2_conc[i]+ 

O2_X3_conc[i])*9.4e-34*exp(484.0/T[i]))/ 

((02_X3_conc[i]*9.3e-12+N2_conc[i]*2.5e-12)+ 

O_3P_conc[i]*2.7e-12+0.55); 

/*  DETERMINE  02(A'3)  CONCENTRATION  (STEADY  STATE) 

/*  SOURCES;  REACTION  1  LOSSES:  REACTIONS  24-27 

02_A3_ch_conc[i]  =  (0.18*O_3P_conc[i]’''O_3P_conc[i]* 

(N2_conc[i]+O2_X3_conc[i])*9.4e-34*exp(484.0/T[i]))/ 

((02_X3_conc[i]*6.7e-ll+N2_conc[i]* 

1.8e-ll)+8.8e-l+8.4e-l); 

/♦  DETERMINE  02(A3)  CONCENTRATION  (STEADY  STATE) 

I*  SOURCES:  REACTION  1  LOSSES:  REACTIONS  20-23 

02_A3_conc[i]  = 

(0 .06  *0_3P_conc  [i]  *0_3P_conc[i]  *(N2_conc[i]+ 

O2_X3_conc[i])*9.4e-34*exp(484.0/T[i]))/ 
((02_X3_conc[i]*l  .0e-10+  N2_conc[i]*2.5e-l  1)+ 
0_3P_conc[i]*l  .5e-l  1+10.7); 


/*  REINITALIZE  COLUMN  DENSITIES 


O_3P_col=0.; 

O_lS_col=0.; 

O2_al_col=0.; 

O2_cl_col=0.; 

O2_A3_col=0.; 

O3_col=0.; 

H_col=0.; 


O_lD_col=0.; 

O2_X3_col=0.; 

O2_bl_col=0.; 

O2_A3_ch_col=0 . ; 

OH_col=0.; 

HO2_col=0.; 


for(i=0  ;i<=60;i++)  { 
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/*  CALCULATE  COLUMN  DENSITIES 

0_3P_col  +=  O_3P_conc[i]*2.0e5; 

0_lD_col  +=  O_lD_conc[i]*2.0e5; 

0_lS_col  +=  O_lS_conc[i]*2.0e5; 

02_X3_col  +=  O2_X3_conc[i]*2.0e5; 

02_al_col  +=  O2_al_conc[i]*2.0e5; 

02_bl_col  +=  O2_bl_conc[i]*2.0e5; 

02_cl_col  +=  O2_cl_conc[i]*2.0e5; 

02_A3_cli_col  +=  O2_A3_ch_conc[i]*2.0e5; 
02_A3_col  +=  O2_A3_conc[i]*2.0e5; 

OH_col  +=  OH_oz_conc[i]*2.0e5; 

H02_col  +=  HO2_conc[i]*2.0e5; 

Hcol  +=  H_conc[i]’'‘2.0e5; 

03_col  +=  O3_conc[i]*2.0e5; 

/*  DETERMINE  CONCENTRATIONS  FOR  NEXT  TIME  STEP 

03_conc[i]  +=  03_new[i]; 

0_3P_conc[i]  +=  0_3P_new[i]; 

H_conc[i]  =  H_new[i]; 

OH_conc[i]  =  OH_new[i]; 

H02_conc[i]  =  H02_new[i]; 

if(0_3P_conc[i]  <=  1.0e-4){  0_3P_conc[i]  =  l.Oe-4; } 
if(03_conc[i]  <=  1.0e-4){  03_conc[i]  =  l.Oe-4; } 
if(H_conc[i]  <=  L0e-4){  H_conc[i]  =  l.Oe-4; } 
if(OH_conc[i]  <=  l.Oe-4) {  OH_conc[i]  =  l.Oe-4; } 
if(H02_conc[i]  <=  l.Oe-4) {  H02_conc[i]  =  l.Oe-4; } 


/*  PRINT  SPECIFIC  COLUMN  DENSITIES 
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ifG  =  0){ 

printfC'Prograin  completed!\n"); 

/*  printfC’Iteration  0(1S)  0(3P)  02[A3]  02[A’3]  02[cl] 

02[bl]  02[al]  OHV');*/ 

/*  printfC  %  3d  %  8.2e  %  8.2e  %  8.2e  %  8.2e  %  8.2e  %  8.2e  %  8.2e 
%  8.2e\n",*/ 

/*  j,0_l  S_col,0_3P_col,02_A3_col,02_A3_ch_col,02_c  l_col,*/ 

/♦  02_bl_col,02_al_col,OH_col);*/ 

1^rintf(timefp;'%84f\t%8.2e\t%8.2e\t%8.2e\t%8.2e\t%8.2e\t%8.2e\t% 

8.2e\t%8.2e\n", 

(float)0+l)/12.,O_lS_col,O_3P_col,O2_A3_col,O2_A3_ch_col,O2_c 

Icol, 

02_b  1  _col,02_al_col,OH_col); 

} 

out^  =  fopen("reac_conc","w"); 
emit^  =  fopen("emit_conc'V'w''); 

printf("\nThe  initial  concentration  altitude  profiles  are  in  the  file 
’init_conc'.\n"); 

prmtf("The  end-of-nm  emitter  altitude  profiles  are  in  the  file  ’emit  conc'.W); 
printf("The  end-of-nm  reactant  altitude  profiles  are  in  the  file  'reac_conc'.\n"); 
printf(”The  time-dependent  vertical  column  densities  and  end-of-run  altitudes 
of  pe^  emissions  are  in  tlie  file  'time_conc'.\n"); 
printf(”\n"); 

^rintf(outj5),\ 

"%8s\t%8s\t%8s\t%8s\t%8s\t%8s\t%8s\n", 

"Altitude",”[0(3P)]",”[02]  ",”[03]  ", 

"[OH]  ","[H]  ","  [H02]\n"); 

^rintf(emit^,\ 

"%8s\t%8s\t%8s\t%8s\t%8s\t%8s\t%8s\t%8s\n”,\ 

”Altitude",”[02(A3)]","[02(A’3)]", 

"[02(cl)]","[02(bl)]",”[02(al)]”,"  [0(1S)]”,”  [OH]  \n"); 
for(i=0;i<=60;i++)  { 

^rintf(out^,”%8d\t%8.2e\t%8.2e\t%8.2e\t%8.2e\t%8.2e\t%8.2e\n”. 
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i+60,0_3P_conc[i],02_X3_conc[i],03_conc[i], 

OH_conc[i],H_conc[i],H02_conc[i]); 

^mitf(eiiiit^;'%8d\t%8.2e\t%8.2e\t%8.2e\t%8.2e\t%8.2e\t%8.2e\t%8. 

2e\n", 

i+60,O2_A3_conc[i],O2_A3_ch_conc[i],O2_cl_conc[i], 

02_bl_conc[i],02_al_conc[i],0_lS_conc[i],OH_oz_conc[i]); 

if(O2_A3_conc[i]>inax[0])  { 
max[0]=O2_A3_conc[i]; 
peak[0]=i; 

} 

if(02_A3_ch_conc[i]>max[l  ])  { 
m2ix[  1  ]=02_A3_ch_conc[i] ; 
peak[l]=i; 

} 

if(02_al_conc[i]>max[6])  { 

niax[6]=02_al_conc[i] ; 
peak[6]=i; 

} 

if(02_c  l_conc[i]>max[2])  { 

max[2]=02_c  l_conc[i] ; 
peak[2]=i; 

} 

if(02_b  l_conc[i]>max[3])  { 
raax[3]=02_bl_conc[i]; 
peak[3]=i; 

} 

if(0_l  S_conc[i]>max[4])  { 
max[4]=0_l  S_conc[i] ; 
peak[4]=i; 

} 

if(OH_oz_conc[i]>inax[5])  { 

iTiax[5]=OH_oz_coiic  [i] ; 
peak[5]=i; 

} 

} 

^rintf(tiinefp,''\nAltitiides  of  peak  emission  (in  kins):\n”); 


fprintf(tiiTie^,''\t\tHerzbergI\t\t%d\n",peak[0]+60); 

^rintf(time^,''\t\tHerzbergII\t\t%d\n",peak[2]+60); 

^rintf(time^,"\t\tChamberlain\t\t%d\n",peak[l]+60); 

i^riiitf(time^,"\t\tOI  5577\t\t\t%d\n",pei[4]+60); 

^rintf(timefi),''\t\tAtmospheric\t\t%d'ii",peak[3]+60); 

]^rintf(timeJ^,"\t\tIR  Atmospheric\t\t%dW',peak[6]+60); 

^rintf(time^,''\t\tMemel\t\t\t%d\n\n'',peak[5]+60); 

fclose(oul4)); 

fclose(time^); 

fclose(emit^); 

} 

void  set_var(a,  max,  set_no) 
double  *a,  set_no; 
intmax; 

{ 

inti; 

for(i=0;  i<max;  i-H-){ 
a[i]  =  set_no; 

} 

} 


END  OF  JOHNSTON-BROADFOOT  MODEL 
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APPENDIX  G;  THE  ATMOSPHERIC  VERTICAL  PROFILE  INPUT  TO 
THE  JOHNSTON-BROADFOOT  MODEL 


Alt 

Temp 

Molec 

Eddy 

m 

[O2] 

[0] 

[O3] 

[H] 

[OH] 

[HO2] 

km 

K 

cm^/s 

cm^/s 

cm'^ 

cm'^ 

cm'^ 

cm'^ 

cm*^ 

cm*^ 

cm*^ 

60 

243. 

1.0e3 

2.2e5 

5.0el5 

1.4el5 

6.0e9 

6.0e9 

2.5e6 

1.0e7 

4.3e6 

61 

240. 

l.le3 

2.2e5 

4.4el5 

1.2el5 

5.8e9 

4.7e9 

2.7e6 

9.0e6 

3.9e6 

62 

238. 

1.3e3 

2.3e5 

3.9el5 

l.lel5 

5.7e9 

3.5e9 

3.0e6 

8.1e6 

3.4e6 

63 

236. 

1.5e3 

2.4e5 

3.5el5 

9.7el4 

5.5e9 

2.7e9 

3.5e6 

7.6e6 

3.2e6 

64 

234. 

1.7e3 

2.5e5 

3.0el5 

8.5el4 

5.3e9 

2.0e9 

4.0e6 

7.2e6 

3.1e6 

65 

231. 

1.9e3 

2.6e5 

2.7el5 

7.5el4 

5.1e9 

1.6e9 

5.4e6 

6.8e6 

3.0e6 

66 

229. 

2.2e3 

2.7e5 

2.4el5 

6.6el4 

4.9e9 

1.2e9 

6.7e6 

6.5e6 

2.9e6 

67 

227. 

2.5e3 

2.8e5 

2.1el5 

5.8el4 

4.8e9 

9.7e8 

8.3e6 

6.3e6 

2.8e6 

68 

225. 

2.8e3 

3.0e5 

1.8el5 

5.1el4 

4.8e9 

7.5e8 

1.0e7 

6.2e6 

2.7e6 

69 

222. 

3.2e3 

3.2e5 

1.6el5 

4.5el4 

4.9e9 

6.1e8 

1.4e7 

6.2e6 

2.7e6 

70 

220. 

3.6e3 

3.4e5 

1.4el5 

4.0el4 

4.9e9 

4.9e8 

1.7e7 

6.2e6 

2.7e6 

71 

219. 

4.1e3 

3.5e5 

1.3el5 

3.5el4 

5.3e9 

4.1e8 

2.2e7 

5.6e6 

2.4e6 

72 

218. 

4.6e3 

3.6e5 

l.lel5 

3.0el4 

5.7e9 

3.3e8 

2.7e7 

5.0e6 

2.2e6 

73 

218. 

5.3e3 

3.7e5 

9.6el4 

2.7el4 

7.3e9 

3.0e8 

3.3e7 

4.1e6 

1.7e6 

74 

218. 

6.0e3 

3.8e5 

8.2el4 

2.3el4 

9.0e9 

2.8e8 

3.9e7 

3.2e6 

1.3e6 

75 

217. 

6.9e3 

3.8e5 

7.1el4 

2.0el4 

l.lelO 

2.4e8 

4.1e7 

2.3e6 

9.5e5 

76 

216. 

7.7e3 

3.9e5 

6.1el4 

1.7el4 

1.2el0 

2.1e8 

4.3e7 

1.5e6 

5.9e5 

77 

214. 

9.0e3 

4.0e5 

5.3el4 

1.5el4 

1.4el0 

1.8e8 

4.2e7 

l.le6 

3.6e5 

78 

213. 

1.0e4 

4.0e5 

4.6el4 

1.3el4 

1.5el0 

1.5e8 

4.2e7 

6.6e5 

2.3e5 

79 

211. 

1.2e4 

4.1e5 

4.1el4 

l.lel4 

2.2el0 

1.7e8 

3.8e7 

4.3e5 

1.4e5 

80 

209. 

1.4e4 

4.1e5 

3.6el4 

1.0el4 

3.0el0 

1.9e8 

3.4e7 

2.0e5 

5.3e4 

81 

207. 

1.7e4 

4.2e5 

3.1el4 

8.7el3 

4.7el0 

2.1e8 

3.0e7 

l.le5 

3.1e4 

82 

206. 

1.9e4 

4.2e5 

2.6el4 

7.3el3 

6.5el0 

2.3e8 

2.5e7 

4.3e4 

8.7e3 

83 

204. 

2.3e4 

4.3e5 

2.2el4 

6.3el3 

9.2el0 

2.2e8 

2.1e7 

2.6e4 

5.0e3 

84 

202. 

2.7e4 

4.5e5 

1.9el4 

5.3el3 

1.2ell 

2.2e8 

1.7e7 

1.2e4 

1.3e3 

85 

200. 

3.2e4 

4.6e5 

1.6el4 

4.5el3 

1.7ell 

2.1e8 

1.5e7 

7.7e3 

8.0e2 

86 

199. 

3.7e4 

4.8e5 

1.4el4 

3.8el3 

2.2ell 

1.9e8 

1.3e7 

3.4e3 

3.0e2 

87 

197. 

4.4e4 

4.9e5 

1.2el4 

3.3el3 

2.8ell 

1.5e8 

1.5e7 

2.8e3 

2.1e2 

88 

195. 

5.2e4 

5.0e5 

1.0el4 

2.8el3 

3.5ell 

1.2e8 

1.7e7 

2.2e3 

1.3e2 

89 

193. 

6.2e4 

5.5e5 

8.5el3 

2.4el3 

3.8ell 

l.leS 

2.6e7 

2.2e3 

1.2e2 

90 

192. 

7.2e4 

6.0e5 

7.1el3 

2.0el3 

4.1ell 

9.0e7 

3.5e7 

2.3e3 

1.2e2 

91 

194. 

8.6e4 

6.2e5 

6.0el3 

1.7el3 

4.6ell 

7.2e7 

4.4e7 

1.9e3 

9.6el 

92 

196. 

1.0e5 

6.5e5 

5.0el3 

1.4el3 

5.2ell 

5.4e7 

5.2e7 

1.5e3 

7.1el 

93 

194. 

1.2e5 

6.7e5 

4.2el3 

1.2el3 

5.4ell 

4.1e7 

5.4e7 

l.le3 

4.9el 

94 

192. 

1.5e5 

7.0e5 

3.4el3 

9.4el2 

5.6ell 

2.8e7 

5.5e7 

7.6e2 

2.7el 

95 

193. 

1.9e5 

7.5e5 

2.8el3 

7.9el2 

6.0ell 

2.0e7 

5.2e7 

5.2e2 

1.9el 
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Alt 

Temp 

Molec 

Eddy 

[N2] 

[O2] 

[0] 

[O3] 

[H] 

[OH] 

[HO2] 

km 

K 

cmVs 

cm^/s 

cm*^ 

cm'^ 

cm'^ 

cm‘^ 

cm'^ 

cm‘^ 

cm‘^ 

96 

195. 

2.2e5 

8.0e5 

2.3el3 

6.4el2 

6.3ell 

1.3e7 

5.0e7 

2.9e2 

l.lel 

97 

197. 

2.7e5 

1.0e6 

1.9el3 

5.4el2 

6.1ell 

9.7e6 

4.5e7 

2.0e2 

7.5 

98 

198. 

3.2e5 

1.3e6 

1.6el3 

4.4el2 

5.9ell 

6.4e6 

4.0e7 

1.2e2 

4.0 

99 

202. 

3.9e5 

1.3e6 

1.3el3 

3.7el2 

5.6ell 

4.4e6 

3.6e7 

8.5el 

2.5 

100 

206. 

4.6e5 

1.4e6 

l.lel3 

3.0el2 

5.2ell 

2.4e6 

3.2e7 

5.0el 

1.1 

101 

209. 

5.7e5 

1.4e6 

8.5el2 

2.5el2 

4.9ell 

1.7e6 

2.8e7 

3.3el 

1.0 

102 

212. 

6.8e5 

1.4e6 

7.1el2 

2.0el2 

4.5ell 

1.0e6 

2.4e7 

1.7el 

0.9 

103 

216. 

8.4e5 

1.4e6 

5.8el2 

1.7el2 

3.9ell 

6.5e6 

2.1e7 

1.2el 

0.6 

104 

220. 

1.0e6 

1.4e6 

4.6el2 

1.3el2 

3.4ell 

3.0e5 

1.8e7 

7.0 

0.3 

105 

224. 

1.3e6 

1.4e6 

3.8el2 

l.lel2 

3.1ell 

1.9e5 

1.5e7 

5.0 

0.2 

106 

229. 

1.6e6 

1.4e6 

3.0el2 

8.3ell 

2.9ell 

9.0e4 

1.2e7 

2.8 

0.11 

107 

234. 

2.0e6 

1.3e6 

2.4el2 

6.9ell 

2.6ell 

5.8e4 

l.le7 

2.0 

0.07 

108 

240. 

2.5e6 

1.2e6 

1.9el2 

5.4ell 

2.3ell 

2.6e4 

9.3e6 

1.0 

3.8e-2 

109 

246. 

3.2e6 

1.4e6 

1.6el2 

4.4ell 

2.2ell 

1.7e4 

8.5e6 

0.7 

2.5e-2 

no 

252. 

4.0e6 

1.7e6 

1.4el2 

3.5ell 

2.2ell 

8.0e3 

7.7e6 

0.3 

1.2e-2 

111 

260. 

5.2e6 

1.9e6 

1.2el2 

2.9ell 

2.0ell 

5.2e3 

6.8e6 

0.16 

l.Oe-2 

112 

268. 

6.3e6 

2.1e6 

l.lel2 

2.4ell 

1.8ell 

2.4e3 

5.9e6 

l.Oe-2 

8.0e-3 

113 

275. 

8.1e6 

2.2e6 

9.7ell 

2.1ell 

1.6ell 

1.6e3 

5.1e6 

8.0e-3 

6.0e-3 

114 

283. 

1.0e7 

2.4e6 

8.5ell 

1.7ell 

l.Sell 

7.6e2 

4.4e6 

6.0e-3 

4.0e-3 

115 

290. 

l.le7 

1.5e6 

7.5ell 

1.4ell 

1.4ell 

5.5e2 

4.0e6 

3.7e-3 

2.0e-3 

116 

302. 

1.2e7 

6.0e5 

6.6ell 

1.2ell 

1.3ell 

3.4e2 

3.5e6 

1.5e-3 

l.Oe-4 

117 

311. 

1.4e7 

4.6e5 

5.8ell 

l.Oell 

1.2ell 

2.6e2 

3.3e6 

l.le-3 

7.3e-5 

118 

320. 

1.5e7 

3.2e5 

5.1ell 

8.0el0 

1.2ell 

1.7e2 

3.1e6 

7.0e-4 

4.8e-5 

119 

330. 

1.7e7 

2.4e5 

4.5ell 

6.7el0 

l.lell 

1.3e2 

2.8e6 

4.2e-4 

2.7e-5 

120 

341. 

1.8e7 

1.6e5 

4.0ell 

5.5el0 

l.lell 

9.0el 

2.5e6 

1.5e-4 

l.le-5 

Data  Sources  (from  Johnston.  1992): 

•Temperature;  1976  US  Standard  Atmosphere  (January)  for  60-90  km;  AFGL  (1985)  for 
92-120  km,  interpolated  such  that  temperatures  matched  at  90  km. 

•Molecular  DiSusion  Coefficient:  Rodrigo  et  al.  (1986) 

•Eddy  Dififiasion  Coefficient:  Rodrigo  et  al.  (1986) 

•[N2]  :  Chamberlain  and  Hunten  (1987),  logarithmic  interpolation. 

•[O2]  :  Constant  proportion  of  N2  number  density. 

•[O]  :  Rodrigo  et  al.  (1991) 

•[O3]  ;  Rodrigo  et  al.  (1991) 

•[H]  ;  Rodrigo  et  al.  (1991) 

•[OH]  ;  Rodrigo  et  al.  (1991) 

•[HO2]  :  Rodrigo  et  al.  (1991) 
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